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CHAPITRE PREMIER 
THÉORÈMES GÉNÉRAUX SUR LE POTENTIEL NEWTONIEN 



Objet du cours. — Nous allons étudier la figure d'équilibre-^ 
d'une masse fluide animée d'un mouvement de rotation : les seules 
forces auxquelles est soumis le système sont des forces intérieures 
dues h l'attraction newtonienne : deux points s'attirent en raison 
directe de leurs masses et en raison inverse du carré de leur 
distance. Nous allons d'abord rappeler quelques résultats connus 
sur le potentiel newtonien. 

Déûnition et propriétés du potentiel. — Un point de masse 
772, soumis a Tattraction de tout le reste du système, est soumis 
à une force dérivant d'un potentiel V, c'est-à-dire que les com- 
posantes, suivant les axes de coordonnées, de la force agissant sur 
ce point, sont les dérivées partielles d'une même fonction 

dV ôV dV 

ôjc oy oz 

Si les masses attirantes sont discontinues, le potentiel a pour 
expression 

en appelant 7\ la distance du point de masse m au point de 
masse nii. 

PoiNCARÉ. — FigTires d'équilibre. i 
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Si ces masses sont continues, le potentiel a pour expression 



V = 




p' étant la densité de l'élément de volume d":' , r étant sa distance 
au point /;/, et l'intégrale étant étendue à tout le volume attirant. 
On sait que l'intégrale a un sens. 

Si la masse attirante est répandue sur une surface, avec une 
densité superficielle ^' ^ le potentiel a pour expression 



Y= fji^ 



Enfin, si c'est un volume infiniment mince d'épaisseur e', alors 
et le potentiel est 



Pour un volume, en dehors des masses attirantes, V est une 
fonction continue, ainsi que ses dérivées premières, dans tout 
l'espace, et on sait que l'on a 



d^V d'Y ^V 
dx'' "^ dtf "*" d,: 



^v = ^-r+^.^ -X-. 



A l'intérieur, en un point de densité p, on a 
AV = — 4^p. 

Les dérivées secondes sont discontinues en une surface de sépa- 
ration de deux milieux différents. 

Pour une surface, V est continue dans tout l'espace, mais ses 
dérivées ne le sont pas sur la surface elle-même. 

Sur la normale en un point A d'une surface S (fig. i), prenons 
deux longueurs 

AB = rf/?c, kC = dni. 
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Si Va est le potentiel en A, le potentiel eu B et C est 
ôV 

dV 

On a crailleurs 

dV . ôV 



hif. d/?, 




[JL étant la densité de la masse attirante. 

Si Ton désigne par Z, m, n les cosinus directeurs de la normale, 
supposée dirigée vers l'extérieur, par 

/dV\ / dV \ /^_\ /_ÔV\ /^\ /dV\ 

Vd^// Ui//e' U>3// \^x)i' \ày/i' \àz/i' 

les valeurs des dérivées partielles à Tintérieur et à l'extérieur de 
la surface, on a 

d/ii \dx / i \ài/ / i \ dz / i 

Théorème de Gauss. — Gauss a démontré que, en dehors 
des masses attirantes, Y n'a ni maximum ni minimum, puisque 
l'on a AV= o. 

D'une façon générale, si AV est positif ou nul, la fonction V 
n'a pas de maximum; si AV est négatif ou nul, il n'y a pas de 
minimum. Dans le cas qui nous occupe, AV est nul ou négatif: V 
n'a pas de minimum et n'a pas de maximum en dehors des 
masses agissantes. 

Donc, si sur une surface fermée ne renfermant pas de masses 
attirantes, V est compris entre deux nombres g et A, il en 
résulte qu'en tout point situé à l'intérieur delà surface, on aura 

s<y<h. 
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Théorème de Green, > — Si U et V sont deux fonctions con- 
ti)iues à l'intérieur d'un volume T limité par une surface S, on a 
l'égalité 

(0 / u4^c/cT= / UAWt 

j. '" i^ 

7/^^ dU OV OU ^V \ , 
\dx dx "^ dtj dy "^ d= d: } ~" 



Dans celte égalité, on peut permuter U et V, et en retranchant 
les deux égalités obtenues, on a 



Cas particuliers de l'égalité (i) : faisons U = i. 

(3) Ji^^^^jAV^.. 

Faisons : 



U = V 

(4) / V-^rfcr= / VAVrfT 



)"H-(-^(V(i^)> 






Si V est une fonction satisfaisant h l'équation de Laplace, 
AV=o, 

(^) /v^-/[(-£)V(f )V(^)]*. 

le second membre étant positif, le premier l'est aussi. 

Si U et V sont des potentiels, sur une sphère de rayon suffi- 
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samment grand, U et V sont de Tordre de -^ , leurs dérivées 

premières de Tordre de -rr^ • On peut alors appliquer la for- 
mule (i) à Tespace extérieur a une surface S et intérieur à une 
sphère de rayon très grand. Si ce rayon est suffisamment grand, 
l'intégrale prise sur la surface de lasphère (p est négligeable : 

c'est-a-dire est de Tordre de -5- ; mais on a la relation 



/"^"^X" 



dV , 

- — d(j = o 



car les éléments de Tintégrale se détruisent deux a deux. 

Quant aux intégrales des seconds membres, leur somme donne 
des intégrales étendues a Tespace entier, et on a, en ajoutant les 
deux égalités obtenues : 



(6) 



/TTw/ ^ T/dV dU ÔV ÙU , dV dU\ - 



les intégrales étant étendues à tout Tespace. 
De même, la formule (2) donnera 



/ 



(UAV— VAU)rfT=o; 

et, en se rappelant que 

— AV = 47tp 
_AV,= 47tp,, 



on a 



/ (V.p~V?.)rfT=0 



(7) 

l'intégrale étant étendue à tout, l'espace. 



u I 
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Cette formule est analogue à celle que Ton trouve dans la 
théorie de Télectriclté, 



2^„,V.-2;«.V = o. 



Dans la formule (7), faisons 
on en conclut la formule 



/ (Vrfo-f. 



(8) / (Vrfp — frfV)^^ = o. 



Travail de raitraction. < — Supposons un système de masses 
m', m'\ ..., etc. attirées par un autre/»',, m'\, ...., etc. 

Soit V'j, V'j,..- le potentiel aux points w', //i",..., etc., du aux 
masses attirantes 772',, m'\, ..., etc. 

V'j V"..., le potentiel aux points m\^ m'\,,,^ dû aux points 



Les masses attirées se déplaçant, on a comme travail 
en appelant 8:c, Sy, S.?, le déplacement du point de masse m ; 

Si je pose 

j'aurai 

Si les masses attirantes se déplacent aussi, mais seules, le 
potentiel V, deviendra 
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et on a 



Si les deux déplacements ont ïieu simultanément, on a 

e=3H+o'II. 
Supposons que les masses attirées forment un volume T, on a 



H= / p\,d^, 



cette intégrale pouvant d'ailleurs s'étendre à Tespace entier, 
puisque p est nul en dehors de T. On a 

SII== I V,SprfT 



=|v.v 



étendue de même à Tespace entier. 

Si la masse attirante est la même que la masse attirée, ou 
a V = V., 



= 1 Wdpd-z, 



ce qui peut s'écrire en tenant compte de (8) 

SH 



= / pSVrf-î, 



ou eacore 



.=/. 



Posons maintenant 



5H= / r'^' ■ d-c. 



w=j4*, 
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rîntégrale étant toujours étendue à tout l'espace; West Ténergle 
du système. 
On a 

8W = £. 

On peut d'ailleurs écrire, puisque AV = — ^n^y 



^/- 



et, en tenant compte de (6), 

rintégrale étant toujours étendue à l'espace entier. 

Conditions d'équilibre. — Ceci étant rappelé, nous allons ccwn- 
mencerrétude du problème que nous nous sommes proposé. 

Considérons une masse fluide, isolée de toute influence exté- 
rieure, tournant autour d'un axe fixe que nous prendrons comme 
axe O^, d'un mouvement uniforme de vitesse angulaire lo; prenons 
les axes Ox et Oy, invariablement liés a la masse fluide. 

Soit/? la pression en un point (^, y, z); p ne dépend que de 
.r, y, s, la force s'exerçant sur l'élément de volume d'z a pour 
composantes 

Ox Oy oz 

Écrivons les conditions d'équilibre relatif, en remarquant que 
l'accélération de Coriolis est nulle. Si on appelle X, Y, Z les 
composantes suivant les axes de coordonnées de la force agis- 
sant sur la molécule (.r, y, z)y on a 

^ ô V 

/ X =^ p -^- 1- iù^ ox\ 

7 ^V 
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et, en posant 



on a 



U: 


=V4- 


~i^' + y% 




do; 


dU 















d'où Ton tire 

-^dx-\ — ^dtj -\--J—dz^=^^ \-rr — dx-\- -r — du-\-——dz\ , 

djj=pd\J, 

Cette relation montre que p est uniquement fonction de /; ; U est 
par suite, aussi uniquement fonction de p et, par suite, de p. 
On arrive alors au théorème suivant: * 

Théorème, — Quand une masse fluide^ animée d'un mouç>e^ 
ment de rotation autour' d'un axe fixe^ est en équilibre relatifs les 
surfaces de nii>eau sont les surfaces d'égale pression^ et aussi 
d'égale densité. 

A la surface de la masse, on a p = o. Donc U est une cons- 
tante. La surface libre est donc une surface de niveau. 

S'il n'y a pas de rotation, on a U = V. 

Dîtns le cas général on a AV = — 4^p. AV est donc simple- 
ment fonction de U. 

Re„,arquons que la force totale -^, '^,^ est en général 

OX OfJ 0.\f 

normale aux surfaces de niveau, U = C. 

Ces conditions d'équilibre sont nécessaires. 

Une condition nécessaire et suffisante d'équilibre est que le 
travail résultant d'un déplacement virtuel soit nul. Ce travail 
comprend le travail de l'attraction, plus le travail di\ à la force 
centrifuge. 
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Le premier a pour expression oW, que nous avons déjà 
étudié, le second a pour expression 



/ 






2 



'■^ ?(•"'' 



■tf)d'Z 



et, si J est le moment d'inertie par rapport à Oz, le travail du à 
la force centrifuge est 

J. 

2 

La condition d'équilibre est donc 

S\V+i:^oJ==o 

2 

pour tout déplacement compatible avec les liaisons. 
Cette condition n'implique pas nécessairement que 

W + ^J 

2 

soit maximum ; il faut encore d autres conditions ; mais je dis 
que si cette fonction est maximum, l'équilibre est stable. 
Soit T la force vive relative, on a l'égalité 



T — (W + -Î^J)=C^ 



Posons W + — J =U, et supposons que pour le maximum 

U soit égale à U^. 
On a 

T-U=T,-U,. 



RELATION ENTRE LA MASSE, LE VOLUME ET LA VITESSE n 

Si on donne maintenant à U une valeur U<Uq et à T^, une 
valeur a, on aura 

T-U=T.-U. 

T = «-(U,-U). 

La force vive dans les instants suivants sera inférieure à Ty, 
Téquilibre est donc stable. 

Relation entre la masse, le volume et la vitesse de rotation. 
— Une condition nécessaire de l'équilibre est, qu'en tout point 
de la surface libre, la force soit dirigée vers l'intérieur, autrement 

AT " 

une partie se détacherait. 11 faut donc — — <o et pa* suite 



/ 



dU ^ 

■r— rf(J<0. 
0//^ 



La formule (3) transforme cette inégalité en la suivante : 

MJd-zKo, 



/■ 



Donc il faut que l'on ait 



/ 



AVtfT+ I 2 0>Vt<0 



c'est-à-dire 

r f 

2W^c/t<0. 



— 4-^1 ?di:-+- I 2 



En désignant le volume par T et la masse par M, il vient 
— 47tM + aw'T<o. 

Dans le cas particulier où la densité est constante et égale à p, 
on a 

2U)^<47îp. 
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Dans Tapplicalion de cette formule, il faut se rappeler que les 
unités sont telles que la force d'attraction newtonienne s'ex- 
prime par la formule 

mm' 



CHAPITRE II 
MASSE HOMOGÈNE FLUIDE 



LA MASSE EST SANS MOUVEMENT DE ROTATION 

Simplification des formules générales. — Etudions le cas 
particulier d'une masse homogène fluide et sans mouvement de 
rotation. 

En appliquant les théorèmes généraux, on trouve que la sur- 
face libre doit être une surface équipotentielle V = Vq. 

Si on prend la densité de la masse pour unité de densité, le 
potentiel en un point a pour expression 




L'énergie est 




V est continu dans tout l'espace ainsi qu'à la surface libre ; il est 
nul à l'infini ; a l'extérieur il n'a ni maximum ni minimum et est 
compris entre o et V^; h l'intérieur il est supérieur a V^ et a un 
maximum. 



L'intégrale 






d's étendue a une surface fermée entou- 
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rant le corps a pour valeur / AV^t, étendue au volume limité 

1-' 

par cette surface, et on a I A\Wt = — 4*^ I P^'^ = — 4'^^^». 

Ce théorème est encore vrai pour la surface limite, et il s'ap- 
plique que le corps soit creux ou plein, et même si l'on a affaire 
à une distribution superficielle. 

Comme enfin nous supposons le corps homogène et de den- 
sité I, on a 

/ -^^7= / AVrfT=— / 47ïprfT= — 47îT 

en désignant par T son volume, et par S une surface qui le 
renferme complètement, en particulier sa propre surface. 

Si on porte sur la normale en M a la surface extérieure un 
segment MM' = rf/ic(fig. 2), le potentiel en M'a pour expression : 

V V ^ ^V . 




Or, on a 



on en conclut que 



ÙV 

<o. 



d«e 



On a vu que cette condition était nécessaire pour l'équilibre : 
ici elle est toujours remplie. 

Soit V(j le potentiel de la masse fluide sur sa surface, et suppo- 
sons maintenant une couche électrique de masse M, en équilibre 
à la surface du corps, et produisant dans l'espace un potentiel 
V'tel que sur la surface et a son intérieur le potentiel soit V'=Vjj. 

A l'extérieur, on a V = V. En effet, sur la surface V' = V^ et 
h rinfini V'= o. Or V et V sont deux fonctions satisfaisant à 
l'équation de Laplace, ayant les mêmes valeurs sur la surface et 
sur une sphère de rayon très grand. Ces deux fonctions coïnci- 
dent donc : c'est le principe de Dirichlet. 
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On a, d'après des formules connues, 



/ 



On 



On a d'ailleurs 

■ dar= — 4'^T. 






Donc la masse M répandue à la surface a pour valeur T et réci- 
proquement si 

M==T, V' = V. 

T , 

Le rapport ,. est ce que Ton appelle la capacité électrique 

du corps. 

Cette notion va nous être utile pour résoudre le problème 
proposé. 

Théorie de Liapounoff. — Parmi les figures d'équilibre, il y 
a une solution évidente qui est la sphère : mais peut-il y avoir 
d'autres figures d'équilibre ? Nous n'en savons rien : mais si Ton 
se propose de chercher les figures d'équilibre stable, le pro- 
blème est différent. 

Le maximum absolu de W est bien atteint pour la sphère, 
mais a priori nous ne pouvons affirmer que W n'ait pas d'autres 
maximums relatifs. Cela a été démontré par LiapounofF; la 
démonstration est assez longue et on va procéder par étapes. 

L'énergie a un maximum, — Le potentiel d'un volume 
quelconque par rapport h un point quelconque, est inférieur au 
potentiel d'une sphère de même volume par rapport à son centre. 

Du point M, comme centre, je décris la sphère S ayant le même 
volume que le corps S (fig. 3). 

Le corps S est partagé en deux volumes : le volume K inté- 
rieur à S, et le volume A extérieur. 

Soit B le volume de S extérieur à S. 



^ 
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On a. par hypothèse B = A. 

Si n est le rayon de la sphère, on a 



Potentiel de A < — - 
Potentiel de B > — jt- = -jt- 



\ 


/ D^oil 
/ Potentiel de B > Potentiel 


*M 


' "^ / 


de A 


y 


Potentiel de S = Potentiel 


IL 


de K + Potentiel de A 


Fig. 3. 


Potentiel de S = Potentiel 




de K + Potentiel de B. 



Donc 

Potentiel de 2 > Potentiel de S. 
Si V est le potentiel au point M, et R le rayon de la sphère 2, 
on a donc 

V<2-R-. 
L'énergie qui est égale a 



est donc inférieure à 



7:R^ / rfT; 



et il vient 

W<7zR-T 

en appelant T le volume de la niasse. 
On a d'ailleurs 

T=-y7TR% 

d'après la définition de R. 
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Donc 

L'énergie W due à un volume T est donc lîmitée supérieure- 
ment ; W a donc un maximum. Je dis que ce maximum n'est 
atteint que si la surface extérieure est une sphère. 

Nous allons d'abord démontrer le théorème suivant : 

La capacité électrostatique a un minimum. — De tous les 
corps ayant même volume, la sphère est celui qui a la moindre 
capacité électrostatique. 

En effet, la formule (9) (p. 8) donne 



W: 



=^/[(-^)'+(^)'+(-ï)>. 



l'intégrale étant étendue à tout l'espace. 

Supposons que sur un corps conducteur de volume T il y ait 
une masse électrique T en équilibre. Si V est le potentiel corres- 
pondant, sur la surface du corps ou a 

CV = T 

par définition de la capacité C du conducteur. 

A l'intérieur du corps, Y' est constant; a l'extérieur, V satisfait 
a Téquation AV = o et est nul h l'infini. 

L'énergie électrique est, d'après une formule connue de l'élec- 
tricité : 

1 I T^ 
W' = — CV'^=— 4r- 

2 2 G 

Elle est d'ailleurs égale à 

w-^/[(^)V(f)V(^)>, 

l'intégrale étant étendue, soit à tout l'espace, soit au volume 
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extérieur seulement, puisque V est constant h Tintérieur du 
volume. 

Nous allons voir qu'elle est inférieure h Ténergie W due au 
potentiel newtonien du corps supposé homogène et de den- 
sité égale à l'unité. 

Soit V = V' + U ce potentiel : je pose 



W: 



^/i[^^^T+[^^]-[^^^]'i-. 



l'intégrale étant étendue à tout l'espace. 






La première intégrale est W. 

La deuxième se décompose en deux : l'une prise à l'intérieur 

. , „ . OV dY dT 

du corps, qui est nulle, puisque — — = - — = — — = o. 
* ^ ùj: Oi/ oz 

L'autre, prise à l'extérieur, se transforme par la formule de 
Green en 






La seconde de ces deux intégrales est nulle puisque AU est nul 
l\ Textérieur. 

La première est également nulle : en effet, sur la surface, V 
est constant et égal à V',,. 
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L'intégrale est donc 



=vi/ Jl.,._v; /-gl.,; 
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or 




rfo- = 47^ 1 , 



^^'' 7 / T 



cette intégrale est donc nulle. 
On a par suite 




^^ ^u 5V au ÔV dU \ _ 

ôx ôy ^y d.3 0.3/ 



La troisième n'est pas nulle, car U n'est pas nul à l'intérieur 
du corps, puisque V y est constant, et que V y varie. 
La troisième intégrale est alors positive, et on a bien 

w>w. 

On a donc, en réunissant les deux résultats précédents : 
7rR^T>W>-^, 



d'oii 



C> 



27tR^' 

et enfin 
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C a donc un minimum. Je dis que ce minimum ne peut être 
atteint que pour la sphère. 

Le minimum est atteint pour la sphère. — Supposons qu'un 
corps conducteur ait la forme pour laquelle ce minimum est 
atteint : on a 

M' 



'=^=^/[(^)+(^)'-(^)i- 



— — désignant la dérivée suivant la normale a la surface équl- 
potentielle passant par le point [Xy y, z)^ 

2W 



'-i/(^)"- 



Déformons le conducteur de façon que l'élément de surface 

d^ devienne d^ (fig. 4). Appelons 
ÎJ la projection sur la normale à la 
surface du déplacement d^ compté 
avec son signe. 

L'augmentation du volume est 
alors 




Fig. 4. 
L'augmentation de l'intégrale est 
WdC 




h- 



Ak^J '^W) -^[sr) J 



Ceci suppose que le potentiel reste constant; mais, par hypo- 
thèse, la figure est en équilibre. Donc il y a un minimum pour 
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M- 

— - ; la variation de Tîntégrale due a la variation du potentiel 

est donc nulle, puisque cette intégrale est minimum» 
Or, si on appelle [jl la densité superficielle, on a 

dV , 



On a donc 



d/i 







Or, il faut que rfC == o : donc 



/• 



et comme on doit avoir en même temps /!^rfcr=o, quel que 

soit IJ, il faut que Ton ait [x = C'^. 

On a d'ailleurs M = [xS en appelant S la surface du conduc- 
teur : on a donc 

Supposons maintenant que le conducteur se déforme en restant 
semblable à lui-même : la capacité est proportionnelle à la 
racine cubique du volume, et on a 

"C"""T T ' 
d'où on conclut 

S2=I27tTC. 

Tous les corps dont la capacité est minimum sont donc tels 
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que cette relation existe entre leur capacité^ leur volume et leur 
surface. 

Mais alors, pour tous les corps de même volume, la capacité est 
minimum si la surface est minimum ; or, parmi tous les corps 
de même volume, la sphère est celui qui a la surface minimum. 
C'est donc celui dont la capacité électrostatique est mini- 
mum. 

Nous arrivons maintenant au théorème que nous nous étions 
proposé de démontrer. 

La sphère est la seule ûgnre d'équilibre. — Je dis que W est 
maximum quand le corps de volume T a la forme d'une sphère. 

Soit T le volume correspondant, donnons-lui un accroisse- 
ment (fig. 4) 



Comme on a 

W 



■/"• 






dvf= h^i+^^d. 



\d0d7. 



Le corps étant supposé de densité égale a i , rfp est égal h — i , 
ou I ; dp n'est différent de zéro que dans le voisinage de la 
surface, qui est une surface équipotentielle, le corps étant sup- 
posé en équilibre ; c'est-à-dire que V est constant. 



rfw=v, ; dpd^ 



/• 
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et l'intégrale est égale à dT, On a donc 

On peut établir une autre relation entre dW et dT. 
Supposons que le corps varie en restant semblable ii lui-même : 

l'énergie varie en restant proportionnelle à la puissance -- du 

o 

volume. On a donc 

\\ ~ 6 T * 

Eliminant dVf et dT entre ces deux dernières relations, on 
trouve 

W = — VT. 
D'après la définition de W, on voit que Ton peut appeler 

2W 



le potentiel moyen. Il est égal h 

Or W est aussi donné par l'intégrale 

étendue à l'espace entier. 

La partie étendue à Fespace intérieur est, d'irprès la formule 
de Green (4) : 
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Sur la surface, V est constant. L'intégrale a donc pour 
valeur 

V / AV w / 



if^^-wj^ 



ce que l'on transfornie immédiatement, puisque — AV =^ 4^p 
= 4^, en 



èf"'-if 



VAVrfx 



+ W 

La partie de l'intégrale étendue au volume intérieur est donc 
h la partie étendue au volume extérieur, dans le rapport cons- 
tant -p, en supposant évidemment le corps en équilibre. 

Mais l'énergie de ce corps est en raison inverse de sa capa- 
cité électrostatique. 

Nous voyons donc bien que la sphère est le volume pour lequel 
W est maximum. 



LA MASSE EST ANIMEE D UN MOUVEMENT DE ROTATION. 

Formules générales. — Supposons maintenant le corps, 
homogène, soumis à une rotation w autour d'un axe fixe. 
Appelons J le moment d'inertie autour de cet axe. 



= 1 d':{x' 



+ y^). 
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Si le corps change de forme, de façon que la projection du 
déplacement de l'élément de surface d7 sur la normale soit !J, 
J varie et on a 



■/ 



rfj= / ^d^Çv^ + tf). 



Nous avons appelé U la fonction VH (•^^ + y^)« 

W représente toujours l'énergie due au potentiel new- 
tonien. 
On a 

dW+^dJ= / !:[v+-î^(.r^ + y^)]rfT 

A la surface d'équilibre U est constant et égal à U^. 
On a donc 



^dJ=V, / ^d^ = V,a 



D'ailleurs si le corps se déforme en restant semblable à lui-même, 

5 
W et J varient comme la puissance -rr- de T. 

On a donc 





rfW+ dJ 
1 


5 STi 




2 


~ 3 T 


et par suite 







4UoT = W + — J = U. 

D 2 



a6 


MASSE MOM 


On a 






AU = 2(o2 + AV, 




= 2W" 4^0, 




= 20>'^ 4^. 



On a déjà vu (page ii) qu'une condition nécessaire pour 
l'équilibre est que 

c'est-à-dire que AU soit négatif. 

Donc U, à rintérieur du corps, ne peut avoir de minimum, et 
comme sur la surface il est constant et égal à U^, à l'intérieur on 
aU>U.. 

Si AU est nul, U est constant et égal à U^. 

Enfin si AU était positif, U serait dans tout le corps inférieur 

Considérons maintenant l'intégrale 



/ 



UrfT=2W+— J. 



Suivant que -AU est négatif, nul ou positif, le premier membre 
est supérieur, égal ou inférieur à 



On a donc dans les trois cas : 
AU 



K / ** \ 2 

AU = o 4-(w + -^J) = 2W+— J W = w''J 

O \ 2 / 2 

AU > o -|- (w H- -^j)> aW+ -^ J W < w^J. 
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Limite de la vitesse de rotation. — Supposons maintenant 
que io varie d'une façon continue. La figure se déforme d'une 
manière continue, et on a 

2 5 

La figure étant supposée en équilibre, W -| J est maximum 

ou minimum, en tous cas 

rfW + — rfJ = o. 



Il reste 



Donc 



3 
taJrfto = -p— T^Uq. 



diù 

JJq croît quand co augmente. 

V étant le potentiel newtonien, on a sur la surface 

Si donc la surface rencontre l'axe de rotation, U^ est le poten- 
tiel newtonien au pale. Si (o croît, ce potentiel croît également. 
\]q ne peut d'ailleurs dépasser airR^ (page 16), R étant le rayon de 
la sphère ayant même volume que la masse fluide. 

En divisant membre à membre, on a 

iùJdiû rfUg 



2 



Si on fait croître w indéfiniment, co^ finira par dépasser 7:. 
W sera inférieur a to^J ; 



— J+W 
2 



'2 
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et par suite 








2 

3 


diù rfU, 



Si to croît indéfiniment, U^ devra croître indéfiniment ; mais 
Uq ne peut dépasser 27rR' ; il est donc nécessaire, ou bien que w 
cesse de croître, ou bien que la surface d'équilibre ne rencontre 
plus Taxe ; on a alors une figure annulaire. 

Nous verrons que parmi les figures d'équilibre, il y a des 
ellipsoïdes de révolution ou des ellipsoïdes a trois axes inégaux. 
Pour les premiers, tu est inférieur à 4*^ X 0,112. Pour les 
seconds, co < 4^^ X 0,098. 

Il est aussi possible qu'il existe une succession de figures 
d'équilibre, pour lesquelles w passerait par une infinité de maxi- 
mums ou de minimums. Ici, le raisonnement ne s'appliquerait 
plus et iù pourrait augmenter indéfiniment. 

Permanence de F axe de rotation, — Jusqu'ici nous ne nous 
sommes préoccupés que d'un mouvement de rotation uni- 
forme. 

On pourrait se demander s'il peut y avoir des figures d'équi- 
libre relatif pour une masse fluide animée d'un mouvement de 
rotation non uniforme. 

Nous allons prouver que c'est impossible. 

Nous pouvons d'abord supposer le centre de gravité du corps 
immobile; puisqu'il n'y a aucune force extérieure, son mouve- 
ment est un mouvement rectiligne et uniforme* 

L'équilibre ne sera pas troublé si on ajoute une nouvelle liai- 
son au corps, si on le solidifie. Alors le mouvement du solide 
autour de son centre de gravité est un mouvement à la Poin- 
sot. 

D'après le principe de Dalembert, le travail virtuel de la force 
totale, pour tout déplacement compatible avec les liaisons, est 
donc nul. Dans le cas qui nous occupe, il n'y a qu'une liaison, 
c'est l'incompressibilité du fluide, s'il s'agit d'un liquide. Si l'on 
a affaire à un gaz en équilibre, l'équilibre subsiste a fortiori en 
supposant le corps incompressible. 

Si Sx, 8y, S.s représentent les déplacements virtuels d'un point 
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x^ y, z, la condition d'incompressibilité s'exprime par la rela- 
tion 

rf % , d ^ , d ^ 
— T- ùx H — - — ou H — - — ùz = o. 
dx dy ^ ' dz 

11 peut y avoir trois sortes de déplacements virtuels : 

I® Les déplacements de toute la masse, considérés comme les 
déplacements d'un corps solide. L'application du principe de 
Dalembert montrerait que le corps doit avoir un mouvement à 
la Poinsot autour du centre de gravité, comme on l'a déjà vu. 

2^ Il y a les déformations du corps. 

3° On peut supposer que la masse soit animée d'uil mouvement 
tel que les molécules se déplaceraient sur des surfaces à densité 
constantes. 

Les surfaces d'égale densité étant aussi des surfaces équîpo- 
tentielles, on aurait 

dp ^ , d? -N , ^? -s 
^x ài/ ^ àz 

Les surfaces d'égale densité restant les mêmes, il en est de 
même de la forme extérieure. Tout se réduit donc à des dépla- 
cements internes dont nous allons étudier le travail virtuel. 

Soit un point M de masse //i, de coordonnées x^ 1/, Zy animé 
d'une vitesse x'^y'^z'y l'accélération du mouvement étant .r", y"^z", 
La force d'inertie a pour composantes mx"y my\ mz" , 

Appelons to^, w^, Wg les composantes de la rotation instantanée 
suivant les axes liés aux corps Ox, Oy, 0^. 

D'après les formules connues, on a 



x' 


== 


'^zU - 


- w,3, 


y' 


= 


■■ Wj V - 


-tù^X, 


^1 


=r= 


tô.,X - 


-^iU' 



L'accélération a pour composantes 

z" = (w^/r' — co^?/') + (toV>f — txi\y). 
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L'équilibre relatif ayant lieu, il y a équilibre entre le travail de 
l'attraction et la force d'inertie m^', mi/\ mz'^. 

Cette force se compose de deux parties : Tune est ce que Ion 
trouverait si l'on supposait le mouvement 
uniforme; l'autre est due h l'accélération 
du mouvement de rotation. 

A un instant t, l'axe de rotation est OP 
(wj, (i)jj, (O3) ; à l'instant t-{-dl l'axe instan- 
tané de rotation est OP' (fig. 5). Le vec- 

PP' 
teur OQ équipoUent au segment 

tend, lorsque dt tend vers zéro, vers 
une position limite, et ses projections sur les trois axes sont 

(W'„ W'2, (O'g). 

Si on prend cet axe OR lui-même comme axe Oz 

{iy'^= {ù[^= o. 

La force d'inertie due à l'accélération angulaire se réduit à 
^3^, — ^3^, o. 

Appliquons maintenant le principe de Dalembert aux déplace- 
ments virtuels que nous considérons. 

Le travail de l'attraction est nul, puisque la forme du corps n'a 

2 
pas changé. Le travail de la force centrifuge est égal à — oJ, en 

posant (0^ == (i)J + a)| -|" ^3- 

Ce travail est nul puisque J ne varie pas. 
Le travail de la troisième partie est 



^^2^ m{ijùx — x^). 



Puisqu'il y a équilibre, il doit être nul. Or, on peut toujours 
choisir un déplacement virtuel tel que 



2jn[ylx — xhy) 



ne soît pas nul; il suffit de supposer un courant autour de 0^. 
La somme aurait pour mesure l'aire du plan des xoy^ limitée 
par la projection du courant sur ce plan. 
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Il faut donc que tUg soilnui, c'est-à-dire que Ton ait un mou- 
vement de rotation unifopme. G. Q* F. D. 

Il faut de plus que Taxe de rotation soit un des axes principaux 
d'inertie ; cela résulte de Tétude du mouvemeut d'un corps 
solide. 

Stabilité de Vèquilibre. — Jusqu'ici on ne s'est pas préoccupé 
des conditions de stabilité de l'équilibre. 

Lejeune-Diriehlet a démoni;ré que la condition nécessaire et 

suffisante de la stabilité de l'équilibre est que W + J soit 

maximum. 

Mais il faut faire la distinction entre la stabilité temporaire et 
la stabilité séculaire. Lord Kelvin qui le premier a fait cette dis- 
tinction^ appelle stabilité séculaire celle qui a lieu même en 
tenant compte du frottement; tandis que la stabilité temporaire, 
qui subsiste quand on ne tient pas compte du frottement, ne 
subsiste plus si on en tient compte. Lord Kelvin a démontré que : 

La condition nécessaire et suffisante de la stabilité séculaire 
de l'équilibre relatif est la condition de Dirichlet. 

Evidemment cette condition est toujours suffisante ; mais sup- 

posons que W H J ne soit pas maximum, il pourra, d'après 

le théorème de L. Kelvin, y avoir équilibre stable, à condition 
qu'il n'y ait pas de frottement. Mais s'il y a un frottement, si petit 
qu'il soit, l'équilibre ne sera plus stable. 

On ne peut appliquer immédiatement le théorème de L. Kelvin, 
car il suppose que tout mouvement détermine un frottement, ce 
qui n'est pas le cas pour cette masse fluide : en effet, si la masse 
isolée dans l'espace se déplace d'un bloc, a la façon d'un corps 
solide, il n'y aura pas de frottement. 

Solide équivalent. — Avant d'aller plus loin, il est nécessaire 
de définir ce que nous appellerons le corps solide équivalent 
à la masse fluide. C'est un solide où, a l'instant considéré, 
les molécules ont la môme position que dans le système fluide. 

La vitesse de son centre de gravité est la même que pour le 
fluide. 
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Les axes principaux d'inertie ayant évidemment la même posi- 
tion, les trois moments de rotation autour de ces axes sont les 
mêmes que pour la masse fluide. 

Son mouvement est donc bien défini à Finstant considéré ; 
mais il faut remarquer que le solide équivalent à l'instant /, n'est 
"pas le solide équivalent a l'instant t\ 

Théorème. — La force vive du fluide est équivalente à la force 
vive du solide équivalent^ augmentée de la force vive du fluide 
dans le mouvement relatif par 7 apport à des axes invariablement 
liés au solide équivalent* 

En effet, appelons $, v), Ç, les composantes selon ces trois axes 
de la vitesse absolue de la molécule m ; $', vj', ^' les composantes 
de lu vitesse de la molécule correspondante du solide équiva- 
lent ; i ', r/', ^" les composantes de la vitesse relative de cette 
molécule par rapport au solide équivalent. On a 

Les trois forces vives considérées ont pour expression : 

T"=2,-r"(^"'+^"'+^"')' 

et il faut démontrer que 

T = T' + T". 
On a 

= T' + T"+ 2 '«{il" +-ri'yi"+%V')- 
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Je vais démontrer que 



2, '«(??' + 'iV' +!:'!;") =o. 



Cette expression représente le travail d'un système de force S'\ 
tel que la force agissant sur la moléciïle m ait pour composantes 
m^'\ mr/\ rriQ' , mais ce système est la différence de deux autres ÎS 
et S', dont les forces sont mÇ, /nv}, inC^i et mW my\'^ rnÇ ; ces deux 
systèmes de forces sont équivalents. La résultante générale est la 
même et les moments résultants des deux systèmes sont les 
mêmes par rapport aux trois axes, par la définition même du 
solide équivalent. 

On a en effet 

puisque la vitesse du centre de gravité du solide équivalent est la 
même que celle de la masse fluide. 
De même on a 

puisque les moments de rotation sont les mêmes autour des trois 
axes de coordonnées. Donc S = S' et S'' est nul; son travail Test 
aussi. 

Par suite 



2]/n(5'$" + VV'+^'C')=o. 



Condition de stabilité de Lejeune-Diricblet. — Revenons 
maintenant à notre masse fluide. 
On a Tégalité 

T — W = C^ 

D'où 

T+T' — Vf = CK 

Cette égalité suppose qu'il n'y a pas de frottement : autre- 
ment T'-|-T'' — W ira constamment en diminuant. 

PoiKCA.BE. — Figures d'équilibre. 3 



34 MASSE HOMOGÈNE FLUWB 

Si on pose 

T'+r'— w=6, " 

dt 
-j- est nul s'il n'y a pas de frottement. 

Mais alors c'est qu'il n*y a pas de mouvement relatif, et par 
suite T'' = o, 

Si T">o, §<o. 

Dans ce mouvement de rotation le théorème des aires s'applique, 
et on peut considérer comme des données les trois moments de 
rotation. 

Je vais démontrer que la condition nécessaire et suffisante de 
stabilité est que T' — W soit minimum. 

Pour une position d'équilibre relatif, on aura 

T' = o. 

Si l'on part de cette position d'équilibre, et que l'on fasse varier 
très peu les éléments des mouvements des molécules, mais de 
façon que les moments de rotation ne soient pas altérés (c'est 
une restriction sur laquelle nous reviendrons), la masse fluide 
ne s'écartera pas beaucoup de la position d'équilibre. 

En effet, soient T' et W, les valeurs de la force vive et de 
l'énergie, on a 

t'-w>t;-w.. 

On a a fortiori 

T+ T" — W = e est d'ailleurs voisin de T^ — Wo. 

Cette quantité ne peut jamais croître, donc elle restera toujours 
voisine de Tq — Wq, 

Les valeurs des variables ne différeront que très peu de leurs 
valeurs pour T^ et W^, donc il y a équilibre stable : la condi- 
tion est suffisante. 
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Elle est aussi nécessaire. En eflFet, si T^ — Wq n'est pas mini- 
mum, je puis chokir T' et W de manière que 



T'_W<T'— W„ 



et mèffîe que 



r+T«_W<T',— Wo, 



T" étant suffisamment petit. 

Mais alors T' -j- T'' — W ne peut que décroître ; T'' lui, tendra 
vers zéro ; T' et W ne pourront alors reprendre leurs valeurs pri- 
mitives, . , 

Il y a une restriction à levçr, Supposons que Ton écarte le 
corps de sa position d'équilibre sans conserver la valeur des 
moniei^ts de ratatixm autour des axes principaux; -on sait alors 
que le solide éqoivalent aurait un mouvement où les maments de 
rotation auraient des valeurs peu différentes de celles qu'ils 
avaient précédemment et Ton est ramené au pTemier cas,. 



Xr'axe de rotation est le petit axe de Vellipsoïde d'inertie. 
— Je vais démontrer maintenant que l'équilibre n'est stable que 
si l'axe de rotation est le 
plus petit axe de l'ellip- 
soïde d'inertie relatif a 
la masse fluide. 

Soient (fig. 6) Oxyz un 
trièdre fixe de coordon- 
nées ; OA,OB, OC, les axes 
principaux d'inertie," tù la 
vitesse instantanée de rota- 
tion dont les composantes 
sont (i>j, (1)^, (i>5 suivant les 
trois axes fixes, /?, y, r sui- 
vant les axes OA, OB,,- 
OC ; on a , 




Fig. 6. 



Soit OR l'axe instantané de rotation, J le moment d'inertie 
par rapport à cet axe. 



3Ô MASSE HOMOGÈNE FLUIDE 

On a . 

Soient u^, l'iiV'sy ^^s moments de rotation par rapport aux trois 
axes de coordonnées. Le segment OM, dont les projections sont 
([*i> [*2> I^s)» ^^* ^^® ^^^^ l'espace (théorème des aîres), et Ton a 

Les projectîons de OM sur les axes OABC sont A/?, B^r, Cr, 
et l'on a encore 

ix^=rAy-f-By+cv. 

La condition de stabilité, comme nous Tarons dît, est que 
T' — W soit minimum. Or, si on change l'orientation sans 
changer la forme du corps, W ne variera pas. Donc il faudra que 
T' soit minimum. OM est fixe dans l'espace. 

Si on pose 

p \/Â= X, 

7 v^r== Y, 

;• v/C = Z. 
On a 

T' = — (X* + Y^ + Z») 

{X* = (AX* +BY^ + CZ^) = C. 

Donc, le point X, Y, Z doit être sur un ellipsoïde donné et sur 
une sphère de rayon minimum. Donc l'axe instantané de rota- 
tion' est le plus petit axe de l'ellipsoïde donné et par suite, de 
l'ellipsoïde d'inertie. 

Si Â est la plus grande des quantités Â, 6, C, l'axe de rota- 
tion sera l'axe des x. Alors on a/> = w, y = o, r=:o. 

A = J, 

■ T' — .ïL 
aJ ■ 



VÉTERMmATlON DE VAXE DE ROTATION 3; 

Il faut donc que 

(i) W — T' = W ^ 

soît maximum ; donc que 

Il suffit aussi quQ 

soit maximum (p, lo), condition qui n'est pas nécessaire. 
Donc 







8(w + if 


')=» 


Les deux 


conditions donnent ; 








8W + i^ 


= 0, 






8W-1- "' 3J 
a 


= o. 



On doit donc prévoir que si la condition (2) est remplie, la 
première le sera également. 

Supposons qu'il en soit autrement, c*est-à-dire que l'on puisse 
avoir en même temps 






2J ^"» 2J/ 

Comme on a jiL=5(oJo, on conclurait facilement 

2 !a J 
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En additionnant membre à membre, on aurait ^ 

d'où . 

J'— aJJ, + V<o, 

(J-J,)»<o, 

ce qui est absurde ; mais la réciproque n'est 'pas vraie ; la 
deuxième condition n'entraine pas la première. 



CHAPITRE ni 

FONCTIONS SPHIÈRIQUES 



Expression de AV dans un système triple ' orthogonal. — 
Soient 

trois familles de surfaces, constituant un sjrstème. triple ortho- 
gonal. 

Considérons le petit parallélépipède rectangle dont les faces 
sont les 6 surfaces. 

Le« arêtes de ce parallélépipède-^ ont pour longueurs 

ad^y bdr\^ cdX^. 

La diagonale est 

ds^ = a^l'' + h^dr^^ + cH'Q, 

Supposons une fonction V ($, tj, Ç) et cherchons a déterminer 
la valeur de T'intégrale 



Jav^t 
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étendue au volume du parallélépipède; elle est égale a 



/ 



— rfcr 



étendue a sa surface. 

Le volume étant infiniment petit, la première intégrale a pour 
expression 

Xy ,a,h,c.d^.dr\,dX^, 

La seconde est la somme de trois intégrales étendues chacune 
aux deux faces opposées du volume. 
Pour la face $ =?= Ç^, on a 

Pour la face opposée $ = Ç^, + rfÇ, on a 

et par suite pour l'en^semble des deux facçs, on a 

On a donc 

AVaic.rfÇ d-n rfÇ 

.e j jv r f» *c ÔV d ca dV . ô ab dV 1 

et finalement on a * 

bc ôV 



AV— ' VA/ilil.\ 



Application aux coordonnées polaires. — Dans ce système 
on a 

* X = r sin 9 cos o, 

1/ :=r sin sin '^, 

z = r cos 0, 
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et, en posant cos 8 = jx, on a 

:jc == r \i — [JL^ cos îp, 
y =r\/i — ^^ sin o, 
z = rjji. 

L'élément linéaire dans ce système de coordonnées a pour 
expression 

ds^ = dx^ + dif + dz' = dr^ + j'^dH^ + /-^ sin^Orfy^ 
Au lieu de Ç, tj, Ç, nous avons/', [a, ? et dans ce système 



a= I, 



■ 


y/.- 


•l*^' 


c = 


7-v/l- 


-l-S 


abc 


= /•*. 




On calcule alors facilement 







Av d»V ^ 2 dV ^ I d r. ,. dV -1 



d»V 



Polynômes sphèriquea. Fonctions spbèriques, — On appelle 
polynôme sphérique un polynôme homogène en x, y, z satisfaisant 
a Téquation de Laplace,»AP =0. 

S il est de degré n en :r, y, z il contient ^ — coetti- 

cients arbitraires. Mais il y a un certain nombre de relations 
entre ces coefficients puisque Ton a AP„ = o. Cette relation est 

de degré n — 2. Il y a donc n ' — • relations entre les coef- 
ficients. 
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Il y a donc 

(p+i)(;i+a) h{n—i) ^^^ ■ ^ 

2 2 

polynômes sphériques, indépendante du degré /i. 

On peut exprimer Xy y^ z en coordonnées, polaires : on aura 

alors 

p «n Y 

Y„ étant uniquement une fonction de <p et de [jl, ~ que nous 
appelons fonction sphérique d'ordre /i. 
On a évidemment la relation 

(i) r^^nP,,. 

Considérons maintenant l'égalité (2) démontrée page 4 

Supposons que la surface S soit une sphère de rayon i, que U 
et V soient deux polynômes spliériques d'ordre rn et w, le second 
membre sera nul, et on aura 



/i 



Donc, ea tenant compte de l'égalité (i), onq: 

{m — n) j P„.Pn<iT = o. 
Sî m est différent de n, on à 

T«P„<f<r = o, • • 



/' 
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l'intégrale étant étendue à une sphère de rayah i ; et on peut 
encore écrire ceci : 

Si 7W = /i, ce résultat n'est plus vraF. 

Fonctions spbèriqnes fondamenUdea. — On a va que Ton a 



2 



./ 2 ./ —2 « +e 

:c = r y I — jjL^ cos ?p = r y i — [jl* 

y = r yi — jjL* sîn © = /• y i — i>? 



e'^—e"'' 



3/ 

Si on reipplace x et y par ces valeurs dans P„^ on aura 



p^n 



P,= 2/*e'^'XS, 



1»=-» 



Xj étant une fonction de [jl seulement ; ce sera un polynôme 
en ]x et en y/ i — jji^. Les polynômes où/? est pair, ne contiendront 
y I — |JL* qu'à des puissances paires. Les polynômes où p est 
impair, ne contiendront^/^—- jji* qu'à des puissances impaires. 

Les premiers seront donc des polynômes en jx, les seconds 
seretnt des polynômes en p. multipliés par y/ i — pi'. 

Substituons la valeur de P dans l'équation A P„ =*? o. 

On obtient : 

I. — . u* J 

En annulant le coefficient de ;•»-* e*?, on obtient l'équation 
à laquelle doit satisfaire Xj. C'est 



d r. ,, <fXS 1 



(/»4->)(i**-^0+y^' y,_^ 



I — fx» , 



44 FONCTIONS SPHÉRiqUES 

L'intégrale générale de cette équation est une fonction trans- 
cendante, mais il y a une intégrale particulière qui est de Tune 
des deux formes exigées par la définition de XJ, et qui est 
déterminée à un facteur constant près. 

On voit que XJ et X;^^ satisfont à la même équation. On a 
donc 

XS = AX;I^ 

A étant une constante quelconque. On a ainsi la forme des 
2/? -|" ï fonctions sphériques cherchées. 

X;, Xie^XJ^^^ , XSe'»'^.., 

Xi6->, X* «-*»>,...., x;;e—>..., 

On peut écrire autrement 2/i + i fonctions indépendantes en 
posant 

XJ, Xicoscp, XJcos2'^ ,..., Xjcos/Ey, 

X^sincp, XJsin 2cp ,..., Xj sin wcp. 

Deux fonctions sphériques difierent donc soit par w, soit par /?, 
soit que n et p soient égaux, mais que Tune contienne le sinus, 
l'autre le cosinus de /?y. 

On a déjà démontré que 






Y„rf'f 



est nul si m=^n; nous allons achever la démonstration et démon- 
trer que l'intégrale 






est nulle si Y et Y' sont deux fonctions différentes. 

Il suffit de le démontrer dans le cas de deux fonctions de 
même ordre. 

On a alors 

»2« 



7 xs Xid^ I , 



cos p'^ cos y^ d'f 
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OU bien 



Xj XJ rfjjL I sin />y sin y(p rf'^, 



/ XSXîrf[x / s 
/ XJ^ rf|x / si 



sin/?'f cosp'fdo. 



Dans les trois cas, la seconde intégrale est nulle. 
Enfin il est évident que 






est différent de zéro. 



Détermination des polynômes XJJ. — Proposons-nous main- 
tenant de déterminer Xj, et pour cela, commençons par le cas 
où p = o. 

L'équation est alors 

C'est l'équation qui définit les polynômes de Legendre, et 
on a 

Xn (Cos 8) est donc une fonction sphérique où 6 représente 
l'angle d'une direction avec l'axe 0^. 

Par raison de symétrie, X^(Cos e), où s représente l'angle de 
la direction considérée avec une autre direction quelconque, mais 
fixe, sera aussi une fonction sphérique ; mais on a 

cos Ê = cos cos 0'+ sin 8 sin 6' cos (o — '^j 
et cela, quels que soient 9 et 9'. 
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Donc 

X^ [ cos 8 cos 6' + sin 8 sin 8' cos (y — 'f ') j 

sera encore une fonction sphérique. 
Si on pose 

sin B' c-'V= ai, 
sin 8'e*f'=: 2/), 

d'où on conclut, eu ajoutant membre à membre après avoir mul- 
tiplié ces deux équations respectivement par e*? et e"" *?, 

sin 0' cos (cp — 'o') =.5e'? H-Tje""*? 

et en multipliant membre à membre, 

sin^e'=4iri; 

d'où 

cos'-6' =1 — 4^ ; 
la fonction 

XSi [v/i — 4irj [X— \/i— y!' (ie'> +r,c-'>) J 

est donc encore une fonction sphérique quels que scient \ et r^. 
En particulier, si on fait 

7) = O, 

est une fonction sphérique. 

Si je développe suivant la formule de Taylor, j'ai un poly- 
nôme : « 

X„ [;x + v/7^^e'']=-X„([.) -+- eh^/T^T^^ + ... 

1.2... p d^l.>' +-""^- I.2.../1. rfjA» 
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Chacun des termes de la somme est une fonètion sphérique. 
On a donc défini les fonctions Xjf a an facteur près. 



xs = (.-,.,-^^^ 



Propriétés des fonctions spbériques. — Soît U une fonction 
de 7% et P„ un polynôme sphérîque : 

Proposons-nous de calculer A (U P„). On a 

L-i ôU ôP,. 



mais 



dU rfU ôr dl] X 
AP„ = o; -——==—— ._— = _—. — , 
ox dr ox dr r 

A(UP„)=P,|-^+--^J 

• , 2 rfU / dP„ ^ dP„ ÔP„ \ 

et, parce que P" est une fonction homogène, 

A(UP..) = P„ (^^ + - -^) + — P» ^ 

Proposons-nous maintenant de calculer 

^(uy„)=a[^xp„] 

puisque P„ =/•** Y„. . 
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Couche sphéroïdale. — Nous allons appliquer ce qui précède à 
Tétude d'une masse fluide sphéroïdale) c'est-à-dire d'une masse 
fluide dont la figure d'équilibre est peu différente d'une sphère, 
ou, d'une manière plus précise, dont la surface est comprise entre 
deux surfaces sphériques de rayons voisins. Une couche sphéroï- 
dale sera une, couche comprise entre deux surfaces sphéroïdales. 

Action d'une couche sphéroïdale, — Soit (fig. 7) rfo- l'élément 
de surface d'une couche de rayon i, e l'épaisseur de la couche 

sphéroïdale, d'z l'élément de volume 
de la couche sphéroïdale ; on a 

On peut supposer que la masse 
contenue dans rfx est sur l'élément 
dTf et est représentée par [jl^t 

Soit un point M dont les coordon- 
nées cartésiennes sont x, y, z^ les 
coordonnées sphériques 

/', [JL = COSO, <p. 

Soit un point M' dont les éléments analogues sont x' , y\ z! \ 
z-' = I , [jl' = cos 8', cp'. 

La distance MM' est égale à i -j- /-^ — 2/* cos (OM, OM'). 

Cos Y = cos (OM, OM') == cos 8 cos 9' + sin 8 sin 6' cos (» — cp') . 

Développement en série de ^^., , — Le potentiel en M, dû à 
la couche sphéroïdale, est 




I MM' ' 



s' est une fonction de |x et de ç ; ;• est plus petit que i . 
Or 

MM" = 1 — 27- cos Y+ '•*= (i — re'i) (i — re-'i), 
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On a 

2.4 



^,_,,.-,)-;_,^.^ 



/•e»ï + 



1.3.5... (2n—i) 



2.4.6... 2/1 



C'est une sérîe absolument convergente pour /• < i , 
De même 

/ A-JL .1 • . . 1.3... f2w — i^ 

\ / 2 ' 2.4... 2/1 

Le produit de ces deux séries est une série encore conver- 
gente, si /• < I ; le coefficient du terme en z-** est 

. 1.3... f2/l — 1.3... f2/l — 3) I , ax 

2.4 .. 2/1 2.4... (2/1 2) 2 

. 1.3... (2/1 — 5) 1.3 ,, ... , , 1.3.5... (2/2 — i) 

2.4... (2/1 4) 2*4 2.4.0... 2/1, 

Le coefficient du terme en e^'^ égale celui du terme en e"^"^ 
On a donc un polynôme h coefficients réels du cosinus de 
Tangle y et de ses multiples jusqu'à n ; comme les coefficients de 
ce polynôme sont tous positifs, ce polynôme atteindra sa valeur 
maximum pour y=o. Mais alors : 

w = -^ = •+'•+••• + '^+- 

Donc les coefficients du développement sont inférieurs ou égaux 
à l'unité. 

On peut donc écrire 

Art étant une fonction de |jl, de cp, de jjJ et de cp'. 

Développement en série du potentiel d'une couche spbéroï' 
dalCp — Le potentiel V sera alors égal à 



^/•" / A„e'c/cr. 



PoiNCARÉ. — Figures d'équilibre. 
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r 

Les termes de V sont moindres que /•" | e'rfo-. 
L'erreur commise quand on prend pour V la somme des n pre- 
miers termes, est moindre que I e'rf<T. 



^/" 



La série sera donc convergente si /•<!. 

peut être développé suivant les puissances croissantes 
des variables x, y, z, 



MM' 



y/{a:' — xy^ijf'-yy^{z'—zy 



z^x^yPz^ j . 



Kpq^'dT^ 



et, si on revient aux coordonnées polaires, on aura un développe- 
ment suivant les puissances de /*"« 

Si j'appelle V„. l'ensemble des termes de degré n dans V, j'au- 
rai 

v=v.+v,+ ...+v„ + ... 

Or V satisfait à l'équation de Laplace : on a donc 

AVo + AV, + ... + AV„ + ... = o. 

Cela exige que chacun des polynômes V^ soit un polynôme 
sphérique; et, si on remplace dans V„, x, y^ z, en fonction de /'jJjl 
et cp, on aura un résultat de la forme /*" Y„ où Y„ est une fonction 
sphérique. On a donc 



^=2''^- 



Cette formule est valable si r est plua petit que i. 
Supposons maintenant r supérieur à i, et considérons le. 
point M^ image du point M, tel que OM^ X OM = i (fig. 8) 'r 
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5i 



j'ai OM == /•, OM, =— . Si M' est un point de la sphère, on a 
les égalités 

M'] 



^M^ ^ Q^ ^ QM' y/OMxOM^ /OM _ 

l'M, OM' OM. -y/oM'xOM, "V ÔM;-'"' 



Or le potentiel en M^ est 



On a donc 

r 

Or le potentiel en M,, point 

situé k une distance — du 
centre, est 

VA. 

Le potentiel en M est donc 

Supposons maintenant que M et M^ se rapprochent de la 
surface de la sphère en restant inverses Tun de l'autre. Les 
valeurs du potentiel M et M^ sont voisines, mais les dérivées 
ont des valeurs différentes de part et d'autre de la surface et 
on a 




Fig. 8. 



dV dV, 

àr àr. 



= 4ne; 



mais 



àr, ~l' dr "^ "^ J dr, 
= —Yr' 



dr 



/•". 
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On a donc lorsque /• tend vers i 

ôV 
Dans les mêmes conditions, puisque — r-f est fini, on a 

-. ôV . . 

lira -57 ('•—0=0' 

enfin 

lImV(i— /•*)== o. ■ 

SI on ajoute ces quatre égalités, on voit que l'on a 



Posons 



^'"[^\4r+^]=4- 



or 



Si r est inférieur à i, ^ peut se représenter par une série de 
la forme 



«ï.= 2(2«+i)/-»Y„ 



<ï> est donc une fonction satisfaisant a l'équation de Laplace, 
puisque c'est une somme de polynômes sphériques. 

Si r tend vers i, <I> tend vers 4""^^. Si on démontre que ce 
développement est encore valable pour /== i, on arrive au théo- 
rème de Laplace : 

Théorème, — Une fonction quelconque de deux variables est 
développable en une somme de fonctions sphériques. 

Potentiel en un point de la couche spbèroïdale. — Soit s la 
fonction donnée sur une sphère de rayon i. 
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Elle définit à Tintérieur de la sphère un potentiel V, prenant 
sur la sphère la valeur 47^6. 

Ce potentiel est développable en une série entière 



(0 ^Kr 



convergente si r est inférieur à i ; et, si on regarde /• comme une 
variable complexe, cette série définit une fonction <ï> convergente 
h l'intérieur du cercle de rayon i, et qui sera obtenue en faisant 
dans V, r=pe'î. 

Dans V, faisons r= e*»; on a une fonction <ï>' définie sur le 
cercle de rayon i, sauf peut-être en un certain nombre de points 
singuliers. Si la série (i) converge pour p = i, elle représentera 
^' d'après le théorème d'Abel. 

Je dis qu'il en est ainsi dans le cas où ^' est développable en 
série de Fourier. 

Nous avons alors 

Je dis que l'on a 

a;,=a, b,; = o. 

En effet, d'après le théorème de Cauchy, la suite des valeurs 
<t>' définit à l'intérieur du cercle de rayon i une fonction qui 
n'est autre que la fonction V (pe*-). 

Cette fonction est développable suivant la série (i), et on a 
d'autre part 



An = 



I / <ï)'rfcp I C^, n XJ 



les deux intégrales étant prises le long du cercle de rayon i. 
Mais la première intégrale se réduit k A'„, la seconde à B'„. On 
a donc 

A'=A„ 



B' = o. 



Donc, si ^' est développable en série de Fourier, la série sera 
encore convergente pour r=ze'^ et en particulier pour r= i. 
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Tout revient donc h démontrer que 4>' est développable en 
série de Fourîer, c'est-à-dire que <ï>' est une fonction de © admet- 
tant une dérivée, sauf en quelques points singuliers, et qu'en ces 



points singuliers l'intégrale / 1 <ï>' | rf'^ a 



un sens. 



Nous allons maintenant démontrer que e est développable en 
une somme de fonctions sphériques. Le potentiel en un point M 
de coordonnées (r, 8, 'f ) est 

v= / -^^ 



/ t'd<j' 
MM" 



e' étant la densité de la couche sphérique, au point (0',»'). 
On a 

rfMM' 



d/- ~ / MM'* 

■/ 



dr 



t e' (cos y — r) d^ 

Comme cela a été démontré, nous avons 

e==$(0, cp)=V+2r-^= ■ — ^^ 



5V_ r t'ii—r") 
d/- ~ / MM'» 



Supposons que OM soit une direction fixe que nous allons 
prendre comme axe O^ : ^ i^ 

9 = 0, ix=i, 
MM'= I — arcosY+r^=i — ajxV-f-r*. 

L'élément de surface est dif =.d^d':^\ et on a pour valeur en M 
de la fonction <ï> 

*,=(,-,•») / ^-^ j 1 e'i?', 
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car e' seul dépend de ^'. Si on pose 

F ([.') = 



on a 



Supposons que la sphère soit partagée en un nombre fini de 
polygones sphériques dont les côtés soient des arcs de courbes 
analytiques, et que dans chacun de ces polygones e' soït une fonc- 
tion continue de \k' et de cp '. 

Dans ces conditions, F ([jl') sera une fonction analytique de jx', 

excepté pour les valeurs de \k' correspondantes aux sommets des 

polygones, ou aux parallèles de la sphère tangents aux côtés des 

F (a') 
polygones. La fonction .. * ^^ sera alors intégrable entre — - i et 

Si on pose r = e*^, on aura 

et si on pose ensuite 

cos^» — [x'==Z% 
on aura 






Nous démontrerons plus loin que <^ est une fonction de i^ 
admettant une dérivée, et est par suite développable en une série 
de Fourier. 

Enfin, si nous développons suivant les puissances de /% nous 
aurons une série convergente pour|r|<i. Elle sera sûrement 
convergente même si |. /• | = i /dans le cas où e ne comprend qu'un 
nombre fini de termes. 
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Voyons quelques exemples qui nous serviront par la suite. 
Si on fait£' = , on aura 



2TC 






F(t^')= / 4!-=^ 



si on fait e' = — , on aura 

On pourra prendre pour fonctions sphériques correspondantes 
Y,=-î-ouYo = — . 

De même, si e' = -i— , on a en appliquant la formule de la page Sa 

4«e' = Y.+3Y.4-5Y. + ... 
Dans le cas qui nous occupe 

Yo = Y, = ... = o 

et 

Théorème de Laplace. — On a toujours le droit d'écrire 

F(f.') = a+i[x' + F,(fx') 
a et i étant des constantes quelconques ; ^ sera alors de la forme 

et si on pose 

f]\ ' — ''^ — isîntt _ ^ 

'nW= 3-= -=^6-'—; 



on aura 



», =, (« 



,rïMx 
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Ceci posé, pour établir que ^ est développable en série de 
Fourier, il faut démontrer : i® que — 77- existe et est fini sauf en 
un nombre limité de points singuliers de la fonction ^ ; a** que 



/ 



|<I>|eZf|) est fini en ces points singuliers. Il y a une difficulté, 



carZ s'annule si [jl' = cos ^. Il faut donc que F ([x') s'annule pour 
jjl'==cos^ et alors on doit étudier vers quelle limite tend la 
fonction lorsque jx tend vers cos ^, 

I® — 77- existe. En effet on a 

Or, je puis disposer de a et i, de façon que Ton ait 

. Fj (cos ^) == o, 
F( (cos tj;) = o. 

Si donc cos ^ — jx' est un infiniment petit du premier ordre, 
F(jJi'), qui est analytique au voisinage de cos ^, sera un infini- 
ment petit du second ordre, et, sous le signe d'intégration, les 

3 5 
dénominateurs sont respectivement d'ordre — et — ; le numé- 
rateur est d'ordre 2. L'intégrale est donc finie, sauf si cos ^ est 
une valeur singulière de F([x'). Il faut qu'en ces points singuliers, 
où <j^ = ^^y l'intégrale 



/ \^\d^ 



soit finie. 
Supposons 

a = F (cos ^q) 
i=o. 



-^^K fllrf^. 
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On a alors 

/ |<ï.|i<^<a / lF(cost{-)l4+ / n4 / l!Mz± 

Le second membre doît rester fini, La première intégrale est 
finie. Il faut démontrer qu41 en est de même de la seconde. 

C'est une intégrale double qui, aux points singuliers, doit être 
un infiniment petit d'ordre inférieur à 2, et aux lignes singulières 
d'ordre inférieur à i . • 

Or, nous avons ici une ligne singulière, la ligne ja' = cos i^^ 

3 
mais pour [x' = cos i^^ le dénominateur est d'ordre — , le numéra- 
teur est d'ordre i, sauf aux points singuliers de F(a'). Donc 
l'élément de l'intégrale devient infini de l'ordre — , l'intégrale 

a un sens. 

Aux points singuliers F([jl') =^ F(cos ^), la fonction est d^ordre 

3 

— , l'intégrale a encore un sens, 

La fonction <ï> remplit donc les conditions nécessaires pour 
être développable en série de Fourier. 

La formule de Laplace est valable pour /'= i, et l'on a 



-iz^' 



(2/l+l)Y„. 

Le théorème de Laplace est ainsi démontré. 



CHAPITRE IV 
MASSE FLUIDE HÉTÉROGÈNE. PROBLÈME DE CLAIRAUT 



Préliminaires. — Supposons le mouvement de rotation de la 
masse suffisamment lent : la surface différera peu de celle d'une 
sphère, le carré de la vitesse de rotation a>^ sera considéré 
comme un infiniment petit du premier ordre. 

On a déjà vu que si le mouvement de rotation est nul, les 
surfaces équipotentielles sont également des surfaces d'égale 
densité. Le potentiel Vq en un point est uniquement fonction 
de r et l'équation de Laplace se réduit à 

La masse contenue entre deux sphères infiniment voisines de 
rayon r et r-^-dr est 

et la masse totale contenue à Tintérieur de la sphère de rayon r 
est 



/ 



47C/'prf/*. 



Si on appelle D la densité moyenne de cette sphère 



/ 



4TO-\s<i/- = -|- TO-»D. 
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D est une fonction de r dont la dérivée D' est négative dans 
Thypothèse d'un équilibre stable. 

En différentiant cette dernière équation, on a 

En supprimant le facteur ^rcr^dr, on a 

(2) 3p==3D + rD'. 

En prenant encore une fois la dérivée, on obtient la rela- 
tion 

(3) 3p' = 4D' + /-D'' 

qui nous sera utile par la suite. 

Faisons encore la remarque suivante. On a évidemment : 

3D^3D + /D'^o. 

La première inégalité est vérifiée puisque D' est négatif; la 
seconde Test, car p est positif. On a donc 

0^—^-3. 

On a -7r=r- =0 SI la masse est homogène. 

On a "Yp = — 3 si toute la masse est concentrée au centre. 

Le potentiel de la sphère de rayon /• en un point de sa surface 
est 

La force qui s'exerce sur une molécule de la surface est 

(4) -^_--y-«,D. 

Puisque p et V,, sont uniquement fonction de r, on peut les sup- 
poser fonctions l'un de l'autre, et, en vertu de la relation (i), 
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AV,, est fonction de V^, et on peut écrire AV^, =f[\^) ; en diffé- 
rentiant 

df{N,) _ d(A\) _ dr _ -47tp' 

dr 6 

En tenant compte delà relation (3), on a 

4D' + rW 



(5). r(Vo) = - 



/•D 



Développement du potentiel en série. — Supposons main- 
tenant un mouvement de rotation lent. La surface se déformera, 
mais peu. Le potentiel en un point d'une sphère sera dévelop- 
pable en une somme de fonctions sphériques^ et Ton aura 

v=v,+2hy, 

Y représentant une fonction sphérique et II étant un coeffi- 
cient qui est une fonction de /• de Tordre de w^, puisque V dif- 
fère peu de Vj,. 

L'équation de Laplace peut s'écrire 

(6) AV = AV,+^A(HY) = -4«p. 

On a d'ailleurs, d'après une formule démontrée plus haut 

en supposant que n soit l'ordre de la fonction sphérique Y. 
Si on pose 

la condition d'équilibre est, comme on Ta démontré plus haut, 
que les surfaces U = constante coïncident avec les surfaces 
d'égale densité. 
On a donc 

AV = — 4^? = fonction de U. 
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0)^ étant très petit, V diffère peu de V,, ; 

AV diffère peu de AV^,. 

La relation entre AV et V diffère peu de la relation entre AV,, 
et V,, et on peut par suite écrire 

AV=/(U)-|-»(U), 

cp étant une fonction dont la valeur est constamment de Tordre 
de (o^. 

Si Ton développe le second membre suivant la formule de 
Taylor, on a 

AV = /-(V.) + (U-V.)/'(V.)4-... 
+ î(V.)+(U-VOî'(V.)+... 

En vertu des hypothèses faites, on peut négliger les termes 
infiniment petits du deuxième ordre et écrire 

AV =/-(V.) + (U - V,) f (V.) + ï (V,) 
= AV. + (U-V.)/'(V„) + cp(V„). 

Enfin, en tenant compte de (6), on peut écrirç 

2A{HY) = (U-V.)/'(V,)+'f(V„). 
On a d'ailleurs 

— [oc^-^y^) est lui-même décomposable en une somme de 
fonctions sphériques 



,2 



~ [s^+f) ^Ji_(^»+y^+.-^)+:^(^«+j,^_..=) 



~3 '""6 1^ 
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Le développement se réduit à la somme des deux fonctions 
sphériques Y = i et Y = ^^ affectées respective- 

ment des coeliicients L. =^ — r; — etL« = — -r. — '. 

Les autres coefficients sont nuls. 
On peut dans ces conditions écrire 

U-V,=2,^H-|-C)Y, 

^A (HY)=2(H + C)Y/'(VJ + 'f (V„). 

Détermination des coefûcients du développement, — Pour 
trouver les valeurs de H, on peut égaler terme à terme les coef- 
cients d'une même fonction Y. 

A(HY) 

contient Y en facteur ; on a donc 

i(HY)=Y(IH-C)/'(V,) 
Si on prend Y= i, on a 

(8) .h_(h+^,-)^2:+«:+«(vj, 

H est donc déterminé par une équation linéaire du second 
ordre. 

Si on prend 

. V x^ + y^ — ^z^ 



H est déterminé par l'équation : 

,9) .(HV) = (H+iî^)Yi£:±i2: 



64 MASSE FLUIDE HÉTÉROGÈNE. PROBLEME DE CLAIRAUT 

Pour les autres, on a simplement Téquation 

Ces dernières équations ont une intégrale évidente, H = o; 
nous allons démontrer que Ton ne doit pas en considérer d'autres. 

Dans le cas où Y est une fonction du premier ordre, on peut 
trouver la valeur correspondante de H. En effet, il y a trois fonc- 
tions indépendantes du premier ordre : 

X y z 
~' "7' 7"' 

Considérons Tune d'elles, -^ par exemple, on a vu que Ton 

avait 

AV,=/-(V.). 

Si je prends les dérivées des deux membres par rapport à Xy 
j'ai : 

•^W)=^/(V„), 
dV.\ „„,,dV„ 



^(•S^)=nv.,^, 

L dr ' /• J dr " ;• /D 

X 

On voit que si Y est égal a —, Téquation 
A(HY) = HY'-"" + ^P' 

est satisfaite pour H = a — r-^. , 11 en est de même si Y = -^ 
^ ar r 

ou 1 = — » 
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On a alors pour les termes du premier ordre 



dr \ r )' 



fl, i, c étant des constantes quelconques. — Nous verrons pro- 
chainement qu'on peut les supposer nulles, à condition de 
prendre le centre de gravité de la masse comme origine des 
coordonnées; 

Les surfaces de niveau sont des ellipsoïdes. — Considérons 
une surface de niveau primitivement sphérique, lorsque la rota- 
tion est nulle. 

Pendant la rotation (fig. 9), 
le point M est venu en M' sur la 
nouvelle surface de niveau. Le 
déplacement normal est 

j; = M]Vrcos(7-, MM'). 

Avant la déformation, on avait - 
en M Fig. 9. 

Après la déformation, la densité en M' est encore p. La den- 

AV 
site en M est 0^= — -- — - ^ en appelant V le potentiel en M, et 

on a la relation 

AV=AV, — i; ^^ 




On a de plus 



dr 



JAV. dAY. rfV.. 



dr ~ d\^ • dr ' 
et, en se servant des relations (4) et (5), 

av=av.+-|Li;(4D'+,-d"). 

PoiNCARÉ. — Figures d'équilibre. 
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Nous avons d'ailleurs déjà écrit 

On a donc 

Revenons maintenant aux équations déterminant les coefli- 
cients II. 

L'augmentation de volume entre la figure d'équilibre au repos, 
et la figure d'équilibre dans le mouvement de rotation, est nulle. 
On a donc 



ou 






AH\W<r = o. 
Ceci peut s'écrire 



2?('') / Yrfcr=o. 



Chaque terme s'annule séparément, sauf toutefois le premier, 
pour lequel Y est égal à l'unité. On a donc, en désignant le coef- 
cient de Y^^ = i , par H^^, 

mais IIjj est uniquement fonction de R si bien que 
I/intégrale générale de cette équation est 



A et B étant deux constantes. 
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Mais remarquons que V = Vo + SHY; a Tinfini, V doit être 
nul comme V^. Donc A est nul. De même B est nul ; autrement 
le potentiel au centre serait infini, ce qui est absurde. 

L'équation (8) qui définit H se réduit dans ce cas à 

/■'(V„)-^ + 'f(V,) = o. 

Cette équation définit o (VJ. 

11 nous reste encore à voir que pour les autres valeurs de Y, les 
seules valeurs possibles de H sont les valeurs H = o. Si nous 
admettons ce résultat, nous voyons que ^ se réduit a un seul 
terme, et que Ton a 



On a donc 



g- "•"-!:-"• i M. 



Ceci montre que les- surfaces de niveau sont des ellipsoïdes. 
En effet, on a 

{jc'-h,/ + z^-) == (r H- Kf = r' -+- o^rX, + ij^ , 
mai^ nous pouvons négliger l^^, et il reste 

C'est Téquation d'un ellipsoïde de révolution autour de Taxe Oz, 

Aplatissement, — L'aplatissement est le rapport de la 
différence entre le rayon équatorial et le rayon polaire, au 
rayon moyen, que nous pouvons supposer égal à r. 

Pour avoir ces rayons nous pouvons calculer Ç pour les valeurs 

(^x = y = Oj z == /•) et pour (z/ = ^ = o, ^ = /), on trouve les 
deux valeurs respectives 



M I 



■v 



68 MASSE FLUIDE HÉTÉROGÈNE. PROBLÈME DE CLAIRAUT 

La différence est 

et l'aplatissement est 

e= Li-^ = 34 (r) 

3/2; _ K 

— x^^tf — iz" rY' 

en désignant par Y la fonction sphérique 





x'^-^if — ^z' • 


On a vu que 


3 2AI1Y 

^~ 411 /•D" + 4D'' 



La somme se réduit ici à un seul terme. L'équation (9) donne 

Donc, 

4™-»DeY = (lH-^)Y, 

Aer^DY= -^[^11 Y + -^Ar^Yl. 

On a d'ailleurs aussi 

A/'Y = o, 

puisque Y est une fonction sphérique du second ordre. 
On a aussi 

-9- AHY = f/D" + 4D') e;Y. 
L'équation devient donc 



\ 
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On a démontré que si P est un polynôme sphérique d'ordre /?, 



on a 



A(up)=pru"+2(/i+i)-^i. 



Si nous posons pour n = 2 

P^r^Y, 
On aura Téquation 



en développant les calculs on trouve, 

e"\) + 2e'D' + ^ eD' + — f^'D = o. 

C'est une équation linéaire du second ordre que Ton peut 
ramener à une équation du premier par le changement de 
variables. 



= 'V+^' — VI- 



e 
L'équation devient 

(12) D(r/i' + V + 5-/i)+2rD'(i +•,)) = o. 

Cest l'équation de Clair aut. 

Limites de -ri. — Nous avons vu (p. 60) que Ton avait 

-3<^<o. 

Dans Téquation (12), substituons ces valeurs extrêmes de rD\ 
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et considérons les deux paraboles (fig. lo), oii l'on a porte en 
abscisses y; et en ordonnées rr\' 



et 



/v/ + V + 5-/i-6(i+70 = o. 




Fig. 10. 

Le point dont les coordonnées sont /v/et rj ne peut se trouver 
qu'entre les deuxparaboles (dans la région hachurée sur la figure.) 

Je vais maintenant démontrer que vj est compris entre o et 3. 

fi ne peut être négatif. En elTet, s'il est négatif pour une 
valeur r^, on pourra trouver deux nombres positifs a et ^ tels que 
l'on ait 

— ?<vi< — a, 
a et p pourront même satisfaire a l'inégalité 
o<a <2 <5 <p. 
Je dis que fj ne peut dépasser — a. 



LIMITES DE fi 71 

En effet, si yj atteignait le valeur — a pour une valeur /• infc- 

rieure a /\„ -^ — serait negatil, /• --7-^ serait negatil. 

Le point dont les coordonnées sont /t/ et yj = a ne se trou- 
verait pus alors dans la région hachurée. 

De même r^ ne peut dépasser — 3. En effet, si pour une valeur 

, . p . , . rfvi . . .« rdT^ 

de /• mlerieure a/-^ on avait /• < /*,,, --7-^ serait positit, —y- serait 

positif aussi et le point représentatif serait dans une région 
du plan où le point [i^\ t)) ne peut se trouver. 

Donc si pour une valeur de r^^ yj est négatif, il est, pour toute 
valeur de r inférieur a r^, compris entre deux nombres — aet — [d, 
on a donc 

~< — a; 



a 
< . 



En intégrant de r^ à r </o, on a 



L — > L 



^0 



(t:)-' 



e > A/— *. 

Or e croîtrait plus que lorsque /* tend vers zéro, ce qui 

est impossible pour /• suffisamment petit, car e serait très 
grand. Donc yj ne peut être négatif. 

Je dis maintenant que */; est inférieur à 3. Le raisonnement est 
le même. 

En effet, supposons que pour /•=: r^ on ait */; = Tj^ > 3, je dis 
que pour /• inférieur à z'^, on aurait encore vi > 3. 

En effet, l'inégalité ne cesserait que si yj devenait inférieur 
à 3, ou si T, devenait négatif en passant par l'infini. Cette dernière 
hypothèse est à écarter d'après ce que Ton a vu plus haut. Il 
reste donc a montrer que y^ ne peut devenir égal a 3. En effet, 
pour yj = 3, 7/ serait positif, /-yj' serait positif, ce qui est absurde, 



7a MASSE FLUIDE HÉTÉROGÈNE, PROBLÈME DE CLAIRAVT 

le point (/-y/, Yi) ne serait pas compris entre les deux courbes. 
On a donc toujours 

Y, <3. 
On a d'ailleurs 

^fî + 'i^ — '^ — ^<^ V^^^ /i > 3. 
A fortiori on a 

/•V + '^i'— -0-6 — 5(y,— 3)<o, 

/V + C/i— 3)^<o 

et en multipliant par dr qui est négatif si r décroit de /^ à zéro, 
on a 

dr , di\ 



c'est-a-dire 



dh h ^ T-<o, 

/' YJ Ô 



r Y,— 3 r^ fïo — 'i 

Le second membre est une quantité finie ; le premier est une 
somme de deux termes, dont le premier seulement devient infini 
lorsque r tend vers o, l'inégalité est donc impossible. 

On a donc bien 

o<vi<3 

C. Q. F. D. 

Les limites peuvent-elles être atteintes? 

Si on a •■ - = o, Téquation se réduit à 

,Y+ 71^+571 = 0, 

et admet la solution 

y, =^ o. 



LIMITES DE T, 73 

Si == — 3, Téquation se réduit à 

'V + (-1-3)(-o+2) = o 

qui admet la solution •/) == 3. 

Le premier cas correspond à un point intérieur de la masse 
supposée homogène ; le second correspond aux points extérieurs 
à la masse, et en particulier au cas où toute la masse est concen- 
trée au centre. 

Je dis que dans ces cas particuliers, les seules intégrales 
admissibles sont les solutions y) = o et v^ = 3, 

En effet : 

'D' ^,, . ,, . , 

bupposons — — - = o. L équation se réduit a 

/V+Vi^ + 5vi=o, 
dr dr^ 



dont la solution est 



5 r, + 3 



Pour /• = o, on aurait 

L 



7, + 5 



donc 71 = — 5 ; or pour /• quelconque, 't\ est toujours positif : 
donc la constante d'intégration ne peut être que — 00 , et l'inté- 
grale se réduit à y; = o. Il en est donc de. même pour 7^ = 3, 
dans le cas de la masse concentrée au centre, A l'extérieur de la 
masse, on a l'équation 

/V+vi^ — VI— 6=0, 
dont l'intégrale générale est 

La constante d'intégration ne peut être que o et on a y» = 3. 
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Ce sont les seuls cas où Tj puisse atteindre ses limites. En effet, 
supposons d'abord que pour une valeur i\y on ait -/; = o, pour r 
différent de zéro, r, est supérieur à zéro, donc y/ (/J = o, et en 
portant ces valeurs dans l'équation de Clairaut, il vient 

D' = o, 
mais comme 

/•D'+3D = 3?, 

on a D= p. La densité en ce point est donc égale à la den- 
sité moyenne, et comme par hypothèse la densité ne peut 
diminuer vers le centre, on a un noyau homogène. 

Supposons de même que pour r = r^ on ait r^ = 3, de 
même que précédemment il faudra que y/ = o ; donc 



D 



p = o en ce point, et pour /• > t\ on a encore p = o ; p ne peut 
augmenter puisqu'il y a équilibre, et par suite y^ = 3 pour /• > /y. 

Il reste à voir ce qui se passe pour /• < r^. Je dis que Ton a 
encore y, =3. S'il en était autrement, y^ serait plus petit que. 3, 
/•^ (y, — 3) serait une quantité qui diminuerait avec /*, sa dérivée 
serait positive. On aurait donc 

/V+5(y, — 3)>o. 

Or, l'équation de Clairaut peut s'écrire 

,V+5(/i-3) + {t.-3)'+[^ + 6]|"i4-vi]=o. 

Les deux derniers termes étant essentiellement positifs, il en 
résulte que l'on devra avoir 

/•'V + 5(y.— 3)<o 

ce qui est contraire à ce que l'on a établi plus haut. 
Pour /• < i\y on a donc encore y^ = 3, et par suite 



[^+a]=„. 
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Donc p = o pour r < r^, 

La masse de la planète est donc entièrement concentrée au 
centre. 

Forme des courbes intégrales. — Etudions maintenant la 
forme des courbes représentées par l'équation * 

,V+V + 5^+^' (' +'^.) = o. 

Traçons les axes de coordonnées Or, et Or. 
Il suffit évidemment d'étudier les courbes comprises dans les 
régions où Ton a 

0<T4<3. 

Traçons la droite D d'équation r^ = 3. 

Par un point du plan passe en général une courbe et une seule, 
à moins que t/ ne soit pas une fonction holomorphe. 
Or, on a 



7/ n'est pas holomorphe si r = o ; la droite r= o est une courbe 
intégrale de l'équation ; il y a sur cette droite deux points remar- 
quables Yj=oet7i= — 5, pour lesquels r\' est indéterminé 

rD' 
à condition que — =-- tende en effet vers zéro avec /•. Or 

rD' 



;r, on a 



r = 3(^-.) 



Le second membre est essentiellement négatif, la densité d'une 
couche étant inférieure à la densité moyenne de la masse inté- 
rieure à cette couche, puisque par hypothèse il y a équilibre. 

Lorsque /• tend vers zéro, p tend vers une limite p^ et l'on a 

P<D<?o, 

P P 



76 
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mais -^ I tend vers zéro avec /*, donc -^ i tend aussi vers 

zéro et par suite — =— tend vers zéro avec /*. 

Aux deux points 

( Y, = o 

( Ti = 5 

les deux courbes intégrales ne sont pas définies d'une façon 
unique. 

Faisons maintenant /• suffisamment grand, alors 



L'équation devient 

rf/ + r,^ — r,— 6=0. 
Elle s'intègre immédiatement, et son intégrale générale est 



/•*(v.-3) 



:C\ 



Les courbes représentatives sont toutes asymptotes à la droite 
Yj =3, mais elles ne sont pas nécessairement asymptotes à la droite 

/•= o, car on a -rr- =p — 3 pour 

/• == o : le raisonnement ne s'ap- 
plique plus. 

Considérons le rectangle OABC 
(fig. Il) formé par les droites 
Y, = G, r = o, Y. = 3, r = r„, /-, 
étant une valeur quelconque. Par 
un point P, pris sur le côté OA, 
passe une courbe intégrale et une 
seule. Pour cette valeur 




Fig. II 



v= 



D 



■>o; 



la portion de la courbe qui pénètre dans le rectangle est donc 



P' 
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dirigée vers le côté AB, quel que soit le point P. Une courbe 
intégrale ne peut donc sortir du rectangle par le côté OA, car si 
elle ressortait en un point P', elle aurait une forme analogue h 
la courbe PMP' (fig. 12), ce qui est impos- 
sible puisque -^ — ne peut être nul dans le 

rectangle. 

Donc cette courbe ne peut sortir du rec- 
tangle que par un point M de la droite AB, 
De même, si l'on considère la courbe inté- Y\^. 12. 

grale passant par un point Q du segment 
GB, la tangente en ce point a comme coefficient angulaire 

dr r 



d'f\ 4 



f[^-^] 



quantité qui est négative puisque Ton a 

2/D' 



D 



6. 



La branche de courbe qui pénètre dans le rectangle est donc 
dirigée vers l'arc AB, le même raisonnement que précédemment 
montre que cette courbe ne peut sortir du rectangle que par un 
point N de la droite AB. Evidemment, d'ailleurs, le point M est 
à gauche du point N puisque deux courbes intégrales ne peuvent 
pas se couper. 

Les points M ont donc une limite M' lorsque le point P se rap- 
proche de G, de même les points N ont une limite N' lorsque le 
point Q se rapproche de G. 

Le lieu des points M' sera" une courbe complètement comprise 
entre les droites 7^ = 0, yj =3. Il en est de même pour les lieux 
des points N'. 

Est- il possible que M' et N' coïncident ? Supposons, ce qui 
sera démontré un peu plus loin que cela arrive. Alors il y 
aura une courbe intégrale et une seule qui passe par le point O. 
Gette courbe sera évidemment la courbe intégrale convenable, 
puisqu'elle restera constamment comprise entre les droites y; = o, 
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r, = 3. Il en résulte que la fonction r» (r) est parfaitement dcter- 

re' 

minée, est donc parfaitement déterminé ; e est alors, en 

e 

désignant par E la base des logarithmes népériens 

La constante e^ sera connue si on se sert de Tune des équations 
précédentes, par exemple 



-ti"-^)- 



Une seule courbe est acceptable. — Démontrons maintenant 
qu'il y a une seule courbe intégrale acceptable. Nous allons, pour 
cela, démontrer qu'il y en a une et une seule qui passe par le point 
/•iz^o, 7i = o, et que toutes les autres passent par le point / ==o, 

Démontrons d'abord qu'il y a au moins une courbe passant 
par le point tj = o, et une passant par t) = — 5. 

Nous avons déjà vu quel était le signe de /-/i' pour diverses 
valeurs de rj ; ry^' est négatif si r^ est inférieur à — 5, ou supé- 
rieur à + 3 et positif si tj est compris entre o et — a. 

Dans les autres intervalles, les deux signes sont possibles. 

Pour ri=z — 5, rvj' peut être nul ou négatif, mais non positif. 

Pour 7j = — 2, /t/ ne peut être que positif ou nul. 

/•D' 
Enfin, nous avons vu que — r;- — est une fonction tendant vers o 

avec r, de façon, qu^étant donné un nombre positif a voisin de 
zéro, on peut déterminer un nombre j\ tel que l'inégalité 

entraîne 

o>-^>-a. 

L'équation de Clairaut peut s'écrire 



LNE SEULE COURBE EST ACCEPTABLE ^9 

^ et Y étant des nombres positifs, tendant Tun vers zéro, l'autre 
vers 5 lorsque a tend vers zéro. 

Ceci va nous donner des indications sur le signe de r//, lorsque r 
est inférieur h /•,. On a le tableau suivant: 

Il y a deux intervalles pour lequel le signe de /t/ est douteux, 
mais ces intervalles sont moindres que ceux précédemment 
déterminés. 

Ceci posé, je vais démontrer que pour /•<r,,T4 est nécessaire- 
ment inférieur a p. Supposons que pour une valeur r^<r^ on ait 
Yj == r^^> p. Je dis que si r < j\ on aura encore r > p. En effet, 
l'inégalité ne pourrait cesser que si t/= p, mais alors dans Tin- 
tervalle, T^'sera positif, ce qui est contraire à la règle des signes. 

Or Téquation peut s'écrire : 

'V+(t.— P)(r. + y)+2(-^ + a) (i +7i) = o; 

. par hypothèse 

;D' 

■> — a, 



D 



on a donc 



A fortiori on a 



'^/+(^-P)(^ + T)<o. 



,V+(71 — P)^<0 



'' I 



et, si on intègre de i\ à r<ï\, 



— <o; 



T I ' T ' I 

— L r>— ^ ô 

L 1 7r<L 



fi '-1 "^-^ii-? 
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Or Tj — j3 serait positif ; on aurait a fortiori 

_ I , I I 

L — <L \ ^. 

Le second membre est fixe, le premier peut grandir au delà de 
toute limite, Thypothèse faite est donc absurde. Par suite 

7^ a pour limite zéro lorsque r tend vers zéro. 
Je dis maintenant que l'équation de Clairaut admet au moins 
une intégrale tendant vers — 5. 

En effet, étant donnée une valeur '/j^ comprise entre — y et — 5 

— 5 <•/),<— Y, 

considérons une courbe intégrale passant par le point (^^,7,,), je 
dis que si i\ < i\y r\^ <y. Autrement r/ serait négatif pour t^ > — y, 
ce qui est contraire à la règle des signes. 
On a donc pour r suffisamment petit 

— 5<o<— y; 

mais y pourra être choisi aussi voisin de 5 que Ton voudra ; il en 
résulte quer, tend vers — 5. 

Soient t,q et rj^ les deux intégrales dont nous venons de démon- 
trer Texistence. Il leur correspond des intégrales e^ et e, de 
Téquation linéaire du second ordre 

4e"D + le'D + -i^ D' + -^ e'D = o. 

L'intégrale générale de cette équation sera 

A et B étant deux constantes. L'intégrale générale de l'équa- 
tion de Clairault sera 

_j^ Ar^Z + Br^î 



e 
A-/) 




A 


^'' +B 
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Or, quand r tend vers zéro, r\^ tend vers — 3, 7^^, tend vers zéro. 
Il en résulte que r\^ — r\^ pourra être inférieur à — 4» On aura 
donc 



'•«/ 


re' 


«1 


e ^' 


< 


e'„ ^ 4 




«.< 2-, 


«1 




i» ^ z> ,.- t . 



c'est-à-dire que la limite de —^ est zéro quand r tend vers 

zéro. Donc quand v tend vers zéro, la limite de r\ est la même que 
celle de rji tant que B est différent de zéro. Si B = o la limite 
est yio. Il en résulte qu'il y a une intégrale et une seule de l'équa- 
tion de Clairaut tendant vers zéro avec r : c'est celle qui corres- 
pond à B == o. 

Relation entre V aplatissement, la pesanteur et la force 
centrifuge à l'équateur, — A l'extérieur de la planète, comme 
cela a déjà été remarqué, on ne peut parler de surfaces d'égale 
densité, mais on peut encore parler .de surfaces de niveau. 

Pour ces surfaces, on a toujours l'équation 

AV=AV, + A(HY)=o; 

maïs AVq = — 4^p est nul à l'extérieur de la masse fluide. On a 
donc 

A(HY) =o 

A(HY) =y| H'^ +4 H' — ^î^^^^ h]== o. 
Faisons n=2, on a 

H''+— H'— -4-H = o. 

L'intégrale générale de cette équation est 
H^Ar2 + B/-3; 

PçiNCARÉ. — Figures d'équilibre. 6 
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inak à l'infini II doit être xiiil, il faut dexic<G|n»e A = o. On a donc 

II = B/-3; 

si M représente la masse fluide, V le volume enfermé par la sur- 
face de niveau considérée, on a 

^ — Y — /•' ' 
en posant C =- 



L'équation 



devient 



471 



D'où 



On a donc 



47; L O' ^ 6C J* 



/ »• 


' -"-■ ■4^- 6C 




i5 , 




iStoî» 




~ 8:tD • 


Comme on a posé 






re'=er, 


en a finalement 








BÉTEBUlNAnOif COMPLETE D£$ COEfFlÇlENT^ 8$ 

Celte relation établie pour Textérieur de la aia^se, ajit encore 
valable a la surface. On peut lui donner une autre forme. 

La force .centiîfuge à l'équateur est égale ii wV ; la pesanteur 

à la surface est — tt/D. Si on désigne par ïï), le rapport de lu 
o 

première a la seconde de ces quantités, on a 

2 

Or 71 est compris entre o et 3 ; il en résulte la double inégalité 

4 2 

Yj == 3 seulement si la planète est homogène. L'aplatissement 
ne peut atteindre la limite -^que si la planète est homogène. 



Détermination complète des coetûcients. — Nous allons 
maintenant démontrer que les équations 

AHY^-^'yP" HY 

qui déterminent les coefficients H du développement 



y=v,-h2^hy, 

n'admettent que la solution H=;o. 

En ejffet, soit n Tordre de la fonction sphérique, posons 

H = Z/"D. 
L'équation devient 

A(Z^DY)=iî^l±:i^Z,»Y; 
/■"Y est un polynôme sphérique ; en appliquant la formule 

A (PU)=ru" + 2:1±1. U' 1 p, 



84 MASSE FLUIDE HÉTÉROGÈNE. PROBLÈME DE CLAIR A VT 

on obtient finalement 

ce qui, développé, s'écrit : 

Z"D + aZ' (d' + .^:t_L d) + ^^^ ZD' = o . 
Si Ton fait le même changement de variables 



Z ' 

on a 

/•£' + e' + (2/l+l)£ + -^^(Ê + /l — i) = o. 

Cette cqiialion se discute de la même manière que l'équation de 

Clairaut ;- — -- — est compris entre les limites o et — o* 

On est donc amené à considérer les deux équations : 

/•Ê'+e^ + (2n+ i)£ = o, 

r£' + £'^ + (2/1 — 5)£ — 6 {n — i) =0. 

On construit les courbes en portant £ en abscisses et rs^ en 
ordonnées (fig. i3); on obtient deux paraboles, et comme pré- 
cédemment on démontre que e ne peut être compris qu'entre 
o et 3. 

A l'extérieur de la planète on a A (HY) =0, ce qui peut s'écrire 

H"+-H'-:^i^i±-^H=o. 
L'intégrale générale de cette équation est 

II devant s'annuler àTinfini, on aa = o, et par suite 

Z = ~ ^'-2(n+i) ^ 
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d'où 



Mais s doit être compris entre o et 3 ; donc b doit être nul. 




Fig. i3. 

H = o est donc la seule solution possible, à l'extérieur delà 
masse. Il en est donc de même à l'intérieur. 

Ce résultat ne subsiste pas si /i = i ou 2, mais nous avons 
déjà examiné ces cas. 

Détermination précise de V aplatissement» — Nous avons 
trouvé 



e= 



/•Y • 



Nous avons donc l'équation 



A(HY) = y izep'/Y, 
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et, en se rappelant la formule 

AC.HY) = y[h"+^ - -<^ h],. 

dans laquelle w = 2, on arrive pour déterminer H à l'équation 
linéaire avec second membre : 

H -I— _-_^„ep/. 

L'intégrale générale de l'équation sans second membre 
étant 

où a et j3 sont des constante» d'intégration, on déterminera l'in- 
tégrale de l'équation avec second membre, en regardant a et ^ 
comme des fonctions de r. La méthode de variation des cons- 
tantes donne 

(i) aV* + ^V-^ = o, 

par suite 

H'' = 2aV— 3pV-'*H-2a+ 12?/-' 
et, portant ces valeurs dans l'équation primitive, on a 

( 2) 2a'/- — 3 |3V-* = -^ ize^'r. 

Les équsltions (1) et (2) déterminent a' et ^' 



'?-^-§ V'-. 



On a donc 



Ja J b 
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Cette f»netvo& eomprend deux constantes arbîtiraire» a et h^ 
Co ai ment faut-il les choisir ? Soit r, le rayon de la planète ; pour r 
supérieur a Tj, p== o, f' = 0, par suite les intégrales se réduisent 
à des constantes 



/■ » /■■ 



Pour /•=oo, II doit s'annuler. Donc le coefficient de /• est 
nul; on a donc a=j\. 

De même pour r = 0,11 doit être fini;, donc le coefficient de 
/•"' doit être nul, si /• = o, on a donc i = o. 

Il est ainsi complètement déterminé. 

Portons cette Taleur de H dans l'équation 



■«-tK4^] 



Il vient 



•=^ / ^p'^^'-w / '^''"^'' 



3 w* 



Sit D * 



Que devient cette expression à la surface pour /• = r,?D Jevient 
D, la densité moyenne de la planète, et on a 






3 r\ 



(3) e = r^ I «?WrH--i-, 



ffl étant comme on l'a déjà indiqué le rapport de la force cen- 
trifuge à l'équateur, à la pesatvtéar.. 
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Moments d'inertie de l'ellipsoïde, — Calculons l'intégrale 
/ Y^rf<T, Y étant la fonction 







r' 




On 


a 








X = r\/i — [x^ cos «p, 








y = r\/l — ^K^sin ç, 








z = /'[X, 








d^ = dikd'-f. 




On 


a donc à calculer 




/ 


-^âic /»+i 


V fî 




/ 


jf' 


— 3ui.^ d\kd'f = 2v: j II- 


-6ix*+9jx» 


< 


^-i 


^-, 




ce qui 


donne 






(4) 




/^■^'=^- 




Calculons les 


moments d'inertie A, B, G de la masse 






k=f?{y'+z')d^. 





w 



On a 



Or on a 



AV = — 4Ttp, 

— 47tB =fx\{z' + x^)d':, 

— 4«C = fxY{x-' + yy': 

AV = AV, + A(HY) 



MOMENTS D'INERTIE DE L'ELLIPSOÏDE Sjf 

et de plus 

rA(HY)rVT = o; 

on effet, on a d'z = r^drd^d^ = drd's. 
Il vient 

C\Wii^d-z = Ç^{Wiydrd<t 

^ j iED_,.2rf,.,Yrfcr= / f{r)dr j Yrfa, 

ce qui est nul, car le second facteur est nul. On a donc. 

27r( A + B + C) = ÇliN^rUi 

=J^yydrd<j, 

et comme AV, ne dépend que de /•, on a 

27:(A + B + C) = 4« f^Y,r'dr. 



Or 



■-^.="^^-^-U^'^l' 



mais on a 

dY, 4^>D 

dr ~ 3 

donc 






K + B + C=^ I •'^lD,')d, 



SI on intègre par parties et que Ton remarque que A = B, on 
voit que Ton a 



aA + C== 
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En négligeant des intiniment petits on a A = C ; et par suite 



:-|l[,'D.-a /,-*D^,-]. 



Calculons maintenant C — A : on a 

47: (C — A) = fikY{z^—x^id-, 

4r:c-B)=yAV(s'-/)rfT. 

Donc 

8-(C— A) =— /'AVz-nWT 



x'-i-y'-t,-J 



en désignant par Y la fonction sphérique 
On a, comme on a vu précédemment : 

La première intégrale est nulle, car c'est le produit de deux 
intégrales : 

f\\ydrj\d:r. 
dont la seconde est nulle. L'intégrale 

= r^,.dr fr^d,; 



et comme on a 



y^v*= i|L, 
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il vient 
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817(0 



" j 



lô f • 


Donc on a 

C — A 



e^'r^dr 



qaantité positive. 

On a donc en se reportant à la formule (3) page 87 

. cp ^ 9 (C-A) 

On aurait pu faire une hypothèse diiTérente de celle de Clai- 
raut, par exemple, supposer un noyau solide non homogène, 
recouvert d'une masse fluide; la relation établie plus haut sub- 
sisterait encore. 



Comparaison de la théorie et des observations. — Posons 

C — A , 

— = J, on a 
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I 



est mesuré par la géodésie, on trouve — ^j-p ; 'f 

290,0 



est mesuré 



par la physique, on trouve— ^^-^ ; J est mesuré par la préces- 
sion des équinoxes, on trouve 0,00.32753 = ^ ^ on a donc 
Tégalité 

2 Çr'l^dr 
_______ jJl-=o,49i45. 

Revenons maintenant à l'équation de Clairaut 

rVD + ('/i^+ 57i)D + 2/-D'(i +r.) = o ; 



on a 



«'■ L J 2v/i+-/i 'fV/i+r, V>+^ 

On peut donc écrire Téquation de Clairaut 

5/*D 



t/i+ïi L 2 10 J 



J'intègre de o à /•, et j'ai 



I \/' + ^ L 2 'O J 

On peut appliquer le théorème de la moyenne, puisque 5/*^ D 
est une fonction positive ; ^ désignant une valeur comprise 
entre o et 3, et correspondant à une valeur de r comprise entre 
o et r^5 j'ai 



V^ + l J 



Ddr. 
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Voyons les variations de 









I 


+ 1- 


10 


K; 






\/H-S 




> 


lorsque 


$ varie 


de 


à3 . 


on trouve : 












l 


= 


K = 


I 








l- 


I 
~ 3 


K=: 


I jOooyS maximum 








V- 


I 
2 


K = 


I ,0002 










= 0,53 
= 0,544 

= 2 
= 3 


K = 

K = 
K = 
K = 


I 

= 0,99954 
■ 0,989 
0,9 

:0,8. 



On voit que K varie très peu et reste compris entre 1,00075, 
et o, 99954. K est donc très voisin de i ; et si dans la formule on 
remplace vi, par sa valeur 0,544? on trouve 



./,.. 



r-D, 



Drf/- 



= 0,497 -j^> 0,49663. 

Précédemment on a trouve pour la même valeur 

0,49145. 

Comment peut-on expliquer cette différence? 
On peut admettre qu'il y a des erreurs sur les quantités 
qui entrent dans le calcul ; c'est-à-dire des erreurs sur J et e,. 

La première équation donnerait alors comme valeur exacte : 



./,v 



Ddr 

= 049143 ^ -j- 8cj + ^'j/^ 8J 



= 0,49145 — 3o5,3i te^ + 155,26 8J. 
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La seconde équation donnerait : 



.liCr'Ddr 



2Kfr'Ddr . ^. 

J =0,49700 — i20,io6ej. 

Supposons 5J = o ; pour que ces formules soient comparables, 
il faudrait que Se^ soit égal à — 0,000027 et la valeur de e^ serait 

Si on suppose que l'erreur est commise sur J, on a 
oJ = o,oooo34- 

J'aurai alors pour valeur -tt . 

^ 002,18 

Cela est-il admissible ? Remarquons d'abord que le J qui 

Q J^ 

figure dans cette formule n'est qu'approximativement — -:; — . 
Cela revient à remplacer J par 



-X' 



2 
L'erreur commise est ûJ=-5- J^ : on trouve 

ôJ == 0^000006. 

Ce n'est pas de Ik que provient l'erreur. 

D'autre part les mesures de nutation sont assez précises. Si 
on désigne par N la constante de la nutation, on a 

N = — ^J [5,36542]; 

L 

en désignant par jjl le rapport ~ de la masse lunaire, à la masse 
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terresU^e ; le nombre entre parenthèse étant non le ooefiicîent 
maïs le logarithme de ce nombre. 

Soit p la préeession des éqïiinoxes : € est la «oanme de deux 
termes 

celui du à l'attraction de la lime 

p' = 866i35 — ^^ — J = 34,38 ; 

celui dû à l'attraction du soleil 

//' = 4871,05 J=i 5,95. 
Posons 



on aura 





{* 


^— ff • 


1 


+ i^ 


~— & • 


oN 


Se 


. ^^J 


N ■ 


s 


+ J ' 


4^ 


-^ 


■ +p" 


(S 

y Û2 


+ P 


SJ 
J ' 



P est mieux connu que N, on peut donc supposer S/> = o ; et 

SJ 
résolvant les deux équations où Ton regarde -j- comme égal 

o,oooo34 X 3o5,3i =o,oio3, on trouve : 

•oe 

= 0,010 

ûN 
-^=o,oo5. 

Cela revient à prendre pour N la valeur 9'' 17 au lieu de9''2i 

et pour UL la valeur 7; — r- au lieu de-^— - Ces chiffres paraissent 
* 01,5 02 

difficiles à admettre. 

On peut alors supposer que l'erreur provient de la somme des 

erreurs commises sur chacun des termes. 
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Une autre explication serait que Thypothèse de Clairaut n'est 
approchée qu'aux termes du second ordre près. M. Callandreau, 
qui a fait le calcul, à montré que considérer les termes du second 
ordre reviendrait a creuser l'ellipsoïde de 9 mètres seulement. 
C'est négligeable. 

On pourrait alors se demander jusqu'à quel point l'hypothèse 
de Clairaut serait exacte : par exemple, on peut supposer la terre 
complètement solide sauf la mer. Au moment de la solidifica- 
tion, les couches obéissaient à la loi de Clairaut, mais depuis ce 
moment, la vitesse de rotation aura pu varier (par exemple, les 
marées diminuent cette vitesse de rotation) ; et de plus l'aplatis- 
sement des ellipsoïdes internes aura pu se modifier par suite des 
tassements qui ont produit les continents et les montagnes. 

Enfin on pourrait supposer qu'une portion intérieure de la 
terre est restée fluide : les deux parties fluide et solide, agiraient 
alors séparément dans le phénomène de la précession : mais 
l'hypothèse de la fluidité intérieure de la terre n'est guère vrai- 
semblable. 

Jusqu'à présent nous avons considéré les planètes comme 
fluides et nous avons montré quelle devait être la distribution 
des densités à l'intérieur. 

Supposons maintenant que l'on ait une masse solide, de forme 
quelconque, où la distribution des densités est quelconque, mais 
entièrement recouverte d'une masse fluide de faible épaisseur, 
soumise à la pesanteur et à la force centrifuge due à la rotation. 

La surface extérieure sera une surface de niveau, sur laquelle 
on aura 

U = V+^(^^+i^2) = C^ 

mais à l'intérieur les surfaces équipotentielles ne sont pas néces- 
sairement les surfaces d'égale densité. 



CHAPITRE V 
MASSE SOLIDE RECOUVERTE D'UNE MASSE FLUIDE 



Théorème de Stokes. — Supposons que Von ait déterminé la 
forme de la surface extérieure 5, la masse M et la rotation w, 
on pourra déterminer en tout point extérieur la {>aleur du poten^ 
tiel V. 

Il y a une infinité de lois de distribution de. la matière dans le 
noyau, qui peuvent donner une même surface extérieure, mais 
quelle que soit cette distribution, le potentiel V, à Textérieur, est 
complètement déterminé par les éléments S, M et o). 

En effet la fonction V satisfait à l'extérieur de la masse h l'équa- 
tion AV = o. A l'infini V = o ; sur la surface on a' 



^(^*4-y) = u==K, 



et enfin l'égalité 



/ 



—— rfo- = — 4^M. 
dn 



Le problème de la détermination de V est complètement déter- 
miné. En effet, soit V une autre solution. Sur la surface, on a 

V — V=U — U' = K — K'== constante. 

On a d'ailleurs A (V — VO = o. 

PoiNCARÉ. — Figures d'équilibre. 7 
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Donc V — V est le potentiel d'une charge électrique répartie 
sur la surface, puisque c'est une fonction constante sur cette 
surface. 

La charge a pour expression 



/ 



dn 



l'intégrale étant étendue à la surface extérieure ; niais elle est 
nulle puisque l'on a par hypothèse 






-- — dT== — 47t;M, 
On 

an 



rf<7= — 47^M, 



V — V',le potentiel d'une charge nulle, est donc nul. C'est ce que 
Ton voulait démontrer. 

L'expérience du pendule qui peut faire connaître les valeurs 
de g à diverses altitudes ne peut donc rien nous apprendre sur la 
distribution des masses intérieures, que ne nous fassent connaître 
les éléments considérés dans le théorème de Stokes. De même, 
nous ne pouvons rien apprendre de plus par les perturbations 
produites sur le mouvement des autres astres. 

Développement de j^^^, en série de polynômes. — Soient 

O l'origine des coordonnées (fig. i4), M' un point attirant de 
coordonnées a;', y', ^', de densité p', M (:c, y, z) un point attiré. 
En prenant les coordonnées polaires on a : 

X = /-y/i — [jL~ cos cp, 
2/ = jyi — [x'^ sin ç, 
z = /-[x, 

x^z= r^^i — - [x'^ cos <p', 
ty = r'\/i — iL" sin'/, 




Fig. 14. \,,y = 7-'[x', 
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MM'2 = r^—2 rr' cos y + r" 
= x' + y' + z' — 2 [xx' + ^jy'+ zz') + x'' + y" + ;]^^ 
I I 



MM' 



. / 2 r' cos Y 7-'^ 



X X' . 7/ 7/' ^ ^' 

cosY = 7+-^ -^H 7* 



Développons suivant les puissances croissantes de ^— ; 

le coefficient de — ^^sera un polynôme n„ (^, y, z -, x' ^y\ 3^), 
On doit avoir les deux relations 






la première étant obtenue en considérant x^ ;y, z comme varia- 
bles, x'y'z\ comme fixes; c'est l'inverse dans la seconde. 



V MM' / ~ 



^' \MM7 



52 _i_. 52 ï 



MM' MM' . MM' 



d^2 ^ ^y' ' 0^^ 



MM' MM' MM' 



On sait que l'on peut écrire 



==.+p.^+... + p„(4)V. 



où P„ est un polynôme en cos y de la forme 
A, cos" Y + Aj cos"" 'y -)-... 
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si on y remplace cos y par sa valeur, on aura 

Le polynôme Q sera un polynôme homogène et de degré n 
par rapport à chacun des systèmes de variables [x^ y, z^ /), 
{x'^ y' , z', 1^). II ne contiendra d'ailleurs ;• et z*' qu'à des puissances 
paires, on pourra donc remplacer r^ et i^^ par leurs valeurs 

^' + y' + ^' et x" + y" + z"', 

on aura donc : 

_^ I _Y ^n[x.y.z;x',y\z') 

I -. ^ -Zj 

MM' ~Zj /•'"■'' ' 

Q' représente un polynôme homogène et de degré n par rapport 
h chacun des systèmes de variables x^y^z) x\ y\ z' , chacun des 
deux ensembles y entrant de la même façon que l'autre. 
On a d'ailleurs 



Donc on a 



A'-Q4 



et, comme /• ne dépend pas de x'^ y' ^ z' , on a : 

A'(Q')=o. 

Q' est donc un polynôme sphérîque d'ordre n par rapport aux 
variables x',y'^ z' et par suite par rapport aux variables [xy y, z). 
Nous savons qu'il y a 2/>+ i polynômes sphériques d'ordre/;, 
de forme 

^pq — l e \^ V- ) t/(JL^ + 9 ' 



CAS OU R EST SUFFISAMMENT GRAND loi 

OÙ q prend les valeurs 

—Py—P+^9 — i^o^i» 2, . . . p—î,p. 

On peut donc écrire 

formule où P est un polynôme fondamental en .r, "y, 3, Q un 
polynôme en [x', î/\ z') ; Q' est d'ailleurs fonction de 

cosY= ul[jl'-|-\/i — jjl"^ \/i — ^'^ cos (cp — cp') 

Q' ne dépend donc que de la difTérence (cp — cp') : il en résulte 
que le coefficient de P^^ne peut être que 

et on pourra écrire : 

MM' ~ Zj r'^'^' 

A^ç est un coefficient numérique que Ton peut déterminer faci- 
lement. 

Application. — Supposons que /* soit suffisamment grand par 
rapport au rayon de la planète. On pourra développer le poten- 
tiel suivant les puissances croissantes de et on aura 

'^ r 

V V^ ^pq "pq 

Zj r'i''^' ' 



D'après ce que Ton a démontré plus haut, la connaissance de 
S,Meta> suffit pour déterminer V, donc ces éléments permet- 
tent de déterminer B^, ^ ; or, on a 

/ MM' Zj r'^^' I P^ ''-'^' 
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il vient donc 



B,, = A,, / p'PV,_,rfT' 



P'p^ g est l'un des polynômes sphériques fondamentaux d'ordre p. 
On doit donc savoir calculer 



/ p'p;,'^-'- 



Considérons d'abord le polynôme d'ordre o : P, , = i ; on a 



/' 



Les polynômes d'ordre i sont ^', y', z' : on connaît donc la valeur 
des intégrales 

Ce sont, à un facteur près, les coordonnées du centre de gra- 
vité de la masse. Supposons que la surface extérieure possède 
un centre de symétrie, et prenons l'origine des coordonnées en 
ce centre ; les intégrales sont nulles puisque ce sont les coeffi- 
cients des termes en -j-, -^, —3-, dans le développement du 

potentiel et que les surfaces équipotentielles doivent admettre 
l'origine comme centre de symétrie. Le centre de gravité de la 
planète coïncide avec le centre de symétrie. 
Les polynômes de second ordre sont 

x'y\ y'^'-, z'^, 

z'^ — j^\ z'^ — y'\ 

Si / x' y' ^'dz est nulle, ainsi que les deux autres intégrales 
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analogues, les axes de coordonnées sont les axes principaux 
d'inertie du système. Si / [x'^ — y'^) pd'z est connu, on connaî- 
tra la différence entre les moments principaux d'inertie. 

On a trouvé une relation entre e^^ '\> Dj et w ; les deux premiers 
éléments sont donnés par la connaissance de S; D^ est connu 
puisque Ton connaît M et S ; o est connu puisque l'on connaît w ; 
C — A est aussi connu. On doit en conclure que la relation déduite 
de l'équation de Clairaut est vraie, quand même la distribution 
intérieure de la matière ne serait pas conforme à l'hypothèse 
faite pour l'établir. 

Relation entre F aplatissement, la pesanteur et la force 
centrifuge à Féquatear. — Proposons-nous de calculer les pre- 
miers termes du développement de ^„^ ' 



On a d'après la formule du binôme 

. — JLI . L l f'*' 

MM' 



ïM^=-L'-Tr— ^"^•^+— ) 

négligeant les termes de l'ordre de —j-\ nous avons 

I I . r' cos Y , î r r'^ , 3 ,2 2 1 



MM' 



^p+,,,+.._(f!±^!±fl] 



^±^{xx'-\.yy'+zz'y 



= — + -3- {xx' + y y' + zz') 
-Ar[{^' + f-{-z'){x'''+-y"+z'')-i[xx'+yy' + zzy]. 

Le coefficient du terme — r- peut s'écrire 

— 6 [xx' y y' -\-xx' zz' + y y' zz'), 



io4 
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On a alors 



-/ 



mmF 



Si M est la masse, a, ,3, y les coordonnées du centre de gravité, 
A, B, C les axes principaux d'inertie, on a 

V= i^ + J^f+Ëf+IL H__L^[3(.^_,^) (B-A) 

Supposons que la surface soit de révolution, que Torigine des 
coordonnées soit placée au centre, que les axes de coordonnées 
soient les axes principaux d'inertie : alors 



B = A 



et 



Calculons maintenant la valeur que Ton déduit pour V de la 
forme extérieure de la surface, que nous supposons un ellipsoïde 
d'aplatissement e^. Ce potentiel est de la forme 

V=-^+IIY + 2QZ, 
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OÙ Y est la fonction 

et Z est une fonction sphérique différente de Y. 

H et Q doivent être d'une forme particulière, puisque le poten- 
tiel doit s'annuler à l'infini, nous avons déjà vu que Ton avait 

11 = A 



IIq et q^ sont des constantes; on a donc 



A la surface libre on a 






Il vient donc 

M w^-^ , w^-^Y , H„Y , v' Q„Z 



,.3 T 7 , ^n+1 ^ • 



Posons 

On a vu (p. 69) que l'on avait 



_3^ 



de plus on a 



^-.[-^]- 



• 
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Tous les termes étant petits, sauf le premier, on peut y négli- 
ger ej et on a 

L'équation U = O" se réduit donc à 
M Me. Y , (o'/f i^hjY H, Y Sy,Z _ , 

Comme par hypothèse nous avons un ellipsoïde : q^ = o; donc 

V se réduit à 

M ^ H,Y 



et on conclut, en comparant les deux valeurs de V, que 

On a donc l'équation 

M /Me. C--A\ , ti'r] wVfY 



__ /Me. C— A\ 6)V; 



i--^V^=c'«-; 



6 

on en conclut que le coefficient de Y est nul, et que l'on a 
Me, _ C— A 





ir^ 2/'î ' 6 . 


On a d'ailleurs 






M = i-uD.rî, 




M 




M 


Il vient donc 






a/1 »' 3 6 ' 




C-A_2;-,^-„/ <p\ 
rî 3 r a^' 
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et, en supposant le rayon de la terre égal a l'unité, on a 



C-A=-^(n-i-). 



Pesanteur en un point de la surface. — Les composantes 

ôU ÔU ôU 
de la pesanteur sont -^ — 7""â~" » "T" > c* oï^ î> 



^=vW^(f7^(¥) 



La composante suivant le rayon OM est, s étant le cosinus de 
l'angle de la verticale avec le rayon OM, 

dU 



comme e est petit, g cos e est sensiblement égal à g, aux quan- 
tités du second ordre près, on a donc 

dU 



e "" . dr ' 
M ^ ^ (C — A) Y 2toV (oVY 

2 r* 3 3 



^ = —+3 „ ,.^ 



Si maintenant on fait r = i -j- Ç, 



-JX= ï —2^' 



on trouve 

nr M oMr_I ^^ ^^^ ^ 

2 O 



^ = M-2M^ + — (C — A)Y— ^-a>% 



3 3 



g=M{i-^q+g,{e,-±yY 

_£__,_iLL ^y Y I ^'f 

g,-' 3 6 ' -*" 3 
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C'est la formule de Clairaut qui donne la pesanteur en un 
point, en fonction de la latitude de ce point. 

Inûnence de F altitude* — Si on s'élève en ballon , la cor- 
rection est facile h faire, on n'a, dans la formule, qu'a rem- 
placer r par i + ^î + /^ au lieu de le remplacer par i + JJ, maïs 
ce n'est pas dans ce cas que l'on a à faire la correction d'alti- 
tude : le plus fréquemment cette correction se présente pour un 

lieu situé sur un plateau, et il faut 
tenir compte de l'attraction due à 
ce plateau. 

11 faudra donc tenir compte des 
aspérités positives ou négatives; il 
faudra d'abord choisir une surface 
définissant le niveau général du 
plateau : cette surface sera la sphère 
décrite sur OM comme diamètre 
(fig. i5), O étant le centre du sphé- 
roïde, M le lieu de l'observation. 
Figr. i5. On calcule l'attraction due au vo- 

lume compris entre les deux sphères, 
l'erreur est peu considérable et d'autant moins que l'attraction 
des montagnes est très faible. On prend, en général, comme 
densité de la masse la densité moyenne des roches a la surface, 

c'est-à-dire Dp en représentant par D^ la densité moyenne 

du sphéroïde. 

Correction de Bouguer. — Supposons une couche sphérique 
homogène d'épaisseur h et de densité p ; le volume de cette 
couche sera 

et sa masse 

47î/t/-?p. 

La correction apportée à la mesure de g sera ^iz^hrl et on a 
hg 4^p/i 3p/t 

't^ 4 ,^,. ^IVT' 



CORRECTION DE BOUGUER 
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D'un point M situé à une certaine distance au-dessus de la 
couche, menons le cône tangent à la sphère (fig. 16) ; le cercle 
de contact partage la couche en deux portions, et nous savons 
que les attractions de chacune de 
ces deux portions sur le point M 
sont égales. 

Appliquons ici ce théorème. 
Nous avons supposé que la surface 
de la planète était limitée à la 
sphère décrite sur OM comme 
diamètre ; si donc on suppose 
que l'attraction de la masse MRS 
se compose de l'attraction d'un 
certain nombre de couches homo- 
gènes et minces, la partie de ces 
couches limitée à la sphère (OM) 
aura sur M la même action que la 
moitié de la couche totale. Il faut donc calculer l'attraction sur M 
d'une couche entourant le sphéroïde de densité p et d'épaisseur h 
égale a l'altitude du point M. La correction à apporter a la 
pesanteur sera la moitié de cette attraction. 

L'épaisseur h de la couche peut se représenter par un dévelop 
pement de la forme 

le potentiel dû à cette couche, d'après une formule connue sera 

VI 47uAY 




Fig. 16. 



j(2/i+i)r" + * ' 



l'accroissement de g sera 
ôSV 



^ ô/' Zj 2/1 + I ' 

car on fait ensuite /• = i dans la formule, et on a 

D7"^2j 2/z+i ^- 



s- _ 3p 



00 
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Si on se borne au terme /i = o 

og _ 3p . 

go~ ^r 

si on en prend la moitié, c'est la correction de Bouguer. 

La correction totale comprend donc deux termes : i** celui 
qui provient de l'augmentation d'altitude sans tenir compte du 
plateau; 2^ celui où Ton tient compte de la masse du plateau, on 
a finalement la formule 

80^ ^h 3 p/i 

Si on applique cette formule on trouve une valeur trop faible : 
la valeur fournie par l'observation est, en général, sur les conti- 
nents, inférieure à la valeur calculée ; elle lui est au contraire 
supérieure sur les océans et dans les îles isolées. 

M. Faye a ainsi trouvé qu'il vaut mieux ne pas tenir compte 
de la correction de Bouguer, ce qui correspond a des masses qui 
devraient se trouver sous les océans, et à des vides qui existe- 
raient dans les continents. 

A la surface des mers on a U = U^. La surface U = Uy est une 
surface qui, se prolongeant sous les continents, doit donner une 
surface de révolution, mais il n'en est pas ainsi. Cette surface 
s'appelle le géoïde. 

Pesanteur en un point de la surface du géoïde. — Considé- 
rons une sphère de rayon voisin du rayon de la terre ; soit t^ l'al- 
titude du continent au-dessus de cette sphère, on peut écrire 

ÎJ' la distance à l'ellipsoïde, peut s'écrire 

^ — IÇf est la distance de la sphère a l'ellipsoïde au lieu consi- 
déré ; 
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Soit h l'altitude du point au-dessus du géoïde, 

/i=V/cY 
Ç' — h est la surélévation du géoïde au-dessus de Tellipsoïde : 



on a 



^-h^'^^ll-q-k)^, 



Î^ — 'Q &e réduit à deux termes : celui qui correspond à Y = i , 

et celui qui correspond à Y = ^^ . 

A la surface du géoïde U = ^^ = constante ; en un point d'al- 
titude A, on a 

On a d'ailleurs 

dU 



d/- *»• 
Donc 

Les coefficients Y sont connus : ils correspondent au développe- 
ment de 

si on rétablit l'homogénéité, on a 

V _ I S(H — ^ — y)Y 

et comme on a : 

ôV 
^ = -17' 
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on en conclut : 

g I s(;»+.)(H-^-y)Y 

• Les opérations géodésiques font connaître A, les opérations 
de nivellement font connaître !J, les observations pendulaires 
font connaître g. Théoriquement, deux séries d'observations suf- 
firaient pour connaître ÎJ, h et g; mais pratiquement les observa- 
tions d'aucune des trois espèces ne sont assez nombreuses, ni 
assez précises pour qu'on puisse se passer d'une série d'observa- 
tions. ^ 

Jusqu'à présent l'on peut dire que l'on a trouvé la surface du 
géoïde peu différente de celle d'un ellipsoïde de révolution. 

Nous allons, dans ce qui suit, montrer qu'il peut exister des 
figures d'équilibre autres que l'ellipsoïde de révolution. 



CHAPITRE VI 
FONCTIONS DE LAMÉ 



Coordonnées elliptiques, — Considérons les quadriquesliomo- 
focales 



K' — a^ 

Par un point de l'espace passent trois de ces surfaces : en effet, 
.r, y, z étant donnés, on a pour déterminer '}} une équation du 
troisième degré dont les trois racines sont séparées par les 
nombres a^, i^, à. 

Il en résulte que Tune des racines, la plus grande, supérieure 
à a^y correspond à un ellipsoïde; la racine moyenne, comprise 
entre a^ et t^, définit un hyperboloïde à une nappe; la troisième, 
comprise entre P et c^, définit un hyperboloïde à deux nappes. 

Soit p^, [JL^, v" ces racines, on a donc 

p2>a2>[JL2>i2>v2>c2. 

Réciproquement si o, »jl, y sont donnés, on a trois surfaces se 
coupant en huit points, placés symétriquement par rapport aux 
divers plans de coordonnées. 

Si on ne considère que les points situés dans un trièdre trirec- 
tangie donné, les trois nombres p, a, v définissent un point et un 
seul. 

On peut écrire 

^' I f I ^' ,._ -(X^-p')().^-ix^)().^-/^) 

\^ —a^~^ X' —0''^ X' — c' ~Çi:' — a') [V —ù') (V — t-j' 

PoiNCARÉ. — Figures d'équilibre. 8 
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Si on multiplie par \^ — a% et qu'on fasse ensuite X^ = a% 
on a 

[a^ — U') [a' — c') 
on obtiendra donc 

et de même pour y et z. 



V [c' — à')[à — b-') 

Le signe qu'il faut prendre devant le radical correspond au 
Irièdre des coordonnées que l'on considère. 
So\if{x) un polynôme, et posons 

N=/-(v^), 

RMN sera un polynôme symétrique des trois quantités p^, jjl^, v'^ 
racines de Téquation 

r2 y^ ^2 



et par suite une fonction entière de x^^ y^, z^. 
Si on pose 

ou bien 
OU encore 
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et que Ton prenne pour M et N les mêmes fonctions de |jl et 
de V, RMN sera un polynôme en ^^, ?/^, ^^, multiplié par^, par^ry 
ou par xyz. 

Supposons que R soit un polynôme entier en p^ ; on peut le 
décomposer en facteurs du premier degré et écrire 

R=n(p^-À?). 

Alors 

RMN = n (p2 — If) ([a'' — 1\) (v^ — X«) 

RMN sera donc, dans tous les cas possibles, un produit de fac- 
teurs de la forme (i) multiplié par un certain nombre de quan- 
tités X, y^ z. 

Élément linéaire en coordonnées elliptiques. — Les sur- 
faces p =C*% [jL=C^% V = C*^ forment un système triple ortho- 
gonal et l'élément linéaire peut s'écrire : 

a, [3, Y étant des fonctions de p, [j., v que nous allons déterminer; 
on a 

ds'=dx'^ + dy^ + dz\ 

dx p<fp ]x.d)x vrfv 

~"^ f — o?- "*" ^x' — o^ "*" ^^ — à' ' 

dy pdp ikd^k v<rZv 

■^"^ p^—ù^ "^ [X^ — b' "^ v^ — b' ' 

dz prfp ^ ikdiL ^ vrfv 



z p* — c^ [x'* 



On a donc 



j , jT x' y^ , z^ ~\ 

[^2 «S ^2 -1 
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et ensuite 



Je vais poser 

Q» = ((x^_p^)(p^_v»)(v'-i.=), 



P-^ 



B^ = 



(i.^-«')(l.'-^^)(t.^-c-') 



.i 



^, (y' — a')(v'' — ^')(v' — c') 

^ - ;5 . 

A^etC sont positifs, B^ est négatif ; donc A et C sont réels, 
B est purement imaginaire. 
Avec ces notations on a 



Q 



P= 



.= Q_. 

y/v^ — p*est purement imaginaire, donc fi est réel. 

Equation de Laplace. — Voyons ce que devient Téquation 
de Laplace en coordonnées elliptiques. Nous avons déjà établi 
(p. 4o)> 4"^ si on désigne Télément linéaire par 

ds' = a^cfp^ + [ihli»} + fdy\ 
on a 

ô r^Y ôv 



a^Y Zj àp La dp J 
L'équation AV = o devient donc : 

Zj dp \ a dp / ^' ' 
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On a 
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iL _ AQ v/ZZZEZZZ 

BC V(,._p.)(p._j,.) 

_ Av/^^(îx'— V») 
BC ^' 

et Téquation de Laplace est alors : 

A seul dépend de p. On a donc 

en multipliant par ABC on obtient l'équation 

Les coefficients de cette équation sont réels, car A^, B^, C^ sont 
réels, et par suite leurs dérivées le sont aussi. 

Fonctions de Lamé. — Parmi les solutions de l'équation de 
l.aplace, il y a des polynômes de la forme RMN : on les appelle 
fonctions de Lamé. 

Remarquons d'abord que nous avons les identités évidentes : 

2(l.^-v')==o; 
Il en résulte que l'équation (2) est équivalente à la suivante 
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Supposons que Ton aitV= RMN et que R satisfasse a Téqua- 
tîon 

M et N satisferont a des équations analogues, et la fonction V 
satisfera à Téquation de Laplace. Cette condition suffisante est 
aussi nécessaire, c'est-à-dire que si V est de la forme RMN et 
satisfait à Téquation de Laplace, R satisfait à l'équation (3), et par 
suite M et N satisfont aux équations analogues obtenues en rem- 
plaçant p par \k et v, et A par B et C. 

En effet si on fait u=v l'équation doit être encore vérifiée; 
R se réduit a M en ce cas, et Téquation se réduit à 

Existe-t-il de telles fonctions ? Considérons la fonction p [u) 
définie par l'équation différentielle 

p'u = 2 ^{pu—e,){pn — e,) [pu — e,),^ 
où 

Posons 





p2=/iW + /l, 




a^= e^ -h A, 




P= ^2 + /*, 




c'=^e,-[-h; 


on en conclut 








On a donc 





:2Ap. 



On a d'ailleurs 



, du 



donc 



et par suite 



L'équation 



devient 
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du 



. rfR rfR dp rfR 



d^ do du du 

i^ + (H?^+K)R = o. 



rf/^^ 



Peut-on choisir II et K de manière que R soit un polynôme 
en p, ou un tel polynôme multiplié soit par yp^ — a^, soit par 

\/f — b\ soit par y/p' — c'? 

Ce polynôme, si on y remplace p^ par p (w), sera une fonction 
doublement périodique admettant les deux périodes de p[u) 
2a)j et 20)^. Elle admettra une double infinité de pôles qui seront 
les pôles de p (?/), c'est-à-dire u == a/ww^ + 2 /ico^, m et n étant des 
entiers positifs, négatifs ou nuls. Ces pôles qui sont des pôles 
doubles de p (u) seront pour le polynôme des pôles d'ordre 2n. 

De plus y/p^ — a^ =ypu — e^, est une fonction doublement 
périodique admettant les périodes 20)^ et 4^^ et ayant comme pôles 
simples les pôles de p {u); donc si on ajoute à u la valeur sw^, 20)^ 
ou 20)3, la fonction R garde sa valeur ou change de signe sui- 
vant qu'elle contient o ou 2 radicaux, ou bien qu'elle est le pro- 
duit d'un polynôme par i ou 3 radicaux. 

Donc R est une fonction paire, ou une fonction impaire, mais 
dans un cas comme dans l'autre on doit avoir le développement 



si on porte cette valeur dans l'équation, en se rappelant que 
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on obtient l'équation : 

r n(n-\-i)a, (n — i){n—2)a, l 

L II"*' u" "^"J 

Pour que cette équation soit vérifiée il faut d'abord que le 
terme en ^ disparaisse, ce qui n'est possible que si Ton a 

\l = — n{n 4- i). 

On a d'ailleurs si q est le degré de R, et qu'il faille lui ajouter p 
radicaux, (p = i, 2 ou 3) : 

n = 2(jf -j-p- 

Formation des fonctions de Lamé. — Ainsi donc les poly- 
nômes de Lamé, doivent s'ils existent satisfaire aux conditions 
suivantes. 

Regardés comme fonctions de x, t/, z ce sont des produits de 
facteurs quadratiques multipliés par un ou plusieurs des facteurs x^ 
t/^z. Si on remplace j;, 7/,^ par leur valeur en fonction de p^,[JL^, v^, on 
obtient une fonction de la forme RMN, où R désigne, soit un poly- 
nôme en p^, soit un tel polynôme multiplié par 1,2 ou 3 des 
radicaux 

v/p=:^ v/?^^^, V7^^> 

et M, N désignent les mêmes fonctions ou est remplacé par [jl 
et V. 

On a pour Q huit formes possibles, à chacune desquelles cor- 
respond une forme pour R. Si n est le degré de a», le degré de R 
en p^ sera suivant le cas 

n n — I n — 2 n — 3 

— ou 



22 2 2 



suivant que R contiendra o, i, 2 ou 3 radicaux. 

Enfin si dans R on remplace p^par /^w, on obtient une fonction 
qui, lorsqu'on remplace u par u -f- 20)^, change son signe ou reste 



J 
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invariable, suivant qu'il y a un nombre impair ou pair des deux 
radicaux y/p- — ^^etyp^ — c-*, car \/pu — e^ ne change pas de 
signe pour w + aw^. 

De même quand on augmente u de aw^, y/^" — e^Qi\pn — e^ 
changent de signe et par suite R changera de signe ou non, sui- 
vant qu'il contiendra un nombre impair ou pair de ces radi- 
caux. 

Ces résultats sont résumés dans le tableau suivant, où la pre- 
mière colonne contient la forme de Q, la seconde la forme corres- 
pondante de R, la troisième le degré du polynôme qui entre 
dans R en fonction du degré de Q, la quatrième et la cinquième 
indiquant si la fonction change désigne pour l'adjonction à m des 
périodes 20)^ et aw^; Q^ désigne uniquement un produit de facteurs 
quadratiques. 



Q.xyz. 



R ne contient pas de radical 
RconlicntV^p- — a^ 

» y/ p'-^ — />^ 



y /p2_c2^p2_ql 
^p2_rt2^p2_/,2 



\ 



n— 3 



+ 

+ 

+ 
+ 



+ 



/i+2 



+ ' - 
-) " 



Détermination de K. — Ceci posé, cherchons à déterminer K. 
Pour les fonctions du premier genre, R est un polynôme 

d'ordre contenant h i coefficients homogènes. 

Substituons dans l'équation différentielle, nous avons un poly- 
nôme de degré — •+ I, qui devra être identiquement nul. Le 
coefficient du terme de degré le plus élevé disparaît de lui-même 
a cause de l'hypothèse II = — /e (/i+ i). Il reste alors — + 1 
coefficients à faire disparaître. 

Les 1- I coefficients du polynôme doivent donc satisfaire 
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à f- I équations linéaires et homogènes, dont les coefficients 



dépendent linéairement de K; l'élimination des coefficients du 
polynôme donnera donc pour K une équation de degré 1- i . 

Si R est du second type, on aura à déterminer coeffi- 
cients d'un polynôme de degré . Le degré de l'équation 

en K sera donc, d'après un raisonnement analogue au précédent, 
/i -}- I 



Pour ceux du troisième, l'équation en K sera de degré — . 
Pour ceux du quatrième, l'équation en K sera de degré 



2 

Ces résultats figurent dans la sixième colonne du tableau pré- 
cédent. 

Il résulte de ceci qu'à chacune des racines des équations en K 
ainsi obtenues correspond un polynôme et un seul. 

Si n est pair, il n'y a de polynômes que de première et de 
troisième espèce ; leur nombre est 

2 2 

Si n est impair il n'y a de polynôme que de seconde et quatrième 
espèce ; leur nombre total 

n — I n-]-i 

h 3 ■ = 2/1 + I . 



Donc il y a au plus in-{- i polynômes de Lamé. Je dis au plus, 
car les équations en K pourraient avoir des racines multiples ; 
nous verrons qu'il n'en est rien. 

Analogie avec les fonctions spbériques. — Supposons que 
j:^ y y z soient infinis du premier ordre ; [jl et v restent essen- 
tiellement finis, mais p deviendra infini du premier ordre aussi. 

L'ensemble homogène des termes de degré le plus élevé dans 
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la fonction de Lamé, satisfait évidemment a l'équation de Laplace. 
C'est donc un polynôme sphérique que Ton peut écrire \\{x,y^z). 
Le terme de degré le plus élevé dansRMN en p sera ap^'MN. 
Passons aux coordonnées polaires 

( x = r sin 8 cos cp, 
< î/ = r sin 6 sin cp, 
V z = /'cosQ, 

si p est très grand il diffère peu de r; u, v dépendront au contraire 
de cp et de 8. 

En effet l'équation qui détermine [x peut s'écrire 



Si Xyij^z sont très grands, l'unité est négligeable; il reste une 

équation du second degré; [i. et vies deux racines non infiniment 

Il z 
grandes ne dépendent que du rapport-^ et — c'est-à-dire de cp 

et de 6. 

On peut alors écrire 

r"n (sin 6 cos cp, sin 8 sin cp, cos 8) = ap"MN, 
et comme /• ==p, on a 

n (sin 8 cos cp, sin 6 sin cp, cos 8) = aMN ; 

le produit MN apparaît donc comme une fonction sphérique, 
et n sera décomposable en un produit de facteurs quadratiques (*). 
M. Moutard a démontré que c'était le seul cas où un polynôme 
sphérique était décomposable en un produit de facteurs quadra- 
tiques et de facteurs linéaires. 

MN, considérée comme fonction de 8 et de co, est donc une fonc- 
tion sphérique : il en résulte que toute fonction de [x et v peut se 
mettre sous la forme d'une somme 

y.aMN, 



(*) Voir Journal de Liouville, t. XI, 1846, p. 278. 
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puisque toute fonction peut se mettre sous la forme d'une somme 
de fonctions sphériques ; mais l'analogie des fonctions sphériques 
et des fonctions de Lamé, va apparaître davantage quand on étu- 
diera le cas où les surfaces homofocales sont de révolution. 

Cas des ellipsoïdes de révolution. — Cela peut se présenter 
de deux manières différentes suivant que Ton a a = b, ou b = c. 

1° Si rt = b, on a des ellipsoïdes de révolution aplatis quand 
Â- est supérieur à a^ et des hyperboloïdes h deux nappes si 
a^ <\^ <c^. Les hyperboloïdes à une nappe se réduisent à des 
plans passant par Taxe 0^. 

2° Si b = c, on a des ellipsoïdes de révolutions allongés autour 
de O.r, pour X^ > a^. Si a" > Â^ > i^ on a des hyperboloïdes à une 
nappe, les hyperboloïdes a deux nappes se réduisant à des 
plans passant par l'axe Ox, 

Dans ces deux cas, soit jjl, soit v ne peut servir de paramètre, on 
est obligé à,e prendre comme paramètre l'angle du plan avec un 
plan fixe, le plan xOz par exemple dans les deux cas. Nous 
avons : 

X V (p:i _ ,j2^ 1^,^^ _ ^^) («2 _ ,^-2^ V ^. _ p. • 

Si Z> tend vers a, jx tend aussi vers a, et l'on voit que l'on a : 






— b' 



a" — [JL^ 



Si cp est l'angle du plan passant par le point (:r, y, z) et Taxe Or. 
avec le plan œOzy on voit que [jl peut être considéré comme fonc- 
tion de cp, et l'on a 



\)? = a^ cos^ îp -4- b^ sin^ cp. 



On voit que l'on a 



limV ^,_^^.. =sincp, 
limy "Jt =cosy. 



a' — b' 
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Dans les mêmes conditions Téquation qui détermine v peut 
s'écrire en supprimant l'unité, négligeable devante, y etz. 



Si on passe aux coordonnées polaires, on a 
/•^ sin^" 



+ 



et ensuite 



- «=v/^ 



sin 



cos 9 : 



V a'—c' 



Dans ce cas on peut immédiatement former les fonctions de 
Lamé. 

Nous savons que MN est une fonction sphérique ; M ne dépend 
que de o, N ne dépend que de ; MN est donc l'une des fonctions 
sphériques fondamentales considérées plus haut (p. 42). 

On peut alors écrire M = cos p'f ou sin jy^ ; et pour N on 
prendra 



en posant •. 



F,Ucose)==F;:[v/^] 






On voit que les zéros de N se déterminent en calculant les 
zéros de Fj, abstraction faite des zéros ±: i . Ces zéros sont tous 
réels et distincts ; car, comme on le voit par le théorème de 
RoUe, les zéros de Ff^ séparent les zéros de F^~K On a 



(-'■)--v/^- 



Donc si p est pair, F{J ne changera pas de signe avec sin 9; au 
contraire sip est impair, F„^ changera de signe avec sin 9. 

Si on remplace ^, 1/, z par leurs valeurs en coordonnées polaires, 



ii6 
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Q changera de signe s'il contient l'un des facteurs x ou y \ 
mais s'il n'en contient aucun, ou s'il en contient deux, il ne chan- 
gera pas de signe; p sera donc impair dans le premier cas, et 
pair dans le second. 

On voit aussi que y et sin p cp changent de signe en même temps; 
donc si Q contient le facteur y, M devra être pris égal à sin/^ï), 
sinon M devra être pris égal a cosp'f. 

De même si b' = c^, on prendra comme coordonnées polaires 

X = rcos 0, 

y = r sin cos cp, 

z= /• sin 6 sin "f. 

La fonction N sera siny^'^ ou cos prs, la fonction M sera 

Fi;(cos6). 

On prendra y; impair si Q contient l'un des facteurs y ou z; 
p sera pair au contraire si Q ne contient aucun de ces facteurs ou 
s'il les contient tous deux. On prendra pour N, sin p (f ou cosy^ cp 
suivant que Q contient ou ne contient pas le facteur z. 

Résumé. — Le résultat de cette discussion est résumé dans le 
tableau suivant : la première colonne contient les 8 formes possi- 



P. . . . 

Pxx. . 

Pxj. 

Fxz. 

Pxv 

^xzx 

Pxxj 

?xxjrz 



K 

Rv/p=* — a-î 

RV/p^ — c^ 

Rt/(p=î — c^)(p-2 — a^j 
Rj/(p*_aV(p2-6*) 



R vV--«^Jp^^^'W-^ 



a- 


= h' 


b^ 


= C2 


o 


Cos 


o 


Cos 


I 


Cos 


o 


Cos 


I 


Sin 


I 


Cos 





Cos 


I 


Sin 


I 


Sin 


o 


Sin 


I 


Cos 


I 


Sin 


o 


Sin 


I . 


Cos 


o 


Sin 


o 


Sin 
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bles de Q, la seconde contient les formes correspondantes de R, la 
troisième et la quatrième contiennent les cas limites pour a^ = b^; 
dans la troisième on a mis le signe o si p est pair, i si p est 
impair ; la quatrième indique si M se réduit à cos /?cp ou sin p'f ; 
la cinquième et la sixième contiennent les indications corres- 
pondantes dans le cas où b'^ = c^. 

Les valeurs de K sont réelles. — Nous allons maintenant 
démontrer que les équations en K ont toutes leurs racines réelles. 

Supposons d'abord a^ = h^. On sait alors trouver les valeurs 
de K. Nous pouvons poser [ji^ = a^ cos^ îp + b^ sin^ cp, puisque [a 
est compris entre a et b et lorsque a = b^ ceci est encore la 
valeur de M ; Téquation qui définit [x 

peut s'écrire, en prenant comme variable cp 

\^ — ^) , 2 + \y — ^j. s"^ ? ^os cp - 



_r,i(«_|-i)a2_KlM =0 

et si on fait a^ = b^ on trouve l'équation 

(^ — ^) -^rr — /2(«H-iK — k|m==o. 

Cette équation doit avoir comme solutions cos/?cp et sin/?'f, on 
a donc 

K — n{n-^\) 2 
^2 ..2 =7^ • 



L'équation en K a donc n -(- i racines réelles, correspondant 
aux valeurs de /?, 

o, i, 2, . . . (ti — i),«. 

On voit que les équations en K ont toutes leurs racines réelles ; 
les valeurs paires de/;, séparent les valeurs impaires. 
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Lorsque b variera de a vers c les racines des équations corres- 
pondant aux deux groupes de fonctions de Lamé, contenant 
zéro et deux radicaux d'une part, sépareront les racines corres- 
pondant aux deux autres groupes contenant un ou trois radicaux ; 
cela ne cessera d'être, que si une des racines d'un groupe d'équa- 
tions, devenait égale à une racine îp d'un autre groupe. 

Soient, K cetle racine, R et S les fonctions de Lamé corres- 
pondantes, appartenant par hypothèse à deux groupes diffé- 
rents : on a par hypothèse 

R^'_(no2+K)R=o, 
S''_(IIp2 + K)R = o; 
on en conclut 

RS'' — S''R=o. 

On peut intégrer immédiatement et écrire 

RS'— S'R==a. 

Si on suppose R et S exprimés en fonction elliptique du para- 
mètre «, on sait que l'on pourra augmenter u de 2 to, de 2 w^ ou 
de 2 o>3 de façon que R par exemple ne change pas de valeur, et 
que S change de signe, et de même pour R' et S'; mais alors le 
produit RS' — SR' aurait changé de signe ; on a donc a = o ; 

on en conclut que -^ est une constante, ce qui est évidemment 

faux. 

Deux équations en K correspondant à une même valeur de n 
ne peuvent donc admettre de racine commune. 

Remarquons que dans l'équation en K le coefficient du terme 
de degré le plus élevé se réduit k l'unité; il est donc impossible 
que le degré de l'équation s'abaisse. 

Les fonctions de Lamé sont linéairement indépendantes. 
— Je vais démontrer maintenant qu'entre les p fonctions de 
Lamé, d'une même forme, et pour une valeur de n donnée il 

n'existe aucune relation linéaire à coefficients constants ;(/^ est 

V j i '* . '^ — I /i-f 
1 un des nombres [- i. , — 



2 
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En effet : supposons que l'on puisse avoir une relation de la 
forme 

en portant R^+i dans l'équation du second ordre, on aura à écrire 
que les p coefficients d'un polynôme seront nuls ; on aura pour 
déterminer a^ aj... a^. /> + i équations du premier degré, dont la 
première sera satisfaite identiquement. Il restera donc p équa- 
tions homogènes a q inconnues. L'élimin.ition de ces inconnues 
donnera un certain nombre d'équations en K de degré au plus 
égal à y et ces équations devront être satisfaites pour les p va- 
leurs de K ; ces équations se réduisent alors à des identités. On 
peut donc déterminer les coefficients a^ ag... a^ quel que soit 
K, et on aurait alors un polynôme comme solution de Téquation 
de Lamé, et cela quel que soit K ; c'est une contradiction. 

Les polynômes de Lamé ont leurs racines réelles, — Je 
vais maintenant démontrer qu'un polynôme R quelconque a 
toutes ses racines réelles et comprises entre a^ et c^. 

On peut partager l'intervalle de — 00 à + 00 en quatre autres 
par les nombres a^, P et c^. 

Je vais d'abord démontrer que le nombre de racines compris 

dans chaque intervalle est constant lorsque b^ varie de a^ à c^. En 

effet, il ne peut en être autrçment que si, pour une valeur de i, 

deux racines réelles devenaient imaginaires. Mais alors, pour 

cette valeur de i, ce polynôme aurait une racine double : R' serait 

nul. On a 

R''=FR = o, 

F étant (H^p + K) ; on a 

R/^/ = FR'+PR, 
Riv =FR''+ 2F'R'+ FR'' etc. ; 

toutes les dérivées seraient nulles, ce qui est impossible. 

Il pourrait encore arriver que Tune des racines devienne égale 
à a^, b^ ou c^. Supposons par exemple qu'une racine devienne 
égale à a^ : comme on a posé 

PoiNCARÉ. — Figures d'équilibre. 9 
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la valeur de u correspandante serait w^, or on a p'ui^ = o, cOj est un 
zéro de /?'co, ; donc ce serait un zéro de R' et l'on revient au cas 
précédent. 

Ceci posé soit s^ le nombre des racines imaginaires de Tensem- 
ble des polynômes correspondant à une valeur de 72 ; 

Êg, le nombre des racines réelles comprises entre a* et b^ ; 

Eg, » » réelles comprises entre b^ et c^ ; 

Ej^y » » réelles non comprises entre a^ et c*. 

Le nombre total des racines est égal à la somme des degrés des 
polynômes 



n [ n \ , 



Lorsque i* est égal à a*, toutes les racines e^ se réduisent à a^ 
et e'g se réduisent à a* ou à c^; lorsque b'^ est égal à c^, toutes les 
racines Sg se réduisent à c^ et s^ se réduisent a a^ ou à c^. 

Mais si a^ = b^ les polynômes sont les dérivées successives 

de (v^ — i)"et ils ont toutes leurs racines réelles ; il y en a un qui 

, . n 

a o racine, un qui en a une, un qui en a deux, ... un qui en a — 

abstraction faite des racines zt 1 ; on a donc 

n [ n \ 

de même on a 

Nous avons donc Tégalité 

d'où Ton conclut : 
et par suite 

Le théorème est donc démontré. 
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Pour savoir maintenant combien il y a de racines dans chacun 
des intervalles a^.., è^, è^.. 6*^,on voit que dans le cas général il j 
en a autant que dans le cas particulier, oii à^ = è^. 

Et on voit qu'il" y a : ' ~ ' 

un polynôme ayant o racine entre a^ et i% et — entre è* et c^. 

n 

» I » 1 » 

2 

n 

» — » o » 

2 

Développement d'une fonction en une somme de fonctions 
de Lamé. — Soit un ellipsoïde de la famille, déterminé par le 
paramètre p ; si je pose 



\/(P"^-l-"/(?-^-vV 
j'ai 

r/MN.M,N,rfT = o. 

MN désignant un des produits déjà désignés, M^N^ en désignant 
un autre, et l'intégrale étant étendue à toute la surface de l'ellip- 
soïde . 
Soit 

V = RMN, y, = R,M,N^. 
Puisque l'on a AV = o, la formule de Green donne 



/(v^-v.^)-=». 



l'intégrale étant étendue à une surface fermée quelconque ; sup- 
posons que cette surface soit l'ellipsoïde ; R et R^ sont alors des 

^ , , dV dY, j , , , , 

constantes. Calculons -7— et —~ i quand on se déplace sur la 
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normale à relHpsoïde, [Ji et v demeurent constants, p augmente 
de d^ ; et on a 

dn dp dn au dp dn 

_ ÔV du 
au dn ' 
maïs nous avons 

— =MN — • 
au du ' 

rintégrale devient donc 



/ 



MN. M,N, (RR', — R.R) ^ rfcr = o ; 



maïs RR'^ — R^R' ne dépend que de p et comme c'est une cons- 
tante, on peut la faire sortir du radical. On a d'ailleurs 

RR', — R,R':^o; 

il vient donc en changeant le signe 



/ 



— MNM.N, ^rfcr = o. 



Tout revient h calculer 

du du rfp 

dn dp dn ' 

nous avons vu que Ton avait, en posant p^ =pu-^ h 



d'où 



^=±V/(P^-«>^_6^)(p^_c^), 



du ' 



lorsque p tend vers l'infinî, pu tend vers oo , u tend vers o, donc 
décroît ; on doit prendre le signe moins. 
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On a d'ailleurs, puisque Télément linéaire est 
dn' = a} df + p^rfjji^ + yVv' 
et que p seul varie sur la normale : 

dn a Q 

On a donc, en posant 

- du 

^~ d^' 

— I 



ce qui est Tégalité annoncée. I/intégrale, d'ailleurs, ne dépend 
pas de p. Car nous avons 

et par suite 



IdfS : 



BCv/(pl-l^o'^-,v^) 
BC 



On a démontré que toute fonction déterminée sur Tellip-^ 
soïde peut se mettre sous la forme d'une somme de fonctions 
de Lamé; 

Le théorème précédent permet facilement de trouver les coef- 
ficients. On a 

r/<ï> Ma Na dr: ^Ç\k^mu Njk M, Nirfcr ; 
le second membre se réduit à A;t / /M|N|ûf<7 : ce qui détermine A^. 
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Les fonctions S. — Nous allons maintenant considérer une 
autre fonction de p satisfaisant aussi à l'équation de Lamé ; nous 
la définirons ainsi qu'il suit. 

Soit R une fonction de Lamé de degré n définie par la relation 

R''4-FR=o; 
la fonction S, associée à la fonction R, sera définie par la relation 

S'' + FS = o 
où F a la même valeur : on en conclut que 

S'H — R"S=o; 
en intégrant on trouve 

S'R— RS'=C^ 

et lions prenons arbitrairement cette constante égale à an -(- i ; 
on a donc 

S^R — RS^ _ an+i 

R^ — JJ2 • 

Donc 



S I an + 1 ' 



on peut choisir arbitrairement la limite inférieure d'intégration 
et on a 



S / an + I 



du. 



Pour patres grand 

R = Ap", 

A étant une constante que nous pouvons prendre égale à T unité ; 
u est alors très petit : puisque 



LES FONCTIONS S 




p^ est approximativement égal à — j 




p » » — . 
^ u 




On a donc 




R=i/-% 




s=«- T-î^ii-*- 


«»+'; 



iJ5 



S est donc sensiblement égal à m" + * = ^^^ ; le produit SMN 
satisfait aussi à l'équation de Laplace. 



CHAPITRE VII 

ATTRACTION DES ELLIPSOÏDES 
FIGURES d'Équilibre 

Problème de Diricblet sur rellipsoïde. — Abordons main- 
tenant le problème de Dirîchlet pour un ellipsoïde, dans lequel 
on se donne, sur un ellipsoïde déterminé par le paramètre p^, 
les valeurs d'une fonction V, en proposant de trouver la fonction 
V satisfaisant à Téquation de Laplace et prenant sur rellipsoïde 
la série des valeurs données. 

Soit sur rellipsoïde 

V = «1»=2a,M.N,; 

on peut poser 

A,= a,R?S?, 

RJ étant la valeur pour p = p^ de la fonction R associée à M et N, 
et S" » » S, » Kj. 

On a alors 

Ceci posé, considérons la fonction définie par l'expression 
V aS^RMN 
à l'intérieur de Tellipsoïde, et 

VaR^SMN 



à l'extérieur du même ellipsoïde. 
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Ces deux séries sont convergentes dans les espaces indiqués, 
la fonction est donc bien déterminée. 
Calculons les quantités 

ÔY 5V 
on a 

= ;y a,.R?S7M,NV 
On a de même 

= -z2a,S?R'?M,N,. 



5V 
d/2 



V est une fonction continue quand on traverse la surface, raais 
ses dérivées ne le sont pas. V est donc le potentiel d'une masse 
superficielle, dont la densité 8 est déterminée par la relation 

_dy ^__ /^g 

Au Heu de supposer la densité variable, on peut supposer une 
distribution de densité égale a i, mais d'une épaisseur Ç, telle 
que Ton ait 

équation qui détermine Ç, car on a 

4«J: = _ ^^^«..(R'.S'i - R'? s?) M,N, 

= iVa(2«-|-i)MN. 

Si l'épaisseur de la couche est déterminée par une fonction de 
la forme 



Yi ^P^N, 
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on en déduira réciproquement le potentiel sur la surface par la 
formule 



2/1+ I 



MN; 



d'où Ton déduit simplement le potentiel à l'intérieur et à Texte- 
rieur de la surface. 

Application. — Voyons maintenant les plus simples des fonc- 
tions de Lamé. 

I® pour /i = o on a 

R. = M,=N„=i, 
S=i/. 

a* pour H = I on a trois fonctions R 

On a 

y = A,R,M,N„ 
;= = /*,R,M,N„ 

en posant 



\/(f-'' — o')(c«— i») 

3" Pour n = 2, on a : 
d'abord deux polynômes du second degré 

p* — h* 
et 

p' - h'\ 



1 
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que l'on peut calculer facilement ; on sait que Ton a 

a^>h^>h':>h">c\ 
et ensuite les trois fonctions 







R, = V/(p-'-A^)(p^-0, 






R. = V^(p'_..^){p^_fl^), 






R. = v/(p^-a^)(p='-A='). 


On a 




Rj = R, Rj, 
R^ = 1^1 R2 5 


le produit R^M^N^^ 


est a un facteur constant près yz. 


» 


R,M,N, 


)) » zx. 


)) 


R,M.N, 


)) )) xy. 



Remarquons encore la signification des fonctions R^, Rj, R3, 
qui sont les demi-axes de rellipsoïde coiTespondant à p. Son 
volume est 

T = -T7- TîR^ ^2^3 ^^^ "î~ ^^1 ^4 • 

Supposons qu^à la surface de rellipsoïde, soit répandue une 
couche de densité constante, et d'épaisseur 

le potentiel à la surface sera 

V==4^R»SÎM,N,= 4««. 

les indices inférieurs indiquant le ipiméro d'ordre de la fonction 
de Lamé, l'indice supérieur indiquant qu'on y fait p = p^. A l'in- 
térieur, le potentiel sera 

le potentiel sera constant à l'intérieur; on peut donc consi- 
dérer la distribution de matière comme la distribution d'une 
couche électrique en équilibre à la surface de l'ellipsoïde sup- 
posé conducteur. 
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Le potentiel sera a Textérieur 

KRJS,M,N, = 47r«. 

Les surfaces équipoteiitielles seront donc les surfaces u = C*®, 
c'est a dire les surfaces p = constante; ce seront des ellipsoïdes 
homofocaux. 

La couche correspondante est celle comprise entre deux ellip- 
soïdes homofocaux comme le montre la valeur de Ç. 

Attraction d'un ellipsoïde homogène. — Considérons (fig. 17) 

un ellipsoïde homogène dont le potentiel soit V. Déplaçons Tel- 

lipsoïde parallèlement à Taxe Ox, d'une 

quantité c; le potentiel en un point (a;, y, z) 

P' va devenir 

V (^ — e, y, ^)=V (^,y,.-)— £ -^' 

et l'ellipsoïde prendra la position E'. La 
différence de potentiel provient de ce qu'il 
y a en plus la région ombrée à droite, et 
Fig. 17. en moins la région ombrée à gauche. On 

peut supposer que c'est le potentiel d'une 

épaisseur Ç répandue à la surface de l'ellipsoïde, ^ étant positif 

à droite et négatif à gauche. 

Calculons Ç : si PP'' est la normale commune aux deux ellip- 

soides, Ç = PP"; PP' = e, et le triangle PPT" est rectangle 

en P'', on a donc 

PP"=£cos(PF,PP'0- 
Comme 

PF' cos (PP, PP'') = projection de PP'' sur Ox, 

et, puisque PP'' est normale a l'ellipsoïde p, on a 

Ç cos P == -^^ — dp 




àx dp 

dp dn 

Oj; dp 

dp dn 



dn 
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de ces égalités on déduit 

p àx dp àx du àx 

do dn du dn au 



Nous savons aussi que 



en posant 



.4^=ARiMA, 



* dp dn 

— — A-£- 



et comme 



on en conclut 

R', = R^, 

et l'on a 

4^=AR.M,N,. 

ou * * * 

Il vient donc 

cosP = AmîM,N„ 

On peut donc maintenant déterminer le potentiel V relatif à 
la couche d'épaisseur î^, en appliquant la formule connue : on a à 
la surface 

V'=— -^ e/iRÎSJRÎM^N,. 


A l'inlérîeur de l'ellipsoïde on a 



V'=_i^eRJS;RjM,N, 
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et, puisque R,MjN, = x , 




V'=— 


sTxSÎ 
Rî ' 


on a donc 




dV 


TxSX 



^X 



r; 



A Textérieur on a 



il=i^R«R;s.MA 



K. 



dV 



A l'intérieur «y" ^st donc proportionnel k Xyk l'extérieur la 

dV 
relation est plus compliquée; -r- représente la composante de 

l'attraction parallèle à Ox, 

Théorème d'Yvory, — Considérons deux ellipsoïdes homo- 
focaux, correspondant, le premier E^ au paramètre p^,, le second E^ 
au paramètre p^; p^ > p^, (fig. i8). Soient P^ et P^ deux points cor- 
respondants, c'est à dire ayant 
respectivement pour coordon- 
nées elliptiques (p^,, u, v) (p^, 
[ji, v) . Les coordonnées rectan- 
gulaires de deux points corres- 
pondants sont proportionnelles, 
c'est-a-dire que l'on a 

X, =px,. 




Fig. 18, 



rz. 



/?, y, /• étant seulement fonctions de p^ et de p^. 

Le théorème d'Yvory consiste à comparer l'attraction de l'ellip- 
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soïde E^ sur P^ à Tattraction de Tellipsoïde E,, sur Pj. La pre- 
mière a pour composantes suivant Ox 

47rARlSÎR;M,N,, 
la seconde a pour composantes suivant Ox 

47rARîR;SîM,Nj, 
le quotient est donc 

on trouverait de même que les rapports des attractions sui- 
vant Oy et Oz sont respectivement 

^' et '^^ 

Ellipsoïde de Mac Laurin, — Considérons un ellipsoïde 
homogène, que nous faisons tourner autour de l'axe Ox^ d'une 
rotation infiniment petite o) ; le point de coordonnées x^ y, ^, 
aura pour nouvelles coordonnées 

X, t/-\-ioz, z — iùy] 

son potentiel sera 

/ dV ôV \ 

V(^, y + -w, z — yiù)=\{x, y, ^)+^\^-j- y-^l' 

V — V est le potentiel dû à la couche superficielle comprise 
entre les deux ellipsoïdes. 

Si [x^ y, z) représente un point intérieur, Faction en ce point 
a pour composantes 



T^s; TijSi izs\ 

Rî ' R'î ' R? 





' 



elle dérive d'un potentiel 



i44 ATTRACTION DES ELLIPSOÏDES 

en appelant V^^ le potentiel au centre ; et cette formule est encore 
vraie à la surface : donc 

— et — 

sont proportionnelles à y et à ^, et on a 

V — V = Aa)yc 

= Aa)R,MA. 

En appliquant la formule démontrée précédemment (p. 137), 
l'épaisseur de la couche en un point est 

5AW M^N^ 



4« Sî 



Pour que Tellipsoïde puisse être une figure d'équilibre, il faut 
qu'à la surface on ait 

Nous allons traiter un problème plus général ; de savoir si un 
ellipsoïde à trois axes peut être une figure d'équilibre, sachant 
qu'il est soumis, outre les attractions mutuelles de ses parties, à 
une force dérivant d'un potentiel 

a, ^, Y étant des constantes. 

Pour que l'équilibre ait lieu, il faut que 

soit constant à la surface de Tellipsoïde dont Téquation est : 

^,f_^^ 

Kj "*" i\i '^ \\i ' ~" °' 

Il faut donc que l'on ait 
(1) aRJ-TS.R, = ?Rt-TSA=TR|-TS3R3. 
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Dans le cas particulier qui nous occupe, c'est-à-dire quand 
Tellipsoïde est en rotation, on a, a = o, p == y = a>^, et nos 
formules deviennent 

TS,R, = 0)^^ — TS,R, = o)2R'^ — TS3R3 ; 

on a donc 

(AIJ T"- P R| -^ Ri ' 



et d'ailleurs 



T--= -y-R^RaRg; 



T est une donnée du problème où les inconnues sont p^ — a*, 

On a une première condition, qui est que w* doit être positif; 
donc il faut que 

RjS, — R^Sj>o, 
R3S3-R,S,>o. 

Ces inégalités sont toujours vérifiées. Démontrons par exemple 
que la première Test, c'est-à-dire que Ton a 

R^S, — R,S,>o. 

Remarquons que R^, R^, S^ et Sg sont des quantités positives, 
puisque p^ > a^ ; il faut donc que l'on ait 

s. > R. ' 



c'est-à-dire 






PoiNCARÉ. — Figures d'équilibre. 
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f du 

s. 



02 » 

tu ^ , > 



et on est ainsi ramené à un théorème d'arithmétique connu. Le 
premier rapport est compris entre la plus grande et la plus 
petite valeur du rapport des quantités placées sous le signe/: la 
plus grande de ces valeurs est i correspondant a m = o, p = 00, 

la plus petite est 1 rj ; le théorème est donc démontré. On 

en donnera plus loin une démonstration plus générale. 

On voit de suite que Taxe de rotation ne peut pas être l'axe 
Oy, car il faudrait que R^S^ soit la plus grande dés trois quantités 
analogues. 

La relation (II) a une solution évidente, c'est R2 = R3, qui con- 
duit à un ellipsoïde de révolution autour de Ox : c'est l'ellip- 
soïde de Mac Laurin. 

Ellipsoïde de Jacobi. — Les relations (II) peuvent aussi être 
satisfaites pour un autre ellipsoïde découvert par Jacobi, mais 
nous devons, avant de le rechercher faire la remarque suivante : 

Si nous nous donnons un ellipsoïde homogène quelconque, 
nous pouvons toujours choisir les constantes a, p, y, de manière 
a satisfaire aux relations (I) c'est-à-dire de façon que cet ellip- 
soïde soit en équilibre sous Faction de la gravitation et d'une 
force supplémentaire dont les composantes sont ox, j3z/, ys ; nous 
allons supposer a = o. 

La condition d'équilibre d'un ellipsoïde homogène E soumis 
h la gravitation et à une force perturbatrice dérivant d'un poten- 



tiel — (Py^ -j-r 'p^) est que, à la surface 

soit une constante. 

Supposons que Tellipsoïde se déforme mais peu, de façon que 
la nouvelle surface S ne diffère pas beaucoup de l'ellipsoïde, et 
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cherchons à quelle condition S sera une figure d'équilibre, en 
continuant à supposer que le fluide reste soumis aux mêmes 
forces, c'est-à-dire à la gravitation et a la force ou:, j3y, y^. 

La condition d'équilibre est que, en désignant par y -\- \f \q 
potentiel dû à la figure S ; 

soit une constante à la surface. 

Soit Ç la longueur de la normale comptée de E à E ; on pourra 
développer !J en une somme de fonctions sphériques : 

et Ton aura à la surface 

^ 2/1+ I 

D'autre part, on aura, à des quantités près de l'ordre de Ç^, 

et, en appelant g l'intensité de la pesanteur à la surface, il vient 

dU 

d'où 

La condition d'équilibre est que U soit une constante et comme 

la condition d'équilibre est que 
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soit une constante a la surface c'est-à-dire que 

Je dis que gl est une constante sur l'ellipsoïde ; en effet, 
déplaçons Tellipsoïde parallèlement h Taxe Ox, l'équilibre ne 
sera pas troublé, car par suite de la condition a = o, le travail 
de la force (ox, ^y, ^z) est nul. Nous avons vu que Ton avait 
pour ÎJ la valeur 

î:=Km,N„ 

K étant une constante ; Téquation précédente se réduit en ce 



cas a 



^ RJSÎM^N, _^;KM,N, = a^ ; 



ceci ne sera constant que si 

4k 



gl=^^^\^\. 



3 
L'équation (i) devient alors 

et Téquation d'équilibre devient, en supprimant l'indice zéro, 
R,S, R.S. 



2/1 -f- I 



= o. 



Si la masse est en équilibre sous l'action de l'attraction et de 
la force centrifuge, elle le sera encore si on lui donne un mouve- 
ment de rotation. Cette rotation est une des transformations que 
l'on peut faire subir à un ellipsoïde de Jacobi, sans troubler 
l'équilibre. 

Nous avons vu que dans ce cas, on avait h un facteur constant 
près (p. i44) 

Ç=M,N,; 
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on en conclut que rellipsoïde de Jacobi doît satisfaire à la relation 

5 — 3 • 

Réciproquement, supposons qu'un ellipsoïde soit en équilibre 
sous Taction de l'attraction et d'une force dérivant du potentiel 

je dis que si on peut le faire tourner d'un angle w sans troubler 
l'équilibre, c'est que P=Y. 

En effet, à la surface on doit avoir, en appelant V le potentiel 
newtonien 

V + ^(Py^ + Y=^) = V, = C''; 

au sommet du petit axe V = Vj; si on fait tourner l'ellipsoïde 
d*un angle o) autour de Ox, le sommet du petit axe restera sur la 
surface de l'ellipsoïde, la constante V^ n'aura pas changé, V n'aura 
pas changé non plus. 
Or 

il faudra donc que sur la surface 
et en tenant compte des relations 

0^ == tùZ 

on devra avoir sur la surface 

(P_y)o>yr. = o; 
donc 
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Discussion de l'équation des ellipsoïdes de Jacohi, — Nous 
avons donc a discuter l'équation des ellipsoïdes de Jacobi 



Cette équation ne convient pas aux ellipsoïdes de révolution^ 
cela provient de la façon dont on a établi cette équation : dans le 
cas d'un ellipsoïde de révolution, on aurait en effet Ç= o. 

Au lieu de discuter Téquation (i),nous allons discuter l'équa- 
tion plus générale 

r= ; ; = 

2772+1 271+1 

OÙ Ton a p^ > a^ ; 772 étant l'ordre de la fonction R^, n étant 

l'ordre de la fonction R,. 

F 
F et -rr^ sont de même signe, de plus R;^ ne peut pas s'annuler 

puisque les racines des polynômes R sont comprises entre «" 
et c^ ; donc l'équation 

p _ F _ S, S, R? 

'- n~ ['^ni+ I) R, (2« + i) R, RJ. - 

a les mêmes racines que F=o. 

F 
Les racines de l'équation -pj^ = o sont séparées par celles 

de la dérivée : 

dl\ 1 d S;, 1 R? d S, 



du 2772-1- 1 du R;^ 271+1 RJ du R^ 

S; d I R/ 



272 + 



\ S, d /Rf\ 

4- I K, du U? / ' 



On a par définition 



d Sj. 2OT+I 
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et par suite 

dV^ ^ 2 S. R; R',R,,— R,R;. 

du ~ 2rt+ 1 Ri ' Ri R| 

^ 2S, r:.r,.-r,r;. . 

mais S^ ne peut s'annuler : donc les racines de \ ^ sont les 

du , 

racines de RfR/, — RiRI. Nous ne nous occupons évidemment 

que des racines supérieures \\ a^. 

Cherchons la dérivée de cette fonction par rapport à m, c'est : 

R'/R,-R,R','; 
et, d'après la définition de R,- et R^, 

R;' = (Hp^+K)R„ 
on a 

R,R,[[«(«+i)-,«(m + i)]p^ + K,-K*), 

Les deux premiers facteurs n'ont pas de racines en p^ supé- 
rieures à a^; le troisième ne peut évidemment en avoir qu'une 

d¥ 
seule : donc — — ne peut en avoir plus de deux, F^ ou F ne peu- 
vent en avoir plus de trois. 

R 
Supposons que -—- soit constamment croissant, sa dérivée ne 
Ra 

s'annulera pas; donc Fi sera constamment décroissant et ne 

, , . » . . Rf 

pourra avoir qu une seule racine. Au contraire, si rr— est toujours 

décroissant, Fi sera toujours croissant. D'ailleurs pour m = o, 

R. 

p = 00 , -^ est a peu près égal h — ; donc F part de zéro et est 

toujours croissant ou toujours décroissant, donc n'a pas de 
racines, sauf la racine p==oo qui n'est pas admissible. 

Dans le cas le plus général, F peut avoir, naus l'avons vu, 
trois racines : mais parmi celles-là se trouva la racine • p =iûOf. , 
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F a donc zéro, une ou deux racines supérieures à a^ \ ^^IT' ^^* 

toujours croissant ou toujours décroissant, elle n'en a aucune. 
Revenons alors à Téquation 

o 2/1 4-1 

Si Ri contient en facteur y/p^ — a*, ("w"/ ®^**^ ""^ polynôme 

dont toutes les racines sont réelles et comprises entre a^ et 
c^; donc la dérivée ne peut s'annuler pour une valeur de p^ supé- 
rieure à a^. L'équation F n'a pas de racine supérieure à a^. Sup- 
posons maintenant que R, ne soit pas divisible par y/p''' — a^ , et 
que n soit plus grand que i, je dis que l'équation a au moins une 
racine. En effet, substituons a^, le premier terme s'annule, le 
second est négatif : on a donc une somme négative ; remplaçons 
p par 00 ; le premier terme de l'équation est a peu près égal à 



JL(A ^V 

P \ 3 2/1+1 /' 



le résultat est positif; il y a donc un nombre impair de racines 
entre a^ et + oo , il ne peut y en avoir plus de deux, il y en a donc 
une et une seule. 

Ce raisonnement n'est pas valable si /i == i , c'est-à-dire si 











R. 


= R,= 


=v/p'^ 


-b\ 


ou bien 








R, 


= R,= 








-\/?' 


-c\ 


dans ces 


deux 


cas, 


la 


quantité 
















Rf 
Rî ~ 




-a' 



et les quantités analogues sont toujours croissantes lorsque p? 
varie de oo à a^. Donc l'équation correspondante n'a pas de 
raciaes : c'est. un théorème déjà énoncé (p. i45). 



f 
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Appliquons ce principe à Tellipsoïde de Jacobi. Il faut consi- 
dérer l'équation 



3 



o; 



R^ n'est pas divisible par R^ : donc il y a une racine, et une seule. 
Donc parmi les ellipsoïdes homofocaux, il y a un ellipsoïde 
et un seul, qui soit une figure d'équilibre ; c'est l'ellipsoïde de 
Jacobi. 

Figures dérivées de l'ellipsoïde de Mac Laurin. — Suppo- 
sons que b^ = c^ et cherchons s'il y a des figures peu différentes 
de l'ellipsoïde de Mac Laurin. Pour qu'il existe une figure d'équi- 
libre S peu différente de cet ellipsoïde, il faut que l'on ait la 
condition 

^'^ -3~— i^+T- 

Nous avons vu qu'en posant 

on avait 

N, = cosyrf, 

h étant une constante. On a posé 

. F(i)=(i— $^)-^"D'»+^(i — iO». 

La condition nécessaire et suffisante pour que l'équation (i) 
ait une racine est que R, ne soit pas divisible par yp'^ — a^; donc 
M, ne doit pas être divisible par \/a^ — [jl% donc n-\-jj doit être 
pair, n etp étant déterminés, une figure d'équilibre voisine d'un 
ellipsoïde de révolution est déterminée par la relation : 

Ç = £ f M, cos po 

ou, d'une façon plus générale, 

^ = ZM* (e cos p'f + e' sin p'f) . 
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On peut supposer £' = o : cela revient k faire tourner les axes 
de coordonnées d'un angle convenable autour de Ox, La figure 
que Ton obtient ainsi est une figure admettant Taxe Ox comme 
axe de symétrie d'ordre/?. F (i) ne change pas si l'on change \ en 
— Ç, donc le plan des yz est un plan de symétrie, Sip = o, Ç ne 
dépend pas de '^, la figure est de révolution. 

On ne doit pas conserver parmi les équations obtenues celle 
où n = Oy car le volume ne serait pas conservé ; de même on doit 
exclure /i = i, car alors R, serait divisible par y p^ — d^ : donc 
la valeur minimum de n est 2 ; à cette valeur de n correspondent 
deux valeurs de/? : p = 0^ et p==2. 

Quelles sont celles de ces figures qui sont ellipsoïdales : con- 
sidérons un ellipsoïde E^ et un autre E, qui en soit peu différent 
mais non homofocal. Le potentiel de chacun d'eux à son inté- 
rieur est une fonction du second degré des coordonnées carté- 
siennes ; la différence des potentiels est donc aussi un polynôme 
du second degré : c'est d'ailleurs le potentiel dû à la couche com- 
prise entre les deux ellipsoïdes. On peut par suite développer ce 
potentiel sous la forme 

VeRMN, 

R, M, N étant des fonctions sphériques d'ordre o, i ou 2. 
A la surface on a 

V=2 eR,MN, 
et à l'extérieur on a par conséquent 

V= V £ -|^ SMN, 
l'épaisseur de la couche produisant ce potentiel est donc 



.=1- 
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Il n'y a que des termes d'ordre o, i ou 2, mais nous avons dit 
que M ne pouvait être d'ordre o ou 1 ; il reste donc /i = 2 et 
^; =0 ou /^ = 2. Si/? = G, la surface considérée est de révolution, 
c'est un ellipsoïde peu différent du premier; sip =2, on a un ellip- 
soïde de Jacobi. 

Variations de w^ avec r aplatissement. — Si on fait varier 
p^ de a^ à + 00 , l'aplatissement des ellipsoïdes considérés varie 
de o h I. 

Lorsque l'aplatissement est nul o>^= o ; l'aplatissement crois- 
sant, w'^ va aussi en croissant; mais il peut avoir des maximums et 
des minimums; lorsque l'aplatissement augmente, le moment 
d'inertie principal augmente aussi. 

On a vu que 

TU, = H = W + -^ J 

ne pouvait pas dépasser une certaine limite. On a de plus 

2 

2 
dans une position d'équilibre on" a 

rfW+— rfJ=o; 
2 

il en résulte que pour une position d'équilibre on a aussi 

rfH==ra> Jé/io* 

Lorsque l'aplatissement augmente, J augmente indéfiniment ; 
donc (O doit diminuer et tendre vers zéro ; tù^ part de zéro, revient 
a zéro, donc a un maximum dans l'intervalle. 

Remarque. — Faisons varier p^ de 00 a a^. On obtient des 
ellipsoïdes de plus en plus aplatis : on en rencontrera pour les- 
quels on aura 

R^S, ^ R,S/ 

3 ~ 2/ï+l ' 
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je dis que la première figure d'équilibre que l'on rencontrera 

T> 

est celle d'un ellipsoïde de Jacobi. En effet, si on suppose que -=^ 
est toujours croissant dans l'intervalle considéré, on voit que 



RyS^ R,St 



5 2/1+1 

Pour p^ = x> on a 



est toujours positif. 



R,S, ^_RA_ R/S. 



>— T-^> 



3 5 ^ 2/1 + 1 ' 

si p décroît, on finira par avoir 

R^Sj Ri^i R|St 

~3~""~5~"^ 2/2+1 ' 

p continuant à décroître, on aura alors 

RjS, R» S, 



3 2/2+1 ' 

mais on voit que l'on rencontre d'abord un ellipsoïde de Jacobi. 

Figures dérivées de l'ellipsoïde de Jacobi. — Il reste à étu- 
dier les figures d'équilibre qui différent peu de l'ellipsoïde de 
Jacobi. Pour cet ellipsoïde on a 

RiSj R4S4 . 

~3~~ 5 ' 

pour une figure peu différente, on aura 

RjS, RjSt 

3 ~ 2/2 + I ' 

la hauteur de la surface au-dessus de l'ellipsoïde étant 

On a donc à considérer les deux équations 

R,S, __ R,S, _ RjSj 
3 5 2/1 + 1 ' 
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On a trois inconnues p- — o}, p^ — J^, p^ — c^ ; et on a pour le 
volume qui est donné 

Pour que Téquation 

RjS^ Ri S, 



6 2/2 + 1 

ait une racine, il ne faut pas que R,- soit divisible parRj=y p^ — d^; 
pour que l'équation 

R.s, ^ R,S, 

5 2/1 + 1 

R 

ait une racine il faudra que -^-^ ne soit pas constamment croissant, 

R, a donc comme seules formes possibles les formes suivantes, où 
P désigne un polynôme en p^. 

p, 

Pxy/ p^ — ^% 
Pxy/ p^— c-% 

pxv/(p'— ^')(p'— ^')- 

Voyons si ces quatre formes sont possibles et rappelons que 

l'on a B.,=s/[?' — l>'){?' — <^'). 

1" Type. Rj est un polynôme ; si a est la plus grande racine 
de R, on pourra poser 

R, = (p^-a) n,. 

« est compris entre a^ et c* et on a 

-^=i/ p' — » t/ p'— « n.. 

R, V p-^-0' -VT^^::^ " 

Ilf croit constamment, 



^/', 
^ 



£; croît constamment, 

p'' — c^ 



^ croit si a est inférieur h t^, décroit dans le cas con- 



1 
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traire ; il ne faut pas que i/ -^ ^ soit décroissant. : donc aL<b' 

et toutes les racines de R^ sont par suite comprises entre b^ et c^. 
Parmi les fonctions de Lamé d'ordre /i, il y en a une et une 
seule ayant toutes ses racines comprises entre b^ et c^ ; c'est celle- 
là que Ton doit prendre pour R^. 

2^ Type. Prenons R, = y/p^ — // P et appelons a la plus 
grande racine du polynôme P ; on a 



R, ^^ V p^— 6'^ 



Les trois facteurs sont constamment croissants, donc toutes les 
fonctions de ce type sont à rejeter. 
3® Tyj?e. Posons 

Ri=v7^=^(p''-«)n. 

n étant un polynôme dont les racines sont inférieures à a; nous 
avons 



Ir-N/p-'^N/IËf-- 



Pour que l'équation puisse avoir une solution, il est néces- 
saire que 0L<^P. Donc toutes les racines de Psont comprises entre 
P et c^ 

4® Type, On voit de même que les fonctions du quatrième type 

- R 

sont telles que -~ soit un polynôme sans racines supérieures à 

a^, donc il n'y a pas à considérer de fonctions du 4*^ typ^- 

Donc pour qu'une figure peu différente d'un ellipsoïde de Jacobi 
soit une figure d'équilibre, il est nécessaire que R,- soit une des 
fonctions que nous avons appelée Ro.„. Je dis que cette condition 
nécessaire est aussi suffisante : en effet, supposons que b^ varie 
de c^ à a^; à chaque valeur de P correspond une valeur de p^ telle 
que 

R,S. _ R,S, . 

3 ~~5~' 
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considérons alors la fonction, 

5 . 2/1 + 1 ^ 

pour P = c^, ~- est constamment croissant, donc F est positif; 
si 6^ = a^, la valeur correspondante de p^ tend vers a^j alors 

R^ =^0 puisque R, ne contient ni y/p^ — d^ ni y/p^ — c*, et que toutes 
ses racines sont comprises entre b^ et c^ : donc 

F<o 

et il y a certainement un système de racines pour les équations 
considérées. 

Nous voyons donc qu' il y a des figures d'équilibre voisines des 
ellipsoïdes de Jacobi. 

Remarque. — Quelle est la première figure d'équilibre que 
Ton rencontre ? c'est-à-dire, pour quelle valeur de n la valeur 
correspondante de p est-elle maximum? 

Considérons deux fonctions Ro,^ et R„„, 

/?2>/i; 

je dis que 

Ro.m Rat 



Ro,n R* 



est constamment croissant; nous allons montrer que la dérivée, 
prise par rapport à w, est constamment positive ; cette dérivée 
est à un facteur positif près 

R'/cRi — R/cRo 

si Rj n'est pas divisible par y/p"^ — a^, R/ l'est; et réciproquement. 
Mais ici R, et R/^ ne sont divisibles ni paryp^ — d^ ni par y p^ — 6^; 
il en résulte que 

R',R,-R,R'; 

s'annule pour p^ == a^ 



n 
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Si on exprime Rj et R;^ en fonction de l'argument », on a 

R, = F(p«), 

RJ=/«F'[/,(«)], 

et on voit que la dérivée n'a pas de racine supérieure à a^. En 
effet la dérivée de cette fonction (prise par rapport à u) est 

(Ri'R,-R,R;') = (?*-?.)R.R* 

comme on Ta déjà vu; R, et R;^ ne peuvent s'annuler, «p, — «p^^ ne 

peut s'annuler deux fois, puisque c'est un polynôme du premier 

degré en p^. La dérivée ne peut donc s'annuler plus de deux fois 

pour a^ ^ p^, c'est-à-dire pour o < u ^ e^, 

La dérivée qui s'annule pour u = o, 11 = e^^ ne peut donc s'an- 
r) 
nuler dans l'intervalle; -~ est donc constamment croissant lorsque 
*^» 

p'^ varie de a^ h -|- oo . Il en résulte que l'on a aussi 



2Ae- 



lorsque p^ varie de 00 à d^^ on rencontre donc d'abord une 
valeur satisfaisant à l'équation 

5 2/1 + I ' 

puis ensuite la valeur de p^ satisfaisant à l'équation 

RA _ R*s, 



5 2/1 + 1 

La première des valeurs de n que l'on trouve est donc /i = i, 
puis /i = 2 : on a déjà examiné ces valeurs ; il faut maintenant 
examiner la surface voisine de l'ellipsoïde de Jacobi et corres- 
pondant à Rç3 ; nous allons étudier d'une façon générale la 
surface correspondant à R(j„. 

L'épaisseur de la couche sera 

t; = £ZMN,; 
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cette épaisseur sera nulle aux points où l'on a 

M,N,==o, 

ou bien 

RiM,Ni = o. 

Ceci est décomposable en facteurs quadratiques de la forme 



^2 TT ~t- 



a — a 



i; 



les a sont d'ailleurs compris entre P et c^ ; les surfaces corres- 
pondantes sont donc des hyperboloïdes a deux nappes coupant 
l'ellipsoïde suivant des 
lignes de courbures 
d'un seul système. 

D'un côté d^une de / / 
ces lignes ^ sera positif; / / 
et de l'autre côté Ç sera \ \ 
négatif. On peut, avec \^ \ 
ces données construire 
la figure dans l'hypo- 
thèse de 71 == 3. On re- 
présente (fig. 19) en 
ponctué le contour apparent de l'ellipsoïde, en traits pleins l'in- 
tersection de la surface avec l'ellipsoïde ainsi que le contour 
apparent extérieur a l'ellipsoïde, en pointillé le contour apparent 
de la surface intérieure à l'ellipsoïde. L'une des lignes de cour- 
bure est d'ailleurs l'ellipse principale, comme on le voit en re- 
marquant que 




Fig. 19. 



R« 



On construira avec les mêmes conventions la figure d'équilibre 
correspondant à R^^^ (fig. 20). 

Je remarque qu'on ne peut prendre pour les fonctions R,- de 
fonction du second ordre, car il faudrait prendre pour la fonction 
Rq2 un polynôme de la forme p^ — A^, ce qui donnerait un autre 
ellipsoïde et nous avons vu qu'il ne pouvait y avoir qu'un ellip- 
soïde de Jacobi pour une vitesse de rotation donnée. 

PoiNCARÉ. — Figures d'équilibre. 11 
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En résumé, si on faît varier b^ de c^ à a^, on obtient des ellip- 
soïdes à chacun desquels correspond une valeur w et une seule ; 
(0 ne peut avoir de maximum ou de minimum, car s'il y en avait 

un (Op à une valeur voi- 
sine de (Oj correspon- 
draient deux figures 
ellipsoïdales ce qui n'est 
pas. Donc w varie tou- 
jours dans le même 
sens. 

Supposons que b'^ et 
c^ tendent vers zéro, la 
Fig^. ao. forme de l'ellipsoïde 

tend alors vers une ai- 
guille allongée; (o a d'ailleurs la même valeur pour deux ellip- 
soïdes semblables. Il suffit donc de voir ce qui se passe dans le 
cas où 6 et c restant fixes, a augmenterait indéfiniment. On 
aurait alors un cylindre elliptique, la force résultante de l'at- 
traction du cylindre sur lui-même et de la force centrifuge 
doit être normale au cylindre, c'est-à-dire parallèle au plan yOz, 
Considérons un point d'une génératrice passant par l'axe de 
révolution : par raison de symétrie, l'attraction doit être parallèle 
au plan yOz^ il en résulte que la force centrifuge qui est per- 
pendiculaire à Os ne peut être que nulle, et par suite w tend 
vers zéro lorsque a augmente indéfiniment. 

On peut d'ailleurs remarquer que le cylindre ne peut être 
que de révolution pour que l'équilibre soit possible. 
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Représentation graphique des résultats précédents. — 
Prenons deux axes de coordonnées (fig. 21) et portons en abs- 

R2 R2 

cisses les valeurs de -p|- et en ordonnées -p|-; ces quantités sont 

les rapports de l'axe moyen et du grand axe au petit. 

On a un premier point A correspondant au cas où ces rapports 
sont égaux à 1, la figure correspondante est une sphère, La partie 



i/ 


i II.' y/H 

i Vk- /^ 




/f 





JC 



Fig^. ai. 
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de la bissectrîce de l'angle xOy partant du point A correspond 
aux ellipsoïdes de révolution. Les ellipsoïdes de Jacobi sont 
représentés par une courbe CBC symétrique par rapport à cette 
bissectrice et la coupant en 
un point B. 

Pour les ellipsoïdes de 
Mac Laurin, (o^ va en crois- 
sant ({u&nd on part en crois- 
sant du point A, jusqu'à ce 
que l'on arrive ea un point D, 
puis il décroît quand on s'é- 
loigne du point D jusqu'à 
l'infini où il arrive h la 
valeur zéro. Sur la courbe 
CBC deux points symétri- 
ques par rapport à AH re- 
présentent le même ellip- 
soïde ayant tourné de 90°; o>^ décroit quand on part de B 
pour aller vers C ou vers C; les points C et C s'éloignant indé- 
finiment correspondent au cas du cylindre de révolution, la 
courbe CBC est asymptote aux droites x = i , y = i , 

Nous marquons sur AH les points E, F, correspondant aux 
figures peu différentes de l'ellipsoïde de révolution, de même sur 
CBC, nous remarquons les points GKLM, G'K'L'M', correspon- 
dant aux valeurs Rj,„ donnant les figures peu différentes d'un 
ellipsoïde de Jacobi, que nous avons rencontrées précédemment. 

11 faut remarquer que les paramètres que nous avons pris ne 
correspondent pas aux vraies données du problème. Les para- 
mètres sont le volume T, la rotation *o et le moment principal 
d'inertie J. Si on pose [ji= w^J, jx et T sont des données du 
problème, puisque w^J doit être constant; [x est du cinquième 
ordre par rapport aux longueurs, T est du troisième ; si on 
pose 

T 3 

M sera une donnée du problème indépendante de l'unité de lon- 
gueur choisie. 
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Liouville a cherché à quelles valeurs de M correspondaient les 
divers points des lignes AH et CBC. 11 a montré que lorsqu'on 
se déplace sur AH, M croît lorsqu'on va de A à D et décroît 
si on va de D vers H. Sur la ligne CBC, M décroît lorsqu'on 
s'éloigne du point B soit vers C, soit vers C. 

Donc si (0 est très grand, il n'y a pas de figure d'équilibre 
possible; si o> est compris entre (Og et (o^, il y a deux figures pos- 
sibles qui sont des ellipsoïdes de révolution. 

Si (0 est inférieur à w^, il y a quatre figures d'équilibre pos- 
sible, dont deux sont des ellipsoïdes de Mac Laurin, et deux des 
ellipsoïdes de Jacobi ; ces deux ellipsoïdes de Jacobi sont d'ailleurs 
égaux comme on l'a déjà remarqué. 

Courbes d'équilibre. — Considérons un système dépen- 
dant de n variables ^\^2 ••• -^n» soumis à un système de forces 
dépendant d'un potentiel F [x^x^.,, x^y X) ; la condition néces- 
saire et suffisante de l'équilibre est que F soit maximum par 
rapport aux Xy \ étant donné ; on a ici 

F=W-^J=W- '^^ 



2 2J 

La condition d'équilibre est que 

dF__^_ _ 5F _ 

pour que l'équilibre soit stable, il faut de plus que la forme 

soit une forme définie négative ; on peut la mettre sous forme 
d'une somme de n carrés : 



^ =2j °^fc^*'- 



Ces coefficients a^. sont appelés coefficients de stabilité ; si l'un 
de ces coefficients s'annule, le discriminant de la forme $ s'an- 
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nule, puisque cette forme est alors réductible à une somme de 
moins de n carrés ; ce déterminant est, en posant 



d^F 



A = 



^1,1 ^1,2 
^2,1 ^2,2 



'''f.fc, 






Si deux coefficients dj^ sont nuls, tous les mineurs du premier 
ordre de A sont nuls ; nous supposerons qu'il n'en est pas ainsi. 
Si A est différent de zéro, on peut résoudre les équations 



5F 

^x. 



5F 



5F 



et on aura 



Différentions la relation 
VI d^F 



5F 

dxi 5^F 



o par rapport à X, on aura 



5^-^50^, 



d\ 



5x,5X 



ce qui peut s'écrire 



dxi 



5T 
5a;, 5).* 



dXi 
Si A est différent de zéro, on peut trouver les valeurs de ,,' , 

puisque l'on a n équations linéaires a n inconnues et les x seront 
des fonctions uniformes de X, dans le voisinage des valeurs con- 
sidérées. 

Supposons A qui, sans que tous les mineurs du premier ordre 
le soient, supposons que l'on ait par exemple 



5A 



^o. 
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Nous pouvons alors résoudre les n — i dernières équations 
linéaires par rapport à 

dx^ djc^ dXfi 

ce qui montre que j:,, jTg, ... x^ seront fonctions uniformes de x^ 
et de X. 

Substituons dans la fonction F les valeurs de x^^x^^ ... x^^ en 
fonction de x^ et \ tirées des équations 

5F _ ^ 

on aura une fonction i^{x^^V)y et on a 

^_ilL_i_^_^_j_ . ^F dx, 
dx^ l 

mais par hypothèse 



dx^ ôXj ^x^ dx^ " dxn dx^ 



ôF _ 5F _ 

5^2 '* ' ^-^n 

quand on y remplace x^x^.,.Xn en fonction de x^ et de X ; on a 
donc 

df^ _ 5F 

dx^ 5^^ 

La condition d'équilibre est alors que l'on ait 

4 _ 

dx^ 
C'est une relation entre x^ et \ qui définit une courbe, x^ sera 

»2 f 

une fonction uniforme de \ si ^ est différent de zéro ; la 

CiX^ 

J. . , , . . . d'ij 

condition A:=o équivaut donc a , ', =o. 

dx^ 

Supposons que cette condition soit remplie sans que l'on ait 
en même temps --7 — ^^-==0, la courbe a alors une tangente ver- 
ticale (fig. 22). 



ÉCHANGE DES STABILITÉS 



167 



Quand X traverse la valeur de Tabscisse de cette tangente, il y a 
pour x^ deux valeurs réelles quî viennent se confondre, puis 
deviennent imaginaires ; de même 
pour ^2, ... ;r„. 



Si 



: o, le point corres- 




pondant de la courbe est un point 
double, à tangentes distinctes ou 
confondues ; Tune des tangentes 
peut être verticale ou non, on a 
deux dispositions principales : 
celle qui correspond à deux tan- 
gentes distinctes et non verticales (fig. aS) et celle qui 
correspond à une tangente verticale et une non verticale 

(fig. 24). 



Fig. aa. 





Fig. a3. 



Fig. a4. 



On a deux cas à distinguer : dans le premier cas, deux valeurs 
réelles de x se confondent et redeviennent réelles quand on tra- 
verse la valeur correspondant au point double. 

Dans le second cas, deux groupes de valeurs réelles se confon- 
dent et deviennent ensuite imaginaires lorque X traverse la valeur 
critique dans un sens convenable, une autre valeur de x est 
réelle avant et après la valeur critique. Les autres dispositions de 
figures peuvent se ramener à celles-là, du moins au point de vue 
qui nous occupe. 

Échange des stabilités. — H y a équilibre stable, avons- 
nous dit, lorsque F est maximum, x^x^ ... Xf^ étant variables. Il 
faut a fortiori que F soit maximum si x^ étant donné on fait 
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varier x^^ jTg, ...,^« d'une façon arbitraire, il faudra donc que 
ôF dF_ dF _ 

On tirera ^2^ ^37 • -j-^n de ces équations dont le déterminant 
fonctionnel n'est pas nul par hypothèse, on portera ces valeurs 
dans F [oc^x^ ••• •^n)i ®t ^^ ^^'^^ ^^® fonction qui devra être maxi- 
mum pour qu'il y ait équilibre. Donc i/ {x^y X) doit être maximum ; 

il y aura donc équilibre si -j^ = o, mais il n'y aura équilibre 

stable que si 1 on a en même temps . ^ < o. 

Un arc de la courbe -7-^ = partage la région du plan où 

il se trouve en deux parties, l'une où /^ > o, l'autre où , ' <o ; 

SI c est un arc sans singularité, , \ ne sera pas nul et sera 

négatif si la région où —7-^— est positif est située au-dessous de 

l'arc de courbe ; et au contraire sera positif, si /^ est positif au- 
dessus de l'arc de courbe considéré. Le premier cas correspond à 
une position d'équilibre instable. 

En couvrant de hachures la région où est positif, on voit 

Q/X . 

que la figure 25 correspond au cas d'un équilibre stable et la 
figure 26 au cas d'un équilibre instable. 





Fig. 25. 



Fig. 26. 



Dans le cas où l'arc BAC a une tangente verticale en A, deux 
cas peuvent se présenter, suivant que c'est la partie convexe 
ou la partie concave qui est hachurée. 
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Dans la figure 27 : 

Tare BA correspond à des positions d'équilibre instable, 
l'arc AC » » stable; 

dans la figure 28 : 

l'arc BA correspond a des positions d'équilibre stable 
l'arc AC » » instable. 



y 




Fig. «7- 




Fig. a8. 



Ce qu'il y a d'essentiel, c'est qu'en suivant l'arc et en passant 
par le point A, les deux systèmes de figures d'équilibre échangent 
leur stabilité. 

Dans le cas d'un point double sans tangente verticale, deux 
cas peuvent aussi se présenter; dans les figures 29 et 3o, les 





Fig. 29. 



Fig. 3o. 



arcs BA et AC correspondent à des figures d'équilibre instable, 
les arcs B'A et AC correspondent à des figures d'équilibre stable. 
Remarquons encore qu'en suivant l'arc BAC ou l'arc B'AC les 
figures d'équilibre, de stables deviennent instables ou récipro- 
quement. 

Enfin dans le cas d'un point double avec une tangente verti- 
cale, on arriverait à la même conclusion, les figures d'équilibre 
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que Ton obtient en suivant un arc quelconque échangeant leur 

stabilité quand on passe par un point où ^-V= ^ (^8^ ^^ ®* ^a). 




jr 


liiiiiiiiifrK. ' 


, 






# 







j& 



Fig. 3i 



Fig. 3a. 



Stabilité de Véquilibre des ûgures trouvées. — Appliquons 
ceci à notre problème et considérons la figure 21; ce que nous 
avons appelé \k représentera ici ce que nous avons appelé X et 
(0 correspondra à x^. 

La fonction F, qu'il faudra rendre maximum est 



W 



2J 



(i>^ augmente lorsqu'on se déplace de A sur D, et décroît de 
D à H; sur AH, M va en croissant, et croît aussi de B à C ou 
à C. 

Il s'agit maintenant de traiter la question de stabilité. 

Partons de A, où w :=: o et où on a une sphère, nous savons 
que dans ce cas l'équilibre est stable : donc tous les coefficients 
de stabilité sont négatifs. Ils cesseront d'être tous négatifs quand 
on rencontrera une figure de bifurcation. 

La première bifurcation que l'on rencontrera en suivant AH 
sera au point B, où Tellipsoïde de Mac Laurin diffère peu d'un 
ellipsoïde de Jacobi. Le coefficient de stabilité s'annulera donc 
en même temps que 



R.S. 



R*s, 



Après cela il ne s'annulera plus ; quand on arrivera en E^ un 
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autre coefficient s'annulera en changeant de signe, mais le 
premier coefficient restera positif. Donc de A à B nous avons un 
équilibre stable, au delà de B il est instable. 

Sur la branche BC, on a au contraire des ellipsoïdes qui sont 
stables jusqu'en G ; en G l'un des coefficients s'annule et au delà 
reste positif, donc à partir de G sur la ligne BC on aura des 
positions d'équilibre instable, mais à partir de G commenceront 
de nouvelles figures d'équilibre et ce sont celles dont il faudra 
étudier la stabilité ; il y à une difficulté, car, pour ces figures, 
M peut croître ou décroître. 

Il est aisé de voir que l'on a en tous cas un maximum ou un 
minimum ; on a en effet pour ces figures, si elles correspondent 
à des formes peu différentes de celle que l'on a en G 

en G même, e:=:o. Si on change e de signe, la figure a une 
forme symétrique de la première, M garde donc sa valeur ; il en 
résulte qu'en G la quantité M passe par un maximum ou un 
minimum. Si M<M(jOn a après une seule figure d'équilibre; si, 
au contraire, M > M^ on aurait trois figures d'équilibre dont une est 
stable. On n'a pas encore fait le calcul complet, bien qu'il semble 
plus probable que l'on a un minimum. 

Conclusion. — Quoi qu'il en soit voyons les deux hypothèses : 
supposons une masse fluide homogène animée d'un mouvement 
de rotation et soumise à un refroidissement assez lent pour 
qu'il ait le temps de se propager dans la masse ; supposons de 
plus que le mouvement de rotation soit le même dans toute la 
masse, et cela en supposant le frottement suffisamment grand. 

La vitesse du mouvement de rotation n'est pas constante, ce 
qui est constant c'est ul = (o^J, en vertu du théorème des aires ; le 
volume T doit diminuer puisque la masse se refroidit, donc M 
augmente. On a des figures d'équilibre stable jusqu'en B ; après 
cela, ce qui devient stable ce sont les ellipsoïdes de Jacobi et cela 
jusqu'en G ; le refroidissement continuant, on ne peut plus avoir 
d'ellipsoïdes de Jacobi, puisqu'ils sont instables. Si, pour l'autre 
série, Mfl est un maximum, cette seconde série ne donnera aucune 
figure, et le corps se dissipera dans l'espace; si, au contraire. 
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ATTRACTION DES ELLIPSOÏDES 



Mq est un minimum, on aura une sorte d'ovale étranglée que 
l'étranglement divise en deux parties inégales. 

On n'a pas fait le calcul jusqu'au bout pour voir ce qui 
arriverait ensuite : peut-être l'étranglement s'accentuera-t-il et 
la masse pourra se diviser en deux parties. Remarquons que ceci 
ne peut s'appliquer directement à l'hypothèse de Laplace, car 
la nébuleuse considérée par Laplace n'est pas homogène, elle 
doit au contraire être fortement condensée au centre. 

Quoi qu'il en soit, on ne sait pas jusqu'à présent si en G corres- 
pond pour M un maximum ou un minimum. 



CHAPITRE VIII 
ANNEAU DE SATURNE 



Etude des ûgures annulaires d'équilibre. — On peut faire 
trois hypothèses sur Tanneau de Saturne : on peut admettre qu'il 
est solide, qu'il est fluide ou qu'il est composé d'une masse de 
solides indépendants ; la troisième hypothèse connue sous le nom 
d'hypothèse de Cassini paraît être l'hypothèse véritable. 

HYPOTHÈSE DE l'aNNEAU SOLIDE 

Historique. — Laplace est le premier qui ait examiné cette 
hypothèse : et il a remarqué que si l'anneau était homogène son 
mouvement serait instable ; pour peu que son centre s'écartât 
du centre de gravité de la planète, Tanneau se précipiterait sur 
celle-ci ; l'attraction de Saturne ne tendrait qu'à accentuer l'ano- 
m.alie. Si donc l'anneau est solide il doit présenter des irrégula- 
rités ; Laplace ne chercha pas à déterminer l'ordre de grandeur 
de ces irrégularités ; Maxwell a fait le calcul et a reconnu 
qu'elles doivent être considérables ; ces irrégularités sont inad- 
missibles. 

Equations du mouvement. — Parmi les forces qui sont appli- 
quées à l'anneau il n'y a pas à s'occuper de l'attraction de l'anneau 
sur lui-même parce que l'anneau étant solide ses diverses parties 
n'éprouvent pas de déplacement relatif; il n'y a donc à s'occuper 
que de l'attraction de Saturne et de celle des autres corps 
célestes ; nous négligerons cette dernière action et nous allons 
nous occuper seulement de l'attraction de Saturne. 

Nous allons supposer que l'anneau a son centre de figure au 
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centre de la planète : l'observation montre qu'il en est ainsi; 
le centre de gravité G de l'anneau est supposé dijfférent de O. 
Un point M de l'anneau, puisque l'on ne parle pas de Tattrac- 
tion de Tanneau sur lui-même, est soumis à l'action de l'attraction 
de Saturne ; l'ensemble des forces d'attraction fait donc équi- 
libre à la force d'inertie. Un point M de l'anneau est soumis à 
l'attraction de Saturne et à la force centrifuge : si l'anneau con- 
serve son centre de figure au même point cette force centrifuge 
et cette attraction sont constantes : M tournera alors avec une 
vitesse égale à celle d'un satellite placé à la même distance ; 
c'est la condition nécessaire et suffisante de l'équilibre. Si 
l'anneau est homogène, cette condition n'est pas nécessaire, car 
la vitesse de rotation peut être quelconque : chacun des points de 
l'anneau serait soumis à la même force, et l'équilibre aurait lieu 
quelle que soit la vitesse de rotation. Dans le cas général le 
centre de gravité G décrira donc d'un mouvement uniforme une 
circonférence autour de Saturne. 

Prenons comme unité de longueur le rayon de l'anneau, 
comme unité de masse la masse de l'anneau ; son moment d'inertie 
par rapport au centre de figure sera égal à i . Il reste à déter- 
miner l'unité de temps : on peut la prendre de façon que, en 
désignant par M la masse de Saturne et par /*le coefficient de 
la gravitation, on ait /'M= i. 

Etudions le mouvement : soient x q\ y les coordonnées du cen- 
tre de gravité, nous avons 

f=r I pcoss.rfs. 



X j prfs= I p 

y I pds = I p sm s. ds^ 



en appelant p la densité de l'élément linéaire de l'arc ds de 
l'anneau supposé étroit. 



çds 
est la masse de Tanneau, 



/• 



STABILITÉ DU MOUVEMENT 
Les équations du mouvement du centre de gravité sont 

-jj- j p«Js=s I p C08 s ds^^x j pds, 
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d^où 

Les intégrales de ces équations sont, en prenant une origine, 
des temps convenable, 

^ = ot cos ^, 
y = P sin t. 

Le point G décrit donc une ellipse dans le cas général ; si 
c'est une circonférence, elle est décrite avec une vitesse angulaire 
égale à Tunité : c'est la vitesse d'un satellite placé a l'unité de 
distance et décrivant un cercle et dans ce cas, a = p. 

Stabilité du mouvement. — Il reste à voir si ce mouvement 
est stable. Supposant négligeable, par rapport à la masse de 
Saturne, la masse de 
l'anneau, nous pou- 
vons supposer Sa- 
turne fixe. 

Par son centre O, 
faisons passer deux 
axés fixes OX et OY 
(fig. 33) ; par le cen- 
tre de figure de l'an- 
neau, faisons passer 
deux axes ox, oy in- 
variablement liés à l'anneau, leur position, et par suite celle 
de l'anneau, sera déterminée, si l'on se donne les coordonnées 
X et y de par rapport à xoy^ et l'angle de l'axe ox avec OX. 

Soit G le centre de gravité de l'anneau, qu'on peut supposer 




Fig. 33. 
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sur ox ; on a oG = a, c'est une donnée du problème ; les coor- 
données X et Y du centre de gravité sont 

X = — y sin + (a — x) cos 9, 
Y = y cos 9 + (a — x) sin 9. 

La force vive de Tanneau est égale à la force vive de trans- 
it 
lation (X'^ + Y'^) augmentée de la force vive de rotation JO'*. 

Le moment d'inertie par rapport au point o est donné par la 
relation I :=: J + a^. La force vive est donc une fonction de x, y, ft 
et de leurs dérivées et on a 

X'^ + Y'^ = y^ + y' '-h r [y' + (a - x)'] - a9' [x'y + (a —x) y']. 

L'énergie potentielle est 

U =— /MV, 

V désignant le potentiel de l'attraction de l'anneau au centre de 
gravité, et comme on a 

on a 

U=— V. 

Supposons que la densité de l'anneau soit une fonction déve- 
loppable en série de Fourier. Si d^ est l'angle sous lequel on 
voit du centre de l'anneau un élément linéaire, on aura, en appe- 
lant p la densité au point déterminé par l'angle ^, 

P = ?o| " + 2«cos^+-^cos2^{/ + -^sin 2^ + ... 1; 

si i et 71 sont les coordonnées d'un point de l'anneau par rap- 
port aux axes mobiles, on a 



I pds=i, 



/' 
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On trouve encore 

r r- 

I p$7) rfs = I psîn2^rf^ = ~-. 

Le potentiel en G de l'attraction de Saturne sur l'anneau 
est, comme on le voit en effectuant le calcul, 



/ 



P''* v_ 



V = — U 



Les équations de Lagrange donnent 

d /dT\ dT . dU 



Q- 



:0, 



dt \ ^q' / ^q àq 

mais comme U ne dépend pas de y' on a 

d Ô(T — U) Ô(T — U) __ 
dt bq' dq 

On peut écrire 

où l'on a posé pour abréger 

A = J+y^+(a-.r)S 
B,= — x'y + y' (a — or) , 
C=x"+y" ; 

PoiNCARÉ. — Figures d'équilibre. 
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remarquons que T — U ne dépend pas de 8 ; la première des 
équations de Lagrange donne donc : 



j^r a(T-u) i 

dt L ôO' J ' 



et, comme U ne dépend pas de 0, il reste 

p étant une constante. Il est d'ailleurs facile de calculer cette 
constante, étant données les conditions initiales ; si on déve- 
loppe l'équation on a 

Ae'+B=/?. 

Supposons x^:^ 0,^0=: o, et par suite x'^ et y'^ très petits, on 
trouve 

Bo = o, 

et par suite p=z^'^\ donc c'est la vitesse de l'anneau, supposé en 
équilibre : on a, comme on l'a démontré plus haut, 0'^, = i ; 
donc p=i : l'équation obtenue est l'équation des aires. 

Transformons maintenant les équations de Lagrange : posons 

H = T — U — /?6'. 

Les équations de Lagrange conserveront la même forme, si on 
y remplace T — U par H ; car, si q est l'une des variables, on a 

5H _ 5(T — U) 5T 59^ 56^ 

dy ~ 0(7 "^ dO' ôy ^ dq ' 



et, comme on a 


ÔT 


on conclut : 


dH_ d(T— U) 
dq dq 
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De même on voit que Ton a 

OH __ 5(T— U) 

On peut faire disparaître 8' de ces équations, on n'a en effet qu'à 
remplacer 9' par 

/>-B . 
A ' 

et comme, ni T ni U ne contiennent Q, H ne sera plus fonction 
que des deux variables x et y, et de leurs dérivées. On trouve 

2A 2 

Supposons maintenant que le centre de l'anneau soit très près 
de Saturne, c'est-à-dire supposons j: et y petits, et voyons si dans 
le mouvement ultérieur x et y resteront petits. 

Les termes de degré zéro qui entrent dans H peuvent être 
négligés, car ils n'interviennent pas dans les équations de 
Lagrange, Les termes du premier degré en ^ et y doivent dispa- 
raître d'eux-mêmes : en effet, dans les équations de Lagrange ils 
donneraient des termes constants, et ces équations doivent 
admettre la solution ar = o, y = o. Quant aux termes du premier 
degré en a/ et y^ ils disparaissent dans les équations de Lagrange. 
Il ne reste donc que les termes du second degré, que nous 
appellerons W, le reste étant négligeable. 

Développons H, en ne conservant que les termes du second 
degré ; on trouve pour W le développement suivant : 

Ex' + 2Cry+Fy' + — + -î^^ y"-x'y +(i - 2a')xy 



en posant 



2 



u 3 2 ■ ** 



4 



C_-, 



4 ^ 4 



M 
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Les équations de Lagrange, comme on Ta remarqué, peuvent 
se réduire à 

d 5W ^ ÔW 

dt àx^ ~ àx ' 

d 5W aW , 
dt ay ~ ày ' 
c'est-à-dire 

(i _ a^) y + 2 ( I — a}) x'= 2Gx + 2Fy . 

Ce sont deux équations linéaires à coefficients constants ; pour 
intégrer on pose 

On aura quatre valeurs possibles de w, et l'équation générale 
du mouvement sera 

x=^l^ e«-*' + ^2 e'"^' + ^8 e»'*»' + Z^e'-*' , 

les l et les m étant des constantes quelconques déterminées parles 
conditions initiales. Pour que le mouvement soit stable, quelles 
que soient ces conditions, il faut que les valeurs de co soient à 
partie réelle négative ou nulle. 

Pour déterminer w, on remplace a: et y dans les équations 
diflférentielles et on a 

Kù^x — 2 { I — a^) (i>y = 2 Ea: -|- 2 Cy, 
(i — a^) (o^y+2(i — a?) (oa: = 2Cj7+2Fy, 

on peut éliminer x ely entre ces deux équations, et on trouve 

co*(i — a^) + a>«(^— iH -^ 1 ^j 

4.J-(9_a2_p2_24^2^8a^a)=o. 

L'équation en co^ doit avoir ses racines réelles : car si elles 
étaient imaginaires, les racines carrées de ces deux racines ne 
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seraient pas des imaginaires pures, et par suite, comme elles 
sont égales et de signes contraires, Tune serait à partie réelle 
positive. 

Or, pour un anneau homogène, où 

a = a= P = o, 

Téquation se réduit à. 

/ »* , .2 I 9 -, . 

les racines en sont imaginaires : si a, a, ^ sont voisins de 
zéro, elles seront encore imaginaires. Il faut donc que a, a et ^ 
soient grands, il y aurait donc de grandes irrégularités dans 
l'anneau : et sa régularité apparente serait inexplicable. 

Maxw^ell a fait le calcul dans Thypothèse où l'anneau serait 
homogène, sauf en un point, où il y aurait une masse supplé- 
mentaire : et il a trouvé que cette masse devrait être au moins 

les -^ de la masse totale de Tanneau, c'est-à-dire que l'on aurait 

a>o,8. 

M. Radau a examiné le cas où l'épaisseur de l'anneau serait 
variable d'un point à un autre, et il a trouvé que la densité 
devrait varier de 2,7 à o,o4, ce qui ne parait pas probable. 

De plus, la section de l'anneau étant mince, il lui faudrait 
une grande résistance pour résister a l'attraction des satellites : 
Hirn a calculé qu'il lui faudrait une résistance mille fois supé- 
rieure à celle de l'acier : donc Vanneau de Saturne n'est pas 
solide. 



HYPOTHESE DE L ANNEAU LIQUIDE 

Potentiel d'une circonférence homogène. — ChercRons le 
potentiel d'une circonférence homogène, en un point P de l'es- 
pace. Soient PA = a, PB = b les distances minima et maxima 
du point P à la circonférence. La circonférence étant donnée 
de rayon Rj aèt b déterminent complètement le potentiel. 
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Soitco Tangle AOM (fig. 34), on a 



PM* = a' cos^ — + *' sin^ —, 




le potentiel 






Fig. 34. 



="/"■ 



=^M^ia,b). 



dtù 



2itXPM 



Comme le potentiel ne dépend que de a, de b et de la masse 
de la circonférence, on peut ramener 

le calcul de V au potentiel d'une cir- ^^ ^û 

conférence, en un point de son plan : 
on décrit une circonférence de dia- 
mètre [a -f- ^), et on calcule le po- et 
tentiel de cette circonférence en un 
point situé à une distance a d*une des 
extrémités du diamètre (fig. 35). 

Le potentiel en F est le même, 
on a 

diù 




Fig. 35. 



'="/ 



2TU 



\/' 



a^ cos^ — + P sin^ — 

2 2 



mais on a aussi 



V= 



Vo ^^^^^ 



4 



Q étaid; le point de la circonférence dont la projection sur le 
diamètre AB est en P, <{/ étant l'angle AOM ; on a donc 

<p{a,i) = <p(OA,QP) 

=<p(-^'v^)- 
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Si donc on forme une série de quantités 

a b 

a, b. 






telles que Ton ait 



on aura 



.,= ^^ 



et cela quel que soit Tentier n. On démontre facilement que a„ 
et in tendent vers une même limite qu'on appelle la moyenne 
arithmético-géométrique de a et de b. 
Soit m cette limite, on a 



V = M 



271 W m^ cos^ -i- + w^ sin* -^ J o 



«p (a, b) est une fonction homogène de a et de i, de 
degré — i , c'est-à-dire que Ton a 

<p().a,Xi) = X-«^(a,4); 

on peut donc écrire 

?(«*)=xt(|-, ') 

et, remplaçant «p (a, i) par «p (a), on a 

î(«.*)=|?(y)- 

• Que devient la formule si P est rapproché de la circonférence? 
a est alors une quantité très petite. On sait que l'on a, en dési 
gnant par L le logarithme népérien 

P 
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p désignant' la distance du point P au point P^ le plus voisin du 
cercle, e étant une quantité tendant vers zéro avec p, jx,, et K étant 
des constantes ne dépendant pas de la position du point P par 
rapport à .P^ non plus que de la position du point P^ sur le cercle. 
En reprenant les notations précédentes on a 





V_a(x,L^; 


et comme 






M = anfXjR, 


on a 






M ttR a ' 



b diffère peu de 2R, à Tordre d'approximation adoptée ; on a 
d'ailleurs 

— = cp(a,i). 
Cejte relation permet de calculer la constante K : on a en efiet 

,(,.)=,(v/a?7:î±i-)=^,(^). 

si on fait ( = i on a 

'■ 2 / i\/a \ 

?w= — i — ? ( — I — )• 

11 vient donc dans le cas qui nous occupe, en négligeant a 
devant l'unité 

(p(«) = 2cp(2v/a); 



mais on a aussi 



ç(a) = L — ; 



donc 



\ « / \ 2v/a / 

LK — La ^ 2 LK — 2 log 2 — La ; 
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d'où on conclut 

K = 4, 

et finalement 

M Tz \a/ 

avec l'approximation donnée. 

On peut chercher une forme plus exacte pour '^ (a). Posons 

+ ^[Bo+B.«+BX + ...]; 

en se servant de la relation 

f \ ^ ( V^ \ 

cp(a) = ; ol — ^^ ); 

^^ ^ i+a • \ i+a/ 

on trouve facilement que l'on a 



A, = o, 


B,=B, = o, 


A — -i 


B,=-^L., 



de façon que l'on peut prendre pour îp (a) la fonction 

4 



^-«(4) 



et cela en négligeant seulement les termes en — ^ , 

Masse annulaire de révolution, — Soit G le centre de gra- 
vité de la section méridienne, menons Gy parallèle à Taxe de 
révolution et Gx dans le plan méridien , perpendiculaire à G// 
et tel que le sens positif soit dirigé de l'axe de révolution vers 
le point G, soit l la distance du point G à cet axe de révolution. 
Un élément de surface rfo-' de la section méridienne, de densité p' 
et de coordonnées x' et y' engendre un anneau dont la masse est 
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Cherchons le potentiel de cet anneau en un point P du plan 
méridien, dont les coordonnées soient x et y; il faut d'abord 
calculer les quantités a et b. Si on suppose le point P du même 
côté de l'axe de révolution que l'élément rfo-, on a 

Supposons l'anneau mince et le point P situé près de G, de 
façon que x^ y, x'^ y' soient négligeables devant l; a est une 
quantité du même ordre ; on peut prendre pour valeur approchée 
de b 

car (y — y'y est négligeable ; il faut maintenant calculer y (a, b). 
On a 

?(«.*)=t?(t) 

^l-\-X'\-x' a 
Le potentiel en P est donc 

V^ ^9'd^{l-hx') ^ /jb . 
b a 

on a 

b 2^ + ^ + 0/ 

comme par hypothèse x et j?' sont petits, on peut écrire approxi- 
mativement 

l+x' 



Il vient ensuite 



I r x'—x'\ 



a a ^ 2I 



MASSE ANNULAIRE DE RÉVOLUTION 187 

ce qui peut s'écrîre, avec le degré d'approximation donné 



2,1 2I 

on peut alors écrire aussi 

Le potentiel de l'anneau sera donc 

les autres termes étant négligeables, et les intégrales étant 
étendues à la surface de la section. Nous allons supposer l'anneau 
homogène et nous écrirons p' = i. 

Les deux premiers termes sont des potentiels logarithmiques 
qui jouent, dans l'attraction des masses planes, le même rôle que 
le potentiel newtonien dans l'attraction des corps célestes ; les 
propriétés du potentiel logarithmique dans le plan sont anulo- 
gues à celles du potentiel newtonien dans l'espace : on démontre, 
par exemple, que le potentiel logarithmique d'un cercle homo- 
gène est le même que si toute sa masse était concentrée au centre. 

Si le point P est à la surface d'un tore, à une distance a^ du 
centre du cercle générateur, il suffit de considérer le potentiel 
logarithmique en P du cercle de rayon a^ ; sa surface est wa*» 
Son potentiel est donc 

„ 81 

Tzal L — ; 

la première intégrale est par suite 

2Tzai L 

La seconde est la somme de deux autres : 

fi/ 



lia fia 
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dont la seconde se calcule immédiatement, on trouve pour sa 
valeur 

X ,^ 81 

la première va se calculer plus tard. 
On trouve de même 

puisque le centre de gravité est à l'origine des coordonnées. 
Il reste à calculer 



/ 






nous allons employer un procédé de calcul dont nous nous 
sommes déjà servi dans une autre question : supposons que la 
densité du cercle soit 



le potentiel sera 



2 2 



lJ=J,'d^L^; 



la densité ne dépendant que de la distance au centre, le potentiel 
logarithmique du cercle est le même que si toute la masse est 
concentrée au centre. Cette masse est facile à calculer, c'est : 



^^d^^l d, ^ 



rdr 
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r <r' r'I 



Le potentiel est donc, en un point Â^ situé à une distance a^ du 
centre du cercle a. > a. 



09 



TZUo 



SI 



Supposons que le cercle C se déplace d'une quantité e parallè- 
lement à l'axe des x, de façon à 
occuper la position C^; le poten- 
tiel en A sera (fig. 36) 

5U ^ 



U- 



la différence de potentiel sera donc 
dU 




Fig. 36. 



Nous allons étudier d'où provient cette différence : en un 
point de l'intérieur du cercle la densité était p', elle est devenue 



p'-e 






Dans la région du cercle C extérieure au cercle C, la densité 
était nulle, elle est devenue égale à p' qui est petit puisque la 
densité p' est nulle sur la circonférence du cercle C, et l'aire de 
la région considérée est de l'ordre de e. Il y a encore à consi- 
dérer l'aire égale à la précédente, intérieure au cercle C et exté- 
rieure au cercle C, la variation de densité en sens contraire de 
la précédente est aussi de l'ordre de e au plus, puisqu'en un 
point de cette aire la densité p', petite sur le bord du cercle, est 
devenue nulle. 

La variation de potentiel due à ces deux aires est donc négli- 
geable, et on a 

dU 
àx 



I Ox' a 
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Or 

il vient donc, en sapprimant le facteur e, 



/ 



a ox 



L'intégrale que nous voulons calculer est donc, au facteur/ près, 

égale à -^ — ; or - — est la composante de l'attraction suivant 

• I • X 

Qx ; cette attraction est , sa composante suivant 0:r est — r-- 

«I ai 



On a donc 



/ 



/ a ^l a^ ' 



et, si on suppose le point situé sur la surface du tore, c'est-à-dire 
si l'on a a^ = a^^ il vient pour la valeur de l'intégrale 



2 ^ 
'^al- 



On a donc pour le potentiel en un point de la surface du tore 

\T 2 T 8^ X ^ 81 ^ 5X 2 

Y = 2Tzaî L y Tzal L 1 77- Tta!, 

^0 ^ ^0 4^ 



en négligeant les quantités de l'ordre de 



P 



Section méridienne de l'anneau. — Si l'on appelle M la 
masse de Saturne, (o la vitesse de rotation de l'anneau, une 
condition nécessaire de l'équilibre est que, à la surface du 

tore, la fonction U -) soit une constante, U représentant 

le potentiel de l'anneau sur lui-même augmenté du potentiel de 
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Saturne, ce dernier peut être représenté, avec l'approximation 
employée, par—— et Ton doit avoir à la surface de l'anneau : 

v+"+J^=c-. 

L 2, 

En posant 

la distance d'un point de la surface au centre du tore est 
ce qui, avec l'approximation que l'on cherche, peut s'écrire 

d'où 
enfin on a 



I X 



p l+x l P ' 

La condition d'équilibre peut s'écrire 

2,Tzal L =- Tzal L h -77- izai 

il faut donc que le coefficient de x soit nul, ce qui donne 
DTraJ Tzal , 8Z M 



L «-+ w^^ = o. 



4^ la 

Cette équation donne w'^, aux quantités du deuxième ordre 
près. 

Comme -y- est aussi petit, on a 

2 M 
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La vitesse de rotation d'un tore engendré par un cercle, de rayon 
suffisamment petit par rapport à sa distance de Saturne, est 
donc la même que celle d'un satellite placé à la même distance. 
S'il n'y a pas de corps central, M = o : on voit alors qu'en 
l'absence de corps central il peut y avoir une figure annulaire 



(t) 



d'équilibre : a^ est alors très petit, mais le rapport tend vers 

l'infini quand a^ tend vers zéro, car L — augmente alors 

indéfiniment. 

Si on voulait pousser plus loin l'approximation, il faudrait 
calculer les termes de la forme de 



■/ 



QJ 



et 



x"^ y"" Ç x"^' y"" dfs 



Supposons maintenant que le tore ait une section elliptique ; 
nous devons alors calculer le potentiel logarithmique d'une 
ellipse. Or le potentiel newtonien d'une droite indéfinie en un 
point, est le potentiel logarithmique en ce point de la projection 
du point sur la droite ; il en résulte que le potentiel logarith- 
mique d'une ellipse en un point de son plan est le potentiel 
newtonien d'un cylindre long par rapport à sa section droite, 
en un point situé k une petite distance de sa surface. On est 
donc amené à calculer le potentiel d'un ellipsoïde très allongé 
pour lequel on supposerait c^ infini. 

Le potentiel intérieur est donc une fonction du second degré 

dont les coefficients sont faciles à calculer (page i4o)- 
Pour le potentiel extérieur on a 
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On peut calculer A^, A^, Ag comme on Ta fait précédemment, 
mais le^ intégrales devront être étendues à la surface d'une 
ellipse dont l'équation est 

2 -2 

•^ _!. '^ — , 

Les quantités du troisième ordre sont négligeables. La fonction 
de force peut de même se développer en série de la forme 

La condition d'équilibre est que l'on ait 

(A. + C.) a- + (A,+C,)^^ + (A, + C,)^' = C» 
à la surface de l'ellipse ; il faut donc que l'on ait 

Aj + Cj = o 

(A,+ C,)a^ = (A3+C,)i^ 

Laplace avait fait le calcul, et a discuté l'équation mais sans 
tenir compte des termes en A^ et Ag. 

^jme jg Kowalewski a été plus loin et a fait voir que la section 
méridienne du tore n'était pas symétrique par rapport à l'axe 
des 1/. 

Nous étudierons plus tard là question de la stabilité ; nous 
allons d'abord discuter l'hypothèse de Cassini. 

HYPOTHÈSE DE CASSINI 

Vraisemblance de l'hypothèse. — Cette hypothèse est la seule 
compatible avec les travaux de Maxwell sur l'anneau de Saturne. 
L'observation confirme cette hypothèse, car l'anneau est trans- 
parent et la lumière le traverse sans trace de réfraction ; il faut 
donc supposer qu'il est composé de corpuscules solides (ou 
liquides), isolées les unes des autres. Les observations spectros- 
copiques montrent d'ailleurs que la vitesse d'une molécule de 
l'anneau n'est pas la même sur le bord interne, ou sur le bord 
externe* 

PoiNCARÉ. — Figures d'équilibre. i3 
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Equations du mouvement — Nous allons commencer par 
étudier un cas très simple : nous allons considérer p satellites 
de même masse (jl, équidistants sur un cercle ; la distance 

de deux satellites sur ce cercle sera = 26 ; le rapport -^- 

n'est pas infiniment grand. 

Un mouvement possible du système est celui où chaque satel- 
lite parcourrait le cercle d'un mouvement uniforme, avec une 
vitesse w, déterminée par l'attraction de la planète, augmentée 
d'un terhie où il faudrait tenir compte de la masse des autres 
satellites ; et nous allons chercher s'il y a un mouvement peu 
différent de celui-là. 

Soit Mfi le satellite de rang /c , son rayon vecteur sera i + p,„ 
son angle polaire sera wi + 2AO -f- o*^; dans le mouvement non 
troublé on aurait 



(7^== o ; 

si le mouvement est stable p^ ^^î* rester petit. 

La force vive T, et la fonction de forces U, sont fonctions de C;^ 
et de (j/,; les équations de Lagrange donnent 



rf /dT \ 

dt U^J" 

^ / 5T \ 5: 

dt [ dp', ; de 



ôT dU _ 

H "^ H -^' 
ôT dU 



en posant 



on a 



dt \ dp',, / dp,, dp,. 



T — U== H, 



m m 



dt dp', dp 

Nous pouvons supposer que Saturne est fixe, c'est-à-dire que 
nous négligerons les perturbations apportées dans le mouvement 
de Saturne par celui des satellites ; la mise en équation du pro- 
blème est alors indépendante du mouvement de Saturne. 
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La force vive du système est 

et, en négligeant les infiniment petits du troisième ordre : 
La fonction de forces est 

R représentant la partie complémentaire du potentiel due à 
l'attraction des satellites les uns sur les autres ; si l'on pose : 

on a 

+ (?.+P.) '-^ cotg<,+ i2î=^colg+ 



K— 'a) 



(i4-2cotg*.|) +...J, 



8 
les termes non écrits étant négligeables. On peut écrire aussi 

— J— =!x(t— pH-py. 

Les termes les plus importants de R sont évidemment ceux 
pour lesquels ^ est petit, sin ij étant un infiniment petit du pre- 
mier ordre, — . ^, est un infiniment £rrand du second, et . ", 
sin^I; ^ sin^j; 

est un infiniment grand du troisième. 

U et T comme on le voit ne dépendent pas de /, et sont du second 

, dT dT 
degré par rapport a p, , or et t ; — — , - ^ sont du premier 
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degré par rapport a ces quantités; enfin —. — ( -^-p j est du pre- 
mier degré avec des dérivées du second ordre. Les équations de 
Lagrange sont donc des équations à coefficients constants de la 
forme 

A + B(p,p',p";.y,0=o, 

A étant une constante, et B étant une fonction linéaire des argu- 
ments ; mais ces équations doivent admettre la solution p = <t = o; 
il faut donc que A disparaisse de lui-même ; les équations sont 
alors des équations linéaires et homogènes, qui s'intègrent par 
les exponentielles, et on a, 



-2 



ae^' 



les X étant des racines d'équations algébriques, et les a étant des 
constantes. 

Le mouvement sera stable si p ne croît pas au delà de toute 
limite ; en mettant en évidence la partie réelle et la partie imagi- 
naire de A, on a 

le mouvement sera stable si \ est négatif ou nul ; mais si nous 
démontrons que l'équation en X a ses racines deux à deux égales 
et de signes contraires, nous verrons que \nQ peut être que nuU 
Or si nous changeons Z en — t les équations ne changent pas, 
puisque les termes du second ordre seuls interviennent dans les 
équations de Lagrange ; il en résulte que si Téquation en \ 
admet une racine, elle admettra aussi la racine de signe con- 
traire, et par suite, la condition nécessaire de la stabilité de 
l'équilibre, est que les racines de l'équation en \ soient purement 
imaginaires. La condition est d'ailleurs évidemment suffisante 
puisque, dans ces conditions, p s'exprime par des fonctions trigo- 
nométriques de la variable réelle t. 

Vous avons 2/? équations du second ordre, donc \ satisfait à 
une équation algébrique d'ordre ^p ; c'est une équation en A^, de 
degré 2/> dont toutes les racines sont réelles et négatives. 
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Les premières des équations de Lagrange sont 

pi-— <^'(ï+?a) — 2W<J\.+ I — 2pjt= [7.^-^j ; 

on peut écrire, en mettant en évidence dans R les termes de 
degrés o, i et 2 

R = R, + Ri + R,; 



et Ton a 



ùR _ ùR, 



dR, 



ôR, dR, r • ,. ^ . 

^ estune constante. ^ est une fonction linéaire. 

Les secondes équations sont 

dR, dR, 

Dans la première équation les termes constants doivent dispa-> 
raître : on a donc 

OR. 



to' — I 



' dp, 
4 ^ sin»} ' 



tj^ prend les valeurs 6, 2O, 39, ..., [p — i) 0. Dans le cas ou [jl=o 
on aurait (o'^= i. 

Les équations sont donc 



p^' — 263<ï;. — ((0^+2)p,= p 



= 1^ 



5R, 
ôR 

05. 



Changement de variables, — Faisons le changement de 
variables, 
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On a par exemple 
Posons de même 



^k=^f\,e'^^iK 



Les \.^ d'une part et les 7).^ de l'autre sont deux à deux imagi- 
naires conjuguées ; d'une façon plus précise je dis que Ç^ et ^p^^ 
sont imaginaires conjuguées. En effet, on peut écrire 



p.=2i?,_./'*(P-T)«=2$,_,,. 



Cette seconde valeur de p/. se déduit de la première en chan- 
geant i en — i et Ç.^ en Çp— ;: il faut donc que $p_.^se déduise de ?., 
par le changement de / en — ï ; donc Ç^ et Çp_.^ sont imaginaires 
conjugués. On démontrerait de même que vj^ et fip^^ sont imagi- 
naires conjugués. Si p est pair, \p est réel, car il est à lui-même 

ï 
son imaginaire conjugué. 

("4^) sera un polynôme du second degré en \\^ qui ne ren- 
fermera d'ailleurs que les termes de la forme Ç'^ i';,-v. En effet, T 
et R sont symétriques par rapport aux p/» et aux o*/^. ; si Ton per- 
mute circulairement les p;;^^ d'une part et les o*^ d'autre part, c'est- 
à-dire si l'on remplace i. par Ç^ e^'r*, \'^ est multiplié par le même 
facteur. 

De même pour tj^ et 'f\\ . Par ce changement ni T ni R ne doivent 
être modifiés : il faut donc qu'un terme Ç'^ 5V ^e transforme en 
lui-même ; donc il faut que le changement de variable multiplie 
ce facteur par e-*'*^' ; il faut donc que 

2/ (y + y') 9 = ^lihp ; - 
donc 



CHANGEMENT DE VARJABLES 199 

En appliquant toujours le même raisonnement, on voit qu'en 
définitive T et R ne doivent contenir que des termes de la forme^ 



i'v^'p-. 


/ / 

U: f\p- Y 


Çy Çp-y 


IfWp-'i 


Çy kp - Y 


il 'f\p ~ Y 


Ç Y '^p - Y 


'^iY ''ip - Y 



Si maintenant on regarde de plus près la formation de la fonc- 
tion II, on voit que Ton peut poser 

ii==ii, + ii, + ... + n, + ...+ n, 

en écrivant II., à la place de 

i; iV _ .,+ 2(0 (i. 'Vp _ .,+ i^ ^ .;Vy) + 0^1, ip _ ..+'Vy 'V/> - Y+ 2^T?P - V 

+ ;^ I — L.^ÇY^P-Y + ^*T(^Y*'^/>-Y + '^Y^P-ï) + ^T'^»Y*'ip-Y T 

L, M, N sont donnés par les formules 

VI / sin^Y^]^ cos^^]^ cos"Y^ \ 

'* ^J \ 4 sin**^ 2 siu ^ / ' 

„ VI sin 2Y^ cos ^ 

^ VI /sin^y^'^os^^ sliryi \ 

' *^ ^ \ 2 sin^i]; 4 sin ^ / ' 

Les équations de Lagrange seront alors 



dt 






et on aura en supprimant les indices 

i" _ Ç ((o^ + 2) - 2co-/i'= ;. [LS + Mr,], 
V'+2«^' =[iL[M$ + Nr.]. 

L'équation en X est 

().2 _ 0,^ — 2 — I-IX) ()>* + Nu) — ([il*M* — 4to*X^) = o. 



200 ANNEAU DE SATURNE 

Discussion des équations, — Quel que soit y, les racines de 
cette équation en X* doivent être réelles et négatives ; on voit 
que cette condition sera remplie si [jl est suffisamment petit; 

2,3 

Maxwell a démontré qu'il fallait que l'on eût [jl < \ — . Dans 

ces conditions ^ipiend vers zéro quand /? augmente indéfiniment; 
c'est là une difficulté, car la masse de l'anneau précisément égale 
à [x/?, serait négligeable ; mais cette difficulté disparait d'elle- 
même si on suppose les satellites répartis non sur une circon- 
férence, mais dans un certain volume. 

Si les racines de l'équation en A^ sont négatives, on pourra 
poser 

$., = Ae^' = Ae"»', 

n étant une quantité réelle, et en prenant 

A=A,e'«% 

on aura 

On aura une solution particulière si à i., r^^ on ajoute les solu- 
tions, 

T,. = 0. 

et on aura alors une solution particulière pour les p;^ en posant 

p/,=Aoe'("'-^^ + 2A--«). 

a cette solution on peut ajouter la solution imaginaire conjuguée; 
on obtiendra ainsi des solutions de la forme 

<7iç= BqCos (nt-\-o-\- aÀ-yÔ). 

Si on suppose un système d'axes entraînés constamment avec le 
satellite M^ dans son mouvement non troublé, et tel que l'un d'eux 
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coïncide avec le rayon OM,., p^^. et d,, seront les coordonnées du 
satellite, par rapport à ces axes, et on voit que le satellite décrit 
par rapport à eux une petite ellipse. 

On passe de la trajectoire d'un satellite à celle du suivant en 
changeant k en A -f- i ; le second satellite est donc en retard sur 
le premier, et ainsi de suite ; il y a ce que Ton appelle un déca- 
lage égal h yO. Si Ayô est un multiple de air le satellite de rang 
A -f- ï sera dans la même position que le premier satellite; celui 
du rang /: + 2 dans la même position que l,e second et ainsi de 
suite. 

Le mouvement sera comme s'il se transmettait d'un satellite à 
l'autre avec un certain retard; on peut aussi dire que le mouve- 
ment se transmet comme une onde, dont la longueur serait . 

Les ondes qui sont les plus dangereuses pour la stabilité 
sont les ondes courtes; les ondes qui compromettent le plus la sta- 
bilité sont donc celles qui correspondent aux grandes valeurs 
de y. Considérons en effet une valeur de n, suffisamment grande ; 
pour une grande valeur de y, l'onde serait très courte, deux satel- 
lites voisins pourraient se rapprocher d'une manière sensible, et 
leur action mutuelle ne serait plus négligeable par rapport à 
l'action de Saturne. 

2 3 
La masse totale ^p doit être inférieure à —V; si on suppose 

p=z 100, on voit que la masse de l'anneau devra être inférieure à 

2 3 

^ , c'est un inconvénient de la théorie de Maxwell ; mais 

10 000 

c'est un inconvénient artificiel, car l'hypothèse que nous avons 
faite est trop simple et il faut supposer une distribution de satel- 
lites occupant un certain volume de l'espace. 

La perturbation apportée au mouvement d'un satellite par l'ac- 

tion d'un autre est de l'ordre de —^ , en appelant S la distance 

des satellites. 

Si les satellites sont répartis sur une circonférence de 

longueur Z, on a ô = — : le mouvement sera stable si la per- 
turbation a une valeur moindre qu'un nombre donné A. Dans 
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rhypothèse de satellites répartis sur une circonférence, il faut 
donc que l'on ait 

ou bien 

Si les satellites sont répartis sur une couronne comprise entre 
deux cercles de distance comparable à la longueur de l'un deux, 
le nombre des satellites sera comparable à p^. S'ils sont dans un 
anneau de volume fini, leur nombre sera comparable à y/, la 
masse de l'anneau sera égale à [jl^' et la condition de stabilité 
sera que l'on ait 

la masse de l'anneau devra donc être inférieure à un nombre qui 
ne dépend pas du nombre des satellites. 

Stabilité de l'anneau supposé liquide. — Nous allons main- 
tenant, d'après Maxwell, examiner de plus près la stabilité du 
mouvement, dans le cas où l'on suppose l'anneau, soit fluide, 
soit composé de corpuscules indépendants. Dans cette partie du 
travail de Maxwell, les raisonnements ne sont pas tout à fait 
rigoureux, ni même très clairs, mais on doit en accepter la con- 
clusion; il est probable qu'un calcul rigoureux conduirait à des 
résultats peu différents de ceux auxquels nous arriverons. 

Soient or, y, z les coordonnées d'un point dans le mouvement 
normal; n, r, ii\ les composantes de sa vitesse ; 



v + l 


y + -^o 


!/ + K, 


, d^^ 


, dyi 


, 'i^ 


"+dt' 


^+ dt ' 


•*'+^'' 



les mêmes éléments dans le mouvement troublé. 
La force vive du système est 

-/[(«+4)'-(-#)'-(-f)] 



odt 
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réiiergie potentielle — U est égale à 



rix±.iii?+?irf., 



en appelant V et p le potentiel et la densité clans le mouvement 
normal, V + V, p + p', ce que deviennent ces quantités dans le 
mouvement troublé. V est une somme de quatre intégrales : 

fl?ÈL+ rVl^+ ClîÉL\^ flÉlL. 

La première intégrale est une constante ne dépendant pas 
de i, 7), Ç; la seconde et la troisième sont égales, comme on Ta 
vu au commencement du cours. Elles représentent la moitié de 
l'énergie potentielle due à l'attraction de la matière répartie, 
comme dans le mouvement normal, sur l'excès de matière qu'il y 
a en chaque point dans le mouvement troublé. 

On peut poser en négligeant le cube de i, yj, JJ, 



1 VI.. .yv v' ^ 



di/dz 



La somme de la seconde et de la troisième intégrale peut donc 
s'écrire 

Il faut maintenant calculer la dernière intégrale; je peux sup- 
poser que le déplacement J, r|, Ç est la somme de deux déplace- 
ments, l'un dans le sens de la propagation de ce mouvement, 
l'autre transversalement à ce sens : on posera ainsi : 
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On sait d'ailleurs, d'après les formules connues des vibra- 
tions longitudinales et transversales que 

çl^dx + ryr.^dy + ^X,^dz = dh 
est une différentielle exacte et que l'on a 

rf.r "•" f^y "•" rfc ~ 
On a d'ailleurs la relation 

o'+^MI+lM+ ^(p^) .....0. 

^ </:r é-Zy rfr 

Dans ces conditions on peut écrire 

Prenons maintenant des axes mobiles tournant avec une 
vitesse co autour de l'axe de l'anneau, l'expression de V ne 
change pas. T devient dans ces conditions 

La condition d'équilibre est que 

soit maximum ; si Ton tient compte de la relation 
pij dx-{- pr,^di/-i-p^^dz = d^. 
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on trouve les équations 






dP 59 

ôj; ^' d^ ' 



P est un polynôme du second degré en (i, 'f\, JJ), 9 est une 
variable auxiliaire liée aux trois premières par la relation. 

Conclusion, — Maxwell est arrivé au résultat suivant : pour 
que l'équilibre soit stable, on doit avoir : P<'n — . 

* * ^ 000 

Ce résultat est vrai pour un amas de poussières cosmiques 
comme pour un anneau liquide ; mais il y a une condition de plus 
pour un anneau liquide : la pression à la surface doit être dirigée 
vers l'intérieur ; or, Laplace a trouvé pour l'anneau à section 

elliptique que la densité doit être supérieure à —, pour l'anneau 

extérieur et à 2 pour l'anneau intérieur. 

On peut alors se demander s'il y a d'autres figures d'équilibre. 
Rappelons que l'on a trouvé qu'il fallait que 

(0^ < 2Ttp. 

Dans cette formule, (o doit être pris égal à la vitesse angulaire 
d'un satellite dont l'orbite coïnciderait avec l'anneau. On a donc, 
en désignant par M la masse de Saturne et par L le rayon de 
l'anneau. 

, M 

Si p, est la densité de Saturne, R son rayon, on a 



M = 4-^R', 



4 
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On trouve donc 



4 R' 






X©' 



Le &econd membre est supérieur à — tt- : si on rapproche ce ré- 
sultat du résultat précédent, on voit que l'anneau de Saturne ne 
peut être liquide. On a d'ailleurs vu plus haut que les observa- 
tions montraient qu'il était probablement composé d'une masse 
de petits corpuscules. 
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Page 3o, ligne i5 en remonlant. — Au lieu de — t»',*, il faut — to',?. 
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Page 37, ligne 10 en remontant. — Au lieu de -~^ — , il faut ., - . 

Page io3, ligne 3 en remontant. — Au lieu de — -, il faut r-- 

Page 104, lignes 6 et 10. — 

Au lieu de — 6 V-^ fx'/p'dz', il faut 3 ^[^ fx'y'p'dx. 



Page io5. — La dernière formule de la page doit être r 
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Page 107. — Les dernières formules de la page renferment quelques erreurs 
de signe faciles à apercevoir, d'où il suit pour la formule finale 






Page ia4, ligne 10 en remontant. — Le numérateur de la fraction sous le 
second radical doit être a^ — c^. 

Page ia8, ligne 12. — Lire S" — (Hp^ -(- K) S = o. La formule suivante 
doit être R"S — S"R = o, d'où en intégrant R'S — S'R = a. 

Page i34, lignes 12 et i5. — Lire S'R — R'S. 
Page i37, ligne 7. — Mettre le signe — devant /. 

Page 187, ligne 7. — Dans le dernier terme de la ligne, il faut x + x' au 
lieu de x — x'. 
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